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I.  INTRODUCTION 


Kinetic  energy  projectiles  and  shaped  charge  jet  warheads  require 
direct  hits  to  inflict  target  damage.  A  warhead  carrying  a  self-forg¬ 
ing  fragment  need  not  hit  the  target,  but  its  orientation  and  fuzing 
must  be  exactly  controlled  for  its  lethal  mechanism  to  hit  the  target. 
When  direct  hits  by  the  warhead  are  unlikely  (as  in  long  range  engage¬ 
ments  or  for  partially  concealed  targets)  or  not  required  (as  in  the 
case  of  a  collection  of  widely  distributed  targets)  fragmentation 
warheads  are  employed.  Since  direct  hits  are  not  contemplated,  the 
flight  and  fuzing  requirements  of  such  warheads  are  less  stringent. 
However,  the  design  of  the  warhead  itself  requires  great  care  if  its 
effect  on  the  target  is  to  be  nearly  optimum.  The  best  design  for  a 
fragmentation  warhead  obviously  depends  on  a  knowledge  of  target 
vulnerability  as  well  as  on  the  engagement  conditions  and  the  various 
options  that  are  available  regarding  fragment  number,  size  and  launch 
characteristics;  The  problem  of  designing  such  a  warhead  for  optimum 
effect  is  rather  complicated  and  is  made  more  difficult  by  the  uncer¬ 
tainty  of  our  knowledge  about  target  vulnerability  and  fragment 
characteristics.  If  each  warhead  is  required  to  address  more  than 
one  type  of  target  (as  is  often  the  case),  the  problem  is  still  more 
difficult.  In  spite  of  these  difficulties,  the  problem  of  optimizing 
or  at  least  of  estimating  the  effectiveness  of  a  design  has  received 
considerable  attention  over  the  last  half  century.  However,  for  a 
variety  of  reasons,  the  trend  has  been  toward  the  inclusion  of  more 
and  more  detail  in  model  descriptions  of  effectiveness  and  away  from 
design  optimization. 

In  current  design  practice  for  fragmentation  warheads,  the  focus 
is  usually  on  the  delivery  system.  Of  course,  this  system  is  impor¬ 
tant  and  inevitably  places  weight  and/or  volume  constraints  on  the 
warhead  designer.  However,  there  is  a  tendency  to  incorporate  tradi¬ 
tional  warhead  designs  into  new  delivery  systems,  with  warhead  effect- 
tiveness  studies  performed  only  after  the  overall  design  has  been  fixed. 
This  tendency  is  understandable  but  regrettable.  It  seems  obvious 
that  the  goal  of  optimal  warhead  effectiveness  should  drive  the  overall 
design  of  a  system,  especially  when  direct  hits  are  not  contemplated. 
The  warhead  designer  should  in  fact  propose  configurations  of  metal 
and  explosive  aimed  at  maximizing  target  kill  probability  even  if 
they  tax  the  ingenuity  of  his  colleagues  to  deliver  them  to  the  target. 
In  this  way  significant  improvements  might  result  from  compromises 
between  an  optimum  warhead  and  an  optimum  delivery  system.  However, 
our  methods  of  describing  warhead  effectiveness  need  to  be  simplified 
first  before  attention  can  be  re-directed  toward  those  factors  which 
most  influence  effectiveness.  The  eventual  object  of  this  report 
series  is  to  employ  simplified  descriptions  as  an  aid  to  design 
optimization.  The  object  of  this  report  is  to  review  the  development 
of  our  present  methods  as  a  prelude  to  suggesting  improvements. 
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II.  EFFECTIVENESS  STUDIES  UNTIL  THE  END  OF  WWII 
A.  Methodology  Prior  to  WWII 

There  seems  to  be  little  extant  literature  on  this  subject  written 
prior  to  the  second  world  war.  Of  course,  devices  which  explosively 
launched  pieces  of  metal  were  in  use,  but  early  treatises  on  ballistics 
emphasized  launch  and  flight  techniques.  Relatively  little  attention 
was  paid  to  measures  of  effectiveness,  either  absolute  or  relative. 

This  was  reflected  in  the  organizational  structure  of  the  Ballistics 
Research  Laboratory  when  it  was  established  in  1938.  There  were 
Interior  and  Exterior  Ballistics  Sections,  but  no  Terminal  Ballistics 
Section  until  1943.  Instead,  a  small  Effect  of  Fire  Unit  within  the 
Interior  Ballistics  Section  conducted  fragmentation  and  penetration 
studies.  Eventually  the  Terminal  Ballistics  Section  was  entitled  a 
Laboratory  and  gave  birth  to  the  Vulnerability/ Lethality  Division. 

The  need  for  an  overview  of  all  aspects  of  cost/effectiveness  led  to 
the  establishment  of  a  Weapons  Systems  Laboratory  within  BRL  and 
eventually  to  the  Army  Materiel  Systems  Analysis  Activity1. 

The  earliest  paper  on  the  subject  of  fragmentation  effectiveness 
appears  to  be  Kent's  report  in  19332.  Kent  took  as  a  measure  of 
effectiveness  the  number  of  casualties  produced  (personnel  put  out 
of  action).  He  commented  on  the  desirability  of  measuring  effectiveness 
by  the  extent  of  damage  produced,  especially  for  materiel  targets  like 
aircraft  where  the  number  of  man-hours  needed  to  repair  the  damage 
might  be  a  useful  criterion.  However,  for  footsoldiers  he  choose  a 
simple  threshold  criterion  to  distinguish  between  a  casualty  and  a 
survivor.  He  directed  his  attention  to  the  simple  case  of  a  vertical 
cylindrical  shell  bursting  at  ground  level  in  an  effectively  infinite 
field  of  uniformly  distributed,  unshielded,  standing  personnel  targets. 
He  wrote  an  expression  for  the. number  of  casualties,  C,  namely, 

00 

c  =  Ta  t  r**  ta  f  (Ap/r2)  (2ir  rdr)  (1) 

rl 

where  T  is  the  constant  number  of  targets  per  unit  area  of  ground. 

A  2 

The  radius,  r^,  defines  an  area  (irrp  on  the  ground,  centered  at  the 

burst  point,  within  which  a  target  will  receive  at  least  one  lethal  hit. 
All  targets  in  this  area  will  be  casualties,  but  there  is  a  loss  in 
efficiency  since  multiple  lethal  hits  are  wasteful.  Consequently, 
the  distance  r^  was  called  the  radius  of  overhitting. 


1  Bailie ticians  in  War  and  Peace.  A  History  of  the  United  States 
Army  Ballistic  Research  LaboratorieSj  Vol.  I,  1914-1956  ,  Edited 
by  John  Schmidt3  Aberdeen  Proving  Ground ,  MD>  pp  123  21a  47  and  98-102. 

2R.  H.  Kentj  "Considerations  on  the  Effect  of  the  Size  of  the  Projectile 
on  the  Efficiency  of  its  Fragmentation ",  BRL  RX-583  Feb  1953. 
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A  target  standing  a  distance  r  from  the  burst  point  was  replaced  by 
an  average  cross-sectional  area  normal  to  r  which  we  call  the  presented 
area,  A^.  Kent  did  not  discuss  this  area  in  detail  and  merely  remarked 

.  2 

that  each  target  subtends  a  solid  angle  (A  /r  )  at  the  burst  point. 

This  is  equivalent  to  replacing  each  target  by  a  spherical  cap  centered 
at  the  burst  point  of  area 


2it  r  9 


-•  n 


r2  sin  9  d  0  d  (Ji  =  2irr2  (1-cos  9)  =  4irr2  sin2  (9/2) 


=  8r 


(2) 


where  9  is  the  co-latitude  measured  from  the  pole  r  and  ij>  is  the 
equatorial  angle.  Thus,  the  solid  angle  subtended  by  this  area  is 

Si  =  Ap/r2  =  4  it  sin  2  (6/2). 

2 

If  Ap  <<  r  ,  the  spherical  cap  can  be  thought  of  as  practically  flat, 

leading  to  the  "poker  chip"  approximation  for  the  target  area.  No 
distinction  was  made  between  the  presented  area  and  that  part  of  this 
area  which  might  be  vulnerable,  since  Kent  was  able  to  eliminate  Ap 

from  Equation  Cl)  by  using  his  definition  of  r^.  The  symbol  stands 

for  the  number  of  lethal  fragments  emitted  per  unit  solid  angle,  so 
that  Si  is  the  number  of  lethal  hits  expected  on  Ap,  In  other 

words,  the  number  of  lethal  hits  on  a  single  target  is  expected  to  be 


NL  -  \  “  ■  NI!  Vr2  (3> 

which  is  unity  for  r<r^,  by  hypothesis  and  less  than  unity  for  r>r^. 

The  number  of  targets  in  a  ring  of  radius  r  and  width  dr  is  T^ 

(2irrdr)  and  N  times  this  number  is  the  number  of  casualties  expected 

L  2 
in  this  differential  area.  If  we  integrate  and  add  T^(ur^  )  we 

obtain  Equation  (1). 

In  order  to  carry  out  the  integration  we  must  make  some  assump¬ 
tion  about  the  r-dependence  of  N  .  This  number  is  expected  to  be  a 

do 

decreasing  function  of  r  since  it  is  observed  that  for  two  fragments 
of  the  same  shape  launched  at  the  same  speed,  a  less  massive  one  will 
lose  speed  more  rapidly  as  r  increases  because  of  air  drag.  If  both 
are  above  some  lethality  threshold  (mass  and  speed  combination)  initially, 
the  less  massive  fragment  will  fall  below  this  threshold  at  a  smaller 
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value  of  r.  Kent  considered  only  the  case  in  which  all  fragments  are 

identical  (same  mass,  shape  and  launch  speed)  so  that  all  are  lethal 
out  to  a  radius  and  ineffective  for  r>r2»  To  represent  this 

situation  he  let  N  be  equal  to  its  value  at  r=r. ,,  namely,  letting 

Npl  in  Equation  (3),  N^r^/Ap  for  rj  <  r  <  r2  with  N^=0  for  r>r2< 

This  eliminates  Ap  from  Equation  (1),  which  becomes 

C  -  Ta  tt  t*  [l+2  In  (r^rp]  .  (4) 


Kent  used  Equation  (4)  to  compare  the  effectiveness  of  two  shells 
with  the  same  geometry,  scaled  so  that  the  maximum  range  of  the 
larger  is  greater  than  that  of  the  smaller,  that  is,  r2L  =  rL  >  rs 

=  r2y.  The  value  of  r^  was  taken  to  be  the  same  for  both  since  the 

fragment  density  is  the  same.  Then  a  measure  of  their  relative 
effectiveness  is  the  ratio  of  Equation  (4)  for  each,  namely, 

CL/CS  =  [l*2  In  (rL/rj)J  /  [l*2  In  t^/rp] 

-  1  *  (2  In  CrL/rs)J  /  [  1*2  In  C^g/r^J 

<  1  +  2  In  (rL/rs)  (5) 


for  r„>r1 .  For  example,  if  the  larger  shell  has  twice  the  diameter 

and  twice  the  height  of  the  smaller,  it  will  weigh  about  eight  times 
as  much.  For  this  case  Kent  estimated  rr/rc  =  3.67,  provided  both 

launch  speeds  =  3,500  ft/sec.  Then  C ,/Cc  <  3.6  in  Equation  (5), 

L  j 

and  on  a  weight  basis  (C^/Wp  /  (Cg/Wp  <  3.6/8  =  0.45  where  W  is  the 


weight 


weight  for  rL  >  rg 


so  the  larger  shell  is  less  than  half  as  efficient  per  unit 
> 


and  CL/CS  = 


rr 

^rll/rlS^ 


If  r..  > 
(3. 67) 2  = 


rL  >  rS*  tlien  =  0  in  Equation  (1) 
13.47  and  on  a  weight  basis  the 


larger  is  more  efficient  by  a  factor  13.47/8  =  1.68.  The  larger  will 
be  less  efficient  only  if  C  <  8C  or  r1T  <  */iTriy, 

Li  b  1L  lb 


In  the  case  just  considered,  the  assumption  was  made  that  x^  > 


'  2S 


since  the  fragments  of  the  larger  shell  are  more  massive  but  not  more 
numerous.  Kent  further  considered  the  case  in  which  r,.  =  roc  =  rc . 

that  is,  the  case  in  which  both  shells  emit  fragments  of  the  same  mass 
so  that  the  larger  emits  more  fragments,  Kent  explicitly  assumed  that 
the  larger  would  emit  eight  times  as  many  fragments  and  implicitly 
assumed  as  before  that  the  launch  speeds  would  be  the  same.  Consequently 
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will  be  eight  times  bigger  for  the  larger  shell  and  its  radius  of 

overhitting  will  be  r1L  =  /8~  r  g  from  Equation  (3)  with  =  1 . 

Now  in  Equation  (4)  r  =  rc  for  both  shells  while  r, ,  =  /8  t.„  -  /8  r.  so 

L  h  1L  1  o  1 

CL/CS  =  s[u2  ln{rE  /  (/Tr,)}]/  [l*2  In  (r^)] 

=  8  1  -  In  8/ |  1+2  In  O^/r^j  .  (6) 

The  relative  efficiency  per  unit  weight  is  one  eighth  of  this  and  is 
necessarily  less  than  unity  if  r  >  r, .  If  r_  <  r,,  Nn  =  0  and 

C  /Cc  =  1  instead  of  8  in  Equation  (1),  since  hitting  each  target  more 
L  j 

often  with  lethal  fragments  does  not  increase  the  number  of  casualties. 

2 

For  each  shell  the  number  of  casualties  will  be  (nr^) .  Thus  on 
a  weight  basis  the  relative  efficiency  is  1/8. 

In  general  then,  Kent  concluded  that  smaller  shells  are  more 
efficient  on  a.  weight  basis,  provided  the  fragments  are  not  so  small 
that  the  effective  length  of  the  fragment  trajectory  is  less  than 
the  radius  of  overhitting.  Kent  realized  he  was  considering  a  highly 
idealized  case  and  did  not  concern  himself  with  the  absolute  accuracy 
of  his  estimates.  Implicitly  he  assumed  that  his  errors  were  about 
the  same  for  homologous  shells  so  that  significant  statements  could 
be  made  about  relative  efficiencies,  at  least  in  the  size  range 
considered  in  practice. 

In  1937  Tolch3  discussed  the  data  available  at  Aberdeen  Proving 
Ground  in  terms  used  at  the  time,  namely,  the  effective  area  or 
effective  radius  of  a  burst.  The  effective  area  was  defined  as  an 
area  within  which  each  target  on  average  received  at  least  one  lethal 
hit  while  the  effective  radius  was  that  of  a  circle  of  equivalent 
area.,  For  a  vertical  cylinder  burst  at  ground  level,  the  effective 
radius  is  the  same  as  Kent's  radius  of  overhitting  in  Equation  (1)  above. 
More  generally,  this  area  will  be  more  or  less  elliptical  when  the 
axis  of  the  shell  is  inclined  toward  the  horizontal  because  of  the 
dominance  of  the  side  spray  over  nose  and  base  sprays.  This  remains 
true  when  account  is  taken  of  the  remaining  velocity  of  the  shell  when 
it  strikes  the  ground.  The  main  effect  of  this  added  velocity  component 
is  to  throw  the  side  spray  forward  and  render  the  base  spray  even  less 
effective.  Tolch  remarked  that  the  attitude  of  the  target  (horizontal 
or  vertical,  face  or  side  on)  will  influence  effective  areas,  but  he 
contented  himself  with  an  average  value  for  Ap.  For  various  munitions 

%.  A.  Tolch ,  "Effective  Area  of  Burst  and  Fragmentation  Efficiency  of 
Hand  Grenades ,  arid  37  and  81  MM  Projectiles ",  BEL  E903  22  Dec  1937. 
(AD#493519) 
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he  estimated  effective  areas  from  the  number  of  perforations  through 
wood  panels  in  "silhouette  fields".  Finally,  he  compared  the  number 
of  shells  required  to  neutralize  a  target  field  as  predicted  by  the 
effective  area  method  with  that  predicted  by  Kent's  method,  using 
C/(TaAp)  in  Equation  (1)  above  as  a  measure  of  efficiency  with  Ap 

chosen  to  be  one  square  yard.  Multiplying  this  measure  by  other 
values  of  or  Ap  would  then  give  values  for  particular  target 

configurations.  Obviously  the  effective  area  method  predicts  that  more 
rounds  will  be  required  than  Kent's  method  predicts  since  it  considers 
only  the  region  of  100%  effectiveness,  the  first  term  in  Equation  (1). 

A  year  later  Tolch  published  a  report4  giving  details  of  more 
elaborate  fragmentation  and  penetration  experiments.  He  used 
horizontal  gun  firings  of  75MM  shells  bursting  in  the  center  of 
semi-circular  arrays  of  one-inch  thick  wood  panels  placed  at  various 
radial  distances  to  measure  the  number  of  fragments  per  unit  solid 
angle  as  a  function  of  co- latitude,  6  (measured  from  the  nose)  and 
distance  from  the  burst  point,  r.  This  quantity,  (r.,  0) ,  was 

studied  as  a  function  of  remaining  velocity  of  the  shell  at  the 
moment  of  burst.  Fragments  were  classified  into  nose  spray  (fuze 
pieces),  side  spray  (about  30  degree  angular  width)  and  base  spray  as 
well  as  into  perforating,  penetrating  and  denting  (<  1/16  inch  pene¬ 
tration)  categories.  From  these  panel  tests  and  the  axial  symmetry 
of  the  shell  it  was  estimated  that  about  5,000  fragments  were  produced. 
These  could  be  roughly  divided  into  700  perforating,  900  penetrating 
and  3,400  denting  fragments.  As  the  remaining  velocity  increased, 
the  nose  spray  panel  effects  became  more  pronounced  (more  perforations), 
the  center  of  the  side  spray  moved  forward  from  6  =  95°  (static  firing) 
to  0=60°  (for  a  remaining  velocity  of  1,100  f/s),  and  base  spray  effects 
decreased,  since  the  remaining  velocity  reduced  the  net  fragment  launch 
velocity  in  this  direction.  Complementary  sand  pit  tests,  in  which 
fragments  above  a  certain  size  were  collected  and  graded  by  sieves, 
produced  about  780  fragments  which  accounted  for  95%  of  the  metal 
weight,  indicating  that  most  of  the  metal  ends  up  in  the  perforating 
category  if  one- inch  thick  spruce  boards  are  used  as  witness  panels. 
Quarter- inch  thick  steel  plates  were  found  to  be  much  tougher  targets 
allowing  only  1/6  the  number  of  perforations  at  15  feet  and  only  one 
perforation  at  36  feet,  Tolch  also  published  reports  on  the  effective¬ 
ness  of  shrapnel  munitions5*6. 


4 N .  A.  Tolch3  n Fragmentation  Effects  of  the  75MM  H.  E.  Shell  T3  (M48) 
as  Determined  by  Panel  and  Pit  Fragmentation  Tests ",  B,RL#R126^^2  Dec  38. 

5N.  A.  Tolch ,  "Effective  Heights  of  Burst  and  Patterns  of  75MM  M37 
Schrapnel  fired  from  the  75MM  Pack  Howitzer  Ml ",  BRL  R100 ,  l  Apr  38. 


SN.  A.  Tolch ,  "Patterns  of  Pack  Howitzer  Schrapnel  at  the  Optimum 
Heights  of  Burst ",  BRL  R102s  16  May  38.  (AD#492929) 
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Two  years  later  Tolch  published  another  report7  in  which  he  noted 
a  number  of  differences  between  the  test  results  of  his  1938  report, 
obtained  under  symmetrical  conditions  with  particular  efforts  to  avoid 
the  influence  of  the  ground,  and  the  results  which  obtain  under  field 
conditions  where  asymmetries  exist  and  the  ground  has  a  large  influence. 
For  example,  when  a  shell  bursts  close  to  the  ground,  many  of  the 
fragments  are  deflected  upward  by  the  ground  and  do  not  follow  straight 
trajectories  even  over  relatively  short  distances  or  they  are  buried 
in  the  ground.  In  addition,  fragments  which  start  on  straight 
trajectories  are  in  flight  because  of  gravity,  pass  over  targets  close 
by  and  hit  targets  farther  away.  Lewy®  later  discussed  such  fragment 
trajectories,  using  the  traditional  methods  of  exterior  ballistics. 

He  pointed  out  that  the  supersonic  portion  of  a  trajectory  is  essen¬ 
tially  straight  since  the  gravitational  force  is  negligible  compared 

2 

to  the  drag  force  which  depends  on  v  .  He  also  gave  methods  for 
calculating  the  point  of  impact  and  terminal  speed  for  all  fragments 
projected  from  a  burst  a  height  h  above  the  ground.  However,  he 
pointed  out  that  fragments  initially  projected  upward  will  reach  a 
terminal  speed  in  air  on  the  downward  part  of  their  flight  which  is 
low  enough  to  render  them  non-lethal  for  common  fragment  masses.  The 
situation  may  differ  for  very  high  altitude  bursts  and  targets. 

Tolch  also  mentioned  the  natural  unevenness  of.  terrain  which  might 
intercept  some  parts  of  a  fragment  pattern  as  well  as  man-made  shield¬ 
ing  like  that  afforded  by  fox-holes.  Such  shielding  has  been  included 
in  later  models.  In  general,  Tolch  pointed  out  that  model  estimates 
of  fragmentation  effectiveness  based  on  idealized  panel  tests  will  not 
accurately  represent  effectiveness  under  service  conditions.  However, 
valuable  insight  into  at  least  relative  efficiencies  of  various  shell 
designs  under  different  engagement  conditions  should  be  obtainable 
from  such  models. 

In  discussing  angle  of  fall,  Tolch  distinguished  between  the  shell 
axis  and  the  tangent  to  the  trajectory  at  the  moment  of  burst.  However, 
he  made  no  use  of  this  distinction  and  assumed  that  they  were  co- linear 
for  all  practical  purposes.  By  using  an  x,  y,  z  =  h  coordinate  system 
with  origin  at  the  burst  point,  Tolch  wrote  an  expression  for  the 
slant  distance  from  the  burst  point  to  a  point  on  the  ground,  namely. 


7N.  A.  Tolch y  " Auxiliary  Curves  for  use  in  the  Computation  of  Fragment 
Density  as  a  Function  of  the  Angle  of  Fall  and  Height  of  Burst" , 

BRL  R178 ,  29  Feb  40.  (AD#491803) 

Lewy,  "Asymptotic  Integration  of  Fragment  Trajectories", 

BRL  RS59,  13  Jul  45. 
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where  h  is  the  height  of  the  burst  above  the  ground.  Unlike  Tolch,  we 
will  take  the  z  axis  postive  in  the  upward  direction  for  convenience  in 
later  comparisons,  so  that  the  ground  point  directly  below  the  burst 
point  is  z  =  -h.  Tolch  defined  the  angle  of  fall,  w,  in  the  x,  h  plane 
between  the  shell  axis  and  the  horizontal  x-axis,  so  that  the  projec¬ 
tion  of  R  on  the  shell  axis  can  be  written  in  terms  of  the  angle  6 
between  R  and  the  shell  axis  as  follows 

R  cos  0  =  x  cos  <ii  -h  cos  (tt/2-w)  =  x  cos  w  -h  sin  w,  (8) 

by  using  the  relation  for  the  angle  between  two  lines  in  terms  of 
their  direction  cosines.  Given  h.-and.^we  may  compute  R  and  0  for 
any  x,  y  pair.  If  we  plot  curves  of  constant  R,  we  obtain  circles  in 
the  ground  plane.  Tolch  also  plotted  curves  of  constant  0,  but  did 
not  point  out  their  general  character,  although  this  is  easily  done. 

If  we  eliminate  R  between  Equations  (7)  and  (8)  we  obtain  the 
quadratic  form 

Kcos  ui\2  ,1  2  „  /cos  i^A/sin  ur\,  2  ,  2["/sin  ur\2  ,1  „ 

-  ‘J*  -  2  *  -y  *  h  |(«rsj  ~T  0  t9) 

which  has  the  discriminant 

5  =  CO 

If  6  <  0,  the  curve  is  an  ellipse,  if  <5  =  0,  it  is  a  parabola  and  if 
6  >  o,  it  is  a  hyperbola. 

If  the  target  area  vector  A  makes  an  angle  (J>  with  R,  then  the 
area  normal  to  R  is  Ap  =  A  Cos  <(>  and  the  total  number  of  lethal  hits 

expected  on  a  target  is 

nl  ■  [Na  tR-  9>  /  r2  ]  V  til) 

where  (R,  0)  is  known  from  Tolch1 s  controlled  experiments,  choosing 

perforations  or  some  other  criterion  as  evidence  of  lethality.  Equa¬ 
tion  (11) is  the  three-dimensional  analog  of  Equation  (3)  and  includes 
a  dependence  on  angle  as  well. 
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In  1938  and  1939  Kent  published  four  reports  concerning  the 
effectiveness  of  fragmentation  shell  against  aircraft9-12.  In  these 
reports  he  considered  costs  and  interior  and  exterior  ballistic 
factors  as  well  as  terminal  effectiveness  in  comparing  proposed  anti¬ 
aircraft  artillery  shells  of  90mm,  105mm,  4.7  inch  and  6  inch.  He 
readily  admitted  that  the  weakest  point  in  his  analysis  was  the  nature 
of  the  assumptions  he  was  forced  to  make  about  the  effectiveness  of 
fragment  hits  on  the  target.  In  general,  he  pointed  out,  small 
fragments  could  be  lethal  to  personnel  and  fuel  cells,  although  large 
fragments  would  be  required  to  damage  engines  or  other  less  vulnerable 
structural  components.  As  a  plausible  compromise  measure  of  effective¬ 
ness  he  initially  used  the  average  of  two  numbers  for  given  shells, 
namely,  their  fragment  number  ratio  (as  determined  by  pit  tests)  and 
their  weight  ratio  (since  larger  shells  should  have  more  larger  frag¬ 
ments)  .  However,  he  soon  substituted  another  criterion,  namely,  the 
product  of  the  number  of  fragments  emitted  by  a  shell  and  their  aver¬ 
age  effective  range.  He  argued  that  this  would  be  a  reasonable 
measure  of  shell  effectiveness  by  considering  the  target  to  be  composed 
of  a  large  number,  n,  of  small  (differential)  areas  of  mean  size,  a, 
such  that  the  total  presented  area  of  the  aircraft  could  be  written  as 
A  =  na.  ,  If  N  is  the  fractional  number  of  hits  on  area  a  such  that 
V  2  a 

N  =  (N„/R  )  a  <  <  1,  analogous  to  Equation  (11)  above,  then  the 
probability  of  not  hitting  the  area  a  is,  by  the  binomial  theorem 


1  -  -H.  *  ^  -  .  ‘(V*2)  (1 

since  n  >  >  1,  N  <  <  1  and  Ap  =  na.  Then  to  a  good  approximation, 

provided  the  target  is  in  the  fragment  spray,  the  chance  of  at  least 
one  lethal  hit  is 


P  =  1 


-(N!1/r2)  AP  -  ^Nfl/R2j 


(13) 


so  that  Equation  (13)  reduces  to  Equation  (3)  for  Ap  <  <  R 


3R.  H.  Kent,  "A  Comparison  of  Antiaircraft  Guns  of  Various  Calibers ", 

BRL  R125,  1  Dec  38.  (AD#492691L) 

10i?.  H.  Kent ,  The  Probability  of  Hitting  an  Airplane  as  Dependent 

upon  Errors  in  the  Height  Finder  and  Director ",  BRL  R127,  14  Dec  3S.(AD#701865) 

ni?.  H.  Kent,  "The  Probability  of  Hitting  Various  Parts  of  an  Airplane 
as  Dependent  on  the  Fragmentation  Characteristics  of  the  Projectile", 

BRL  R132,  6  Mar  39.  (AD#491787) 

12fl.  H.  Kent,  "The  Probability  of  Disabling  an  Airplane",  BRL  R15S, 

5  Sep  29.  (AD#491799) 
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This  argument  of  Kent's  certainly  captures  the  spirit  of  the 
argument  generally  presented  today  for  the  Poisson  distribution  as  an 
approximation  to  the  binomial  distribution.  If  we  consider  a  (space 
or  time)  interval  A  divided  into  n  parts  of  size  a  =  A/n  with 

equal  to  the  constant  average  rate  of  occurrence  of  a  random  event  in 
the  interval  A,  then  the  probability  of  one  event  occurring  in  the 
sub-interval  A/n  is  N^A/n  =  p,  and  the  probability  of  no  events 

occurring  is  1  -  NAA/n.  If  these  events  are  randomly  distributed, 

the  probability  of  finding  exactly  s  of  them  in  n  intervals  (or  trials) 
is  given  by  the  binomial  distribution 


8  (S;  n)  =  si (n-s)!  (NAA/n) s  (1  -  NAA/n)n‘s  (14) 

for  s  =  0,  1,  2...  n  and  g  =  0  for  all  other  s.  In  the  limit  as  n  °» 
and  p  ->■  0  but  np  =  NAA  <<n  remains  finite,  we  have 


Um  g  =  XMll  r  lim 

n-*»  s  *  n-H* 


n! 

(n-s)i  ns 


j  £lim  (1- NAA/n) 


(N^ 

s! 


-N.A 
e  A 


=  P(S) 


(15) 


which  is  the  Poisson  distribution.  The  third  limit  in  Equation  (15) 
is  obviously  unity  while  the  first  can  be  shown  to  be  unity  as 
follows.  Divide  the  numerator  and  denominator  by  (n-s) 1  and  obtain 


lim 

n-w*> 


n(n-l) (n-2) . . , (n-s+1) 

s 

n 


=  lim 
n-*30 


1(1-  -)  (1-  -).. 
v  n'  v  n 


(16) 


If  we  let  z  =  -  NAA/n,  the  second  limit  in  Equation  (15)  becomes  the 
C-NaA)  power  of  the  limit  which  defines  the  base  of  the  natural 
logarithm,  namely. 


lim  [(!♦,) 13,1 1  ("NAA>  -  e-NAA  ' 

Finally,  the  probability  that  one  or  more  events  will  occur  in  the 
interval  is 
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S  PCS)  =  2 


P(S)  =  Cu  PCS)  -  P(0) 

S=1  s=o 


l  q‘NaA 

=  1  -  e  A 


(18) 


since 


s=o 


PCS) 


s=o 


(N^Ay 
s ! 


•N  A  +N. A  -N.A 
A  A  A  . 

e  =  e  e  =1 


Equation  (18)  is  the  same  as  Equation . (13)  if  we  identify  =  N^/R 

and  A  =  Ap,  and  gives  the  probability  of  at  least  one  lethal  hit  if 

N  is  the  number  of  effective  or  lethal  fragments  per  unit  area  of 
A  2 

fragment  front,  such  that  =  N^A  =  Ap/R  is  the  expected  number 

of  lethal  hits  and  is  much  smaller  than  the  number  of  trials  (total 
fragment  number) . 


(19) 

2 


? 

Kent  went  on  to  let  GR  dR  represent  the  probability  that  the 
aircraft  lies  in  an  elementary  volume  of  the  fragment  spray.  Here  G 
is  a  function  expressing  the  aiming  errors  of  a  particular  gun  versus 
aircraft  encounter.  Then  the  probability  that  a  shell  burst  will  be 
effective  can  be  written  as 


P 


(1  -  e 


-(N, 


ft 


/R2)Ap 


)  GR  d  R 


(20) 


where  R^  is  the  effective  range  of  a  fragment  (a  single  value  if  all 

fragments  are  alike  or  an  average  value  in  the  case  of  natural 

2 

fragmentation).  When  Ap  <  <  R  ,  then  Equation  (20)  becomes 


P  =  G  Ap  (Nn  R2)  t21) 

provided  G  Ap  is  independent  of  R.  For  comparable  aiming  systems 

against  the  same  target,  G  Ap  will  be  the  same,  so  that  (Nfi  R2)  or 

even  ft  R2  =  (NR2)  can  be  used  to  compare  the  relative  effectiveness 

of  different  size  shells.  Kent  expanded  the  exponential  in  Equation 
C20)  and  retained  the  first  few  terms  in  order  to  estimate  P  for  the 
shells  being  considered.  He  also  suggested  that  G  have  the  form  of  an 
error  function,  but  did  not  carry  out  any  calculations  employing  this 
suggestion.  Finally,  he  estimated  that  the  probability  of  a  direct 
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hit  on  the  aircraft  by  a  shell  was  as  great  as  the  probability  of 
killing  an  unprotected  essential  component  like  the  pilot.  In 
order  to  do  this  he  replaced  the  aircraft  by  two  intersecting  cylinders, 
one  representing  the  fuselage  and  the  other  representing  the  wings. 

He  concluded  that  prompt  point  detonating  fuzes  should  be  considered 
in  addition  to  time  fuzes,  especially  if  aiming  techniques  were  improved. 

In  the  fourth  report  of  this  series,  Kent  analyzed  the  results  of  a 
Navy  report13  which  gave  estimates  of  two  types  of  disability  or  "kill" 
based  on  experimental  firings.  Type  A  was  expected  to  cause  mission 
abandonment  and  an  immediate  forced  landing,  while  type  B  should  cause 
an  eventual  forced  landing  short  of  base  without  mission  abort.  These 
two  types  of  kill  were  given  for  hits  on  the  engine,  lubricating 
system,  fuel  system,  engine  controls,  engine  bearers,  pilot  and 
surface  controls.  Kent  suggested  that  these  disability  estimates  be 
treated  as  conditional  probabilities  of  kill  (type  A  or  B)  given  a 
hit  on  the  Ath  component,  namely,  *  If  each  component  had  a 

presented  area,  (Ap)^  ,  and  p^  is  tne  probability  of  hitting  a  unit 

area  of  the  aircraft,  then  (Ap) ^  p^^  is  the  probability  of  hitting 

the  Hth  component.  The  product  of  this  probability  with  (Pk/h) l 
then  gives  the  probability  of  disabling  the  airplane  by  a  hit  on  a 
given  essential  part.  Since  all  these  probabilities  are  small,  an 
approximation  for  the  net  probability  of  killing  an  aircraft  hit  by 
a  single  burst  is 


K/H 


=  ^  l  l  PH1  ”  ^  Avt 


■HI 


=  A 


v  p 


HI 


(22) 


Here  we  have  used  the  symbols  Ay^  and  Ay  to  represent  the  component 

and  total  vulnerable  areas  of  the  target.  Kent  did  not  use  the  symbols 
or  the  phrase  vulnerable  area,  but  he  computed  the  sums  for  both  kill 
types  using  the  Navy  data.  The  Navy  report  gave  the  number  of  frag¬ 
ments  hitting  each  part  and  the  number  of  estimated  A  and  B  kills. 

As  an  estimate  of  one  might  use  the  ratio  number  of  kills  N^ 

divided  by  number  of  hits  NH.  Instead,  Kent  considered  it  better 
practice  to  apply  Bayes*  rule  and  take  p^H  =  (NR  + 


1)/(Nh  +  2). 


Since  the  Navy  airplane  had  one  engine  and  one  pilot,  Kent  could 
not  directly  transfer  these  results  to  multi-engined  Army  bombers 
with  pilot  and  bombardier.  He  recommended  that  the  best  method  of 
obtaining  estimates  of  kill  probabilities  would  be  to  subject  such 
an  airplane  to  fragment  attack.  Because  the  Navy  tests  did  not 
involve  a  large  number  of  hits  on  vital  parts,  Kent  estimated  that 
the  values  he  calculated  could  easily  be  in  error  by  50%  or 


1  ^ Navy  Report  .25,  "A.  A.  Gunnery ,  Damage  Effect  as  Determined  by  tests 
against  XF  7B-1  Airplane  Navy". 
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more.  Since  no  tests  had  been  made  on  a  multi-engined  bomber,  P-k/h 
estimates  based  on  the  Navy  report  would  be  even  less  reliable.  Of 
course  the  (Pv/^H  values  used  by  Kent  were  for  an  average  over  all 
fragments  emitted  by  a  particular  shell  at  a  particular  location 
relative  to  the  target. 

About  the  same  time  in  England  Cunningham  and  coworkers llf  considered 
a  similar  problem  in  which  they  analyzed  the  effectiveness,  of  machine 
gun  fire  by  an  RAF  defender  attacking  a  twin-engined  bomber.  In  their 
report  they  gave  an  admirable  summary  description  of  the  goal  they 
were  trying  to  achieve: 

"The  volume  of  data  required  to  describe  with  complete  fidelity 
all  the  changing  conditions  which  occur  in  any  one  actual  combat  is 
clearly  too  large  to  admit  either  of  specification,  or  of  useful 
analysis.  This  difficulty  seems,  at  first  sight,  to  preclude  the 
possibility  of  a  mathematical  theory  of  combat  for  which  much  realism 
can  be  claimed.  It  has  nevertheless  been  proved  that  the  statistical 
results  of  a  large  population  of  similar  combats,  on  which  strategical 
interest  is  mainly  centered,  must  depend  only  on  mean  values  of  the 
varying  influences..." 

In  this  report  reference  is  made  to  a  previously  published 
"Mathematical  Theory  of  Air  Combat"  by  L.  B.  C.  Cunningham  without 
date  or  report  number.  This  report  is  now  unavailable,  but  modifica¬ 
tions  to  it  appear  in  Appendix  I  of  the  report  we  are  now  considering. 
This  appendix  speaks  of  "the  whole  potentially  vulnerable  area",  S, 
of  a  bomber,  and  the  set  of  fractions,  0  ,  contributed  by  thp 
&th  vital  component  to  this  area  as  determined  by  firing  trials 
such  that  the  actual  total  vulnerable  area,  \  ,  is 


The  pilot,  gunner  and  fuel  tanks  were  considered  "singly  -  vulnerable" 
components,  while  each  engine  consisted  of  a  "duplicated  pair,  whose 
contribution  to  the  bomber’s  vulnerability  may  be  analyzed  by  the 
standard  method".  Such  components  are  now  called  multiply  -  vulnerable. 
Unfortunately,  Cunningham's  other  reports  are  now  unavailable  so  we 
are  unable  to  explain  how  he  used  Equation  (23) .  However,  he  seems  to 
have  influenced  a  number  of  later  writers  as  we  shall  see. 


B.  C.  Cunningham  y  E.  C.  Comford ,  W.  Rudoe  and  J.  Knox ,  "An 
Analysis  of  the  Performance  of  a  Fixed-gun  Fighter  Armed  with  Guns 
of  Different  Calibres  in  Single  Home-defence  Combat  with  a  twin- 
engined  Bomber ",  Feb  40. 
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B.  THE  VARIABLE  TIME  PROXIMITY  FUSE 


In  1942,  the  U.  S.  Office  of  Scientific  Research  and  Development 
issued  a  significant  report, OSRD  738,  authored  by  a  group  consisting 
of  Professor  Garrett  Birkhoff,  Dr.  Ward  F.  Davidson,  Dr.  David  R. 

Inglis,  Professor  Marston  Morse,  Professor  J.  Von  Neumann  and 
Dr,  Warren  Weaver lb.  In  Appendix  I  of  the  present  report  brief  summaries 
are  given  of  the  background  of  some  of  these  men  and  some  others  who 
were  involved  in  this  type  of  work  in  England  and  the  United  States. 

The  OSRD  report  we  are  now  considering  compared  the  effectiveness  of 
time  and  proximity  fuzing  for  shells  fired  by  guns  on  ships  at  sea 
defending  themselves  against  high  level  bomber  attack.  It  is  worth 
quoting  part  of  the  introduction  to  this  report  since  what  is  said  is 
still  true  today  in  many  cases,  although  sometimes  not  so  clearly 
recognized: 

"It  was  well  understood  that  many  items  of  experimental  knowledge 
are  now  lacking]  and  that  it  would  therefore  be  necessary,  at  various 
points  in  the  discussion,  to  use  estimates  based  not  on  established 
fact,  but  based  rather  on  judgment  or  even  guesswork.  A  serious 
attempt  has  been  made  to  call  attention  to  all  such  soft  spots  in 
the  argument;  bracketing  figures  are  frequently  used  to  indicate  the 
possible  range  from  'average*  estimates  down  to  'low'  estimates  and 
up  to  'high'  estimates;  and  in  the  cases  of  major  importance,  there 
is  some  indication  of  the  sensitivity  with  which  the  final  results 
respond  to  changes  in  estimated  quantities". 

Proximity  fuzes  were  relatively  new  and  the  aim  was  to  find  the 
ratio  of  kill  probabilities  for  the  same  shell  with  two  different 
fuzes.  Even  large  absolute  errors  should  not  seriously  influence 
estimates  of  relative  effectiveness.  As  these  authors  remarked,  the 
advantage  of  the  proximity  fuze  was  found  not  to  be  affected  impor¬ 
tantly  by  assumptions  made  concerning  target  vulnerable  area,  initial 
fragment  velocity  or  the  nature  of  velocity  loss  due  to  air  drag. 

Their  general  method  was  to  multiply  the  probability  that  a  burst  occur 
at  a  certain  point  relative  to  the  target  by  the  conditional  probability 
that  a  burst  at  that  point  would  inflict  specified  damage  and  to  sum 
these  products  for  all  possible  locations  of  the  burst  to  obtain  the 
total  probability  that  a  given  round  inflict  the  specified  damage. 

The  ratio  of  these  probabilities  for  proximity  and  time  fuzes  was 
called  the  advantage  ratio. 

A  cartesian  coordinate  system  was  located  with  origin  at  the 
effectively  stationary  target  which  is  being  attacked  in  its  forward 
hemisphere  while  in  level  flight.  The  z  axis  was  chosen  anti  parallel 
to  the  shell  trajectory  at  burst  and  the  x-axis  was  chosen  horizontal. 

A  Gaussian  distribution  with  circular  transverse  (x,  y)  error  a  and 


15  G,  Birkhoff ,  V .  F.  Davidson ,  Dr.  R.  Inglis ,  M.  Morse,  J.  Von  Neumann 
and  W.  Weaver,  " The  Probability  of  Damage  to  Aircraft  through  Anti¬ 
aircraft  Fire.  A  Comparison  of  Fuzes  When  Used  against  High  Level 
Bombers  Attacking  a  Concentrated  Target ",  OSRD  738,  16  Jul  42. 
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range  error  a  was  used  to  describe  the  probability  that  a  burst  occur 
in  a  volume  near  the  target,  that  is,  for  a  time  fuze 


P..  = 


,,  ,3/2  2 
(,2tt)  a  a 


-fffa  [-  *  I 


,  2  2, 
(.*  x/J. 


dxdy  dz  (24) 


2  r  2  2  2 

or  simply  (2tt a  )  exp  I  -  (x  +y  )/(2o  )J  dxdy  for  a  proximity  fuze 

which  possessej^  burst  surface  described  by  z  -  z  (x,  y) .  If  R  = 

(x2  +  y2  +  z2)  and  ^  is  the  latitude  measured  from  the  x,  y  plane 

toward  the  z-axis  (ifi  =  ir/2  -  0  where  0  is  the  co-latitude),  then  the 
expected  number  of  lethal  hits  on  the  target  will  be  a  function 
N^(R,  40  because  of  symmetry  about  the  z-axis.  These  authors  used 


the  poisson  distribution  to  represent  the  conditional  probability  that 

at  least  one  hit  be  lethal  to  the  target,  that  is  D  =(l-e_NL),  as  in 
Hquation  (18)  above,  following  the  practice  of  Pearson  and  Welch  in 
a  British  report  now  unavailable.  They  also  used  the  usual  air-drag 
law 


m  dv/dt  =  m  v  dv/dR  =  -C  p  A  v2  (25) 

D  Si 

where  m  is  the  fragment  mass,  C  is  the  drag  coefficient,  p  is  the 

u  3. 

air  density  and  A  is  the  average  area  presented  by  the  fragment  in 
flight.  Following  Pearson  and  Bishop1”  they  represented  this  area  for 
a  randomly  rotating  fragment  by 


A  =  k  m2/3  (26) 

where  k  is  a  constant  shape  factor.  If  we  integrate  Equation  (25)  and 
use  Equation  (26),  we  obtain  for  straight  trajectories 


1/3 

v  =  v  e-CD  pa  a/”  =  v  e-CD  pa  kR/n,1/3  =  V'”  R/”  (27J 

0  o 

We  discuss  the  derivation  of  Equation  (26)  in  Appendix  II.  For  the 
number  of  lethal  hits  they  wrote 


lb  E.  S.  Pearson  and  D.  J,  Bishop ,  " The  Derivation  of  a  Retardation  Law 
for  Shell  Fragments" ,  A.  C.  247 8 j  5  Aug  42. 
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which  is  of  the  form  N  =  — —  used  in  Equations  (18),  (13)  and  (3) 

KM  R 

above  with  Na  =  N  w  and  A  =  AytR)  rather  than  Ap.  Here  K(t)i)  is  the 

ratio  of  the  number  of  fragments  per  unit  solid  angle  to  the  average 
number  per  unit  solid  angle  over  the  hemisphere  forward  of  the  shell’s 
equatorial  plane.  The  function  K(ip)  was  determined  for  2 .5 0  intervals  from 
experimental  shell  firings  and  depends  on  altitude  because  of  decreasing 
remaining  velocity  of  the  shell.  Here  N  is  the  total  number  of 
fragments  larger  than  some  minimum  size  (.04  oz)  and  Ay(R)  is  "a  sort 

of  vulnerable  area"  which  the  target  presents  to  the  burst  and  is  defined 
by  Equation  (28) .  In  order  to  estimate  this  area,  its  value  wa3  assumed 
to  be  divided  into  three  component  areas,  A c  =  Ay/2,  representing  the 

crew  (pilot  and  bombardier),  and  =  Ap  =  Av/4,  representing  the 

motors  and  fuel  cells  respectively.  As  they  remarked,  there  was  "some 
evidence  that  two  fragments  hitting  a  fuel  tank  at  the  same  point  (or 
near  together)  at  short  intervals  of  time  are  especially  effective  in 
starting  fires".  Threshold  lethality  criteria  were  chosen  for  the 
three  component  areas,  namely,  120  foot  -  pounds  of  fragment  energy 
for  A^,  the  ability  to  perforate  3/16  inch  of  mild  steel  for  A^  and 

the  ability  to  perforate  1/8  inch  of  mild  steel  for  Ap.  The  fragments 
were  divided  into  fairly  homogeneous  groups  with  the  l-th  group 
containing  n.  fragments  of  average  mass  m. ,  residual  speed  v  (m.,  R) 
from  Equatioft  (27)  and  vulnerable  area  A^fnu,  v(itk,  R)].  Here1 

A.  =  A  .  =  A  .  +  A„.  +  A„.  =  52  A.,  for  all  groups  on  the  assumption 

l  vi  ci  Mi  Fi  “  Ai 

that  they  all  have  some  sufficiently  high  launch  speed  and  mass.  However, 
as  R  increases,  one  or  more  of  the  component  vulnerable  areas  will 
drop  out  of  the  sum  for  (A^  =  A^  (R)  because  of  loss  in  velocity,  so 

that  a  plot  of  each  A^  versus  R  will  be  a  step  function  approaching 

zero.  Consequently,  the  total  vulnerable  area 


*.«  '  <E  \  V  7  (£  =  k  E  ",  z  t29) 

i  i  i  i 

will  also  be  a  step  function.  Figure  1  reproduces  their  figure  for 
effectiveness  versus  R  for  a  Navy  5  inch  shell  with  initial  fragment 
speed  vq  =  3500  ft/sec.  Near  the  target  for  example  the  contribution 

from  1/2  oz.  fragments  is  given  by  a  vertical  line  segment  between 
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Figure  1.  Fall-off  in  Effecti 
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ss  Subdivided  According  to  Fragment  Weights 
from  OSRD  738) 
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the  1/2  oz.  and  the  1  oz.  curves.  These  1/2  oz.  fragments  show  losses 
in  effectiveness  near  150  feet  and  300  feet.  Their  vertical  scale  is 
arbitrary,  but  if  they  had  plotted  Av(R)  /  Av(0)  the  curve  would  approach 

unity  as  R  -*■  0 .  One  also  sees  why  they  could  approximate  the  cumulative 
effectiveness  curve  by  an  exponential,  namely. 


\W  =  Av(0)  e"CR  (30) 

where  C  depends  on  p  ,  Cn  and  k  as  well  as  fragment  mass.  The  ratio 

Ay(R)  /  Av(0)  was  introduced  into  Equation  (28)  for  convenience,  since 

it  could  be  computed  for  any  Av(0),  assuming  the  same  division  into 

component  areas.  Then  the  value  of  Ay(o)  could  be  chosen  later  to 

fit  a  particular  size  target  of  the  same  type.  The  above  scheme  for 
estimating  Ay(R)  and  Nl(R,i|>)  is  straightforward  and  similar  schemes 

have  been  widely  used  since.  It  is  interesting  to  note  the  authors' 
remarks  on  this  scheme  since  they  are  often  true  even  today: 

"It  must  be  confessed  that  the  computation  of  N  (R,ip)  is  a  some- 

L 

what  shaky  affair.  But  it  is  comforting  to  note  that  the  final  results 
of  the  whole  study  are  rather  surprisingly  resistant  to  alterations  in 
the  method  of  computing  N T".*7 

If  we  recall  that  their  study  was  of  relative  calculated  kill  probabi¬ 
lities  rather  than  of  absolute  values,  then  this  result  is  not  too 
surprising.  However,  they  also  noted  close  agreement  between  the  kill 
probabilities,  p  , -for  time  fuzes  which  they  calculated  and  those 

calculated  by  the  British: 

"Taking  into  account  the  variety  and  complexity  of  the  different 
assumptions  and  the  different  methods  of  calculation,  the  above 
comparisons  of  p^  seem  almost  miraculously  close;  while  the  essential 

coincidence  of  the  ranges  of  value  for  the  advantage  ratios,  as  just 
stated,  is  also  reassuring". 18 

If  we  put  Equation  (30)  into  Equation  (28),  we  find 

^^v(0)  -PR  2  -PR  2 

NL  =  -  K0I0  e  LK/R^  ■=  K  -010 e  LK/R  (31) 


17 Reference  153  p.  B-7. 
lBReferenoe  15 ,  p.  46. 
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which  for  large  R  is  approximately  equal  to  the  damage  function,  that 

is,  D  =  1  -  e  NL  N.  .  It  is  clear  that  R2D  *“  R2N.  ->-0  as  R  -*■  0 

and  that  R^D  ■+  0  as  R  -*■  o .  For  a  proximity  fuze  we  can  express^^R,^) 

in  cylindrical  coordinates,  namely,  D(r,z)  where  r  =  (x  2  +  y2) . 

Their  calculations  showed  that  D(r,z)  has  its  principal  maximum  for 
z=z(r),  a  surface  of  revolution  which  is  the  optimum  burst  surface  for 
a  proximity  fuzed  shell.  From  their  tables  it  is  apparent  that  this 
surface  is  approximately  a  cone  with  the  target  near  the  apex.  Such 
a  burst  surface  would  place  the  target  in  the  center  of  the  side  spray. 
In  general  they  concluded  that  proximity  fuzes  were  anywhere  from 
twice  as  effective  to  an  order  of  magnitude  more  effective  than  time 
fuzes,  depending  on  circumstances.  The  advantage  ratios  were  found  by 
evaluating  numerically  or  by  power  series  integrals  over  the  product 
probability  distribution  of  D  and  the  Gaussian  error  function,  namely. 


P1  = 


1 

J/2 


(2i0 

4d 

for  the  time  fuze  and 
1 


fff-  z2' 

1 — JJh  2  02  — 2  d  c*.  y. 

do  L  °  a  ^ 


z)  dx  dy  dz  (32) 


P(<0  = 


2nff 


JJ°  2a2 


D  [x,  y,  z(x,  y)  ]  dx  dy 


(33) 


for  the  proximity  fuze.  Here  p^  or  p  (o)  are  the  expected  values  of 

the  probability  function  D  when  x  »  y  and  z  are  random  variables  with 
Gaussian  probability  distributions. 

It  is  interesting  to  note  that  the  concept  of  vulnerable  area 

presented  in  Equation  (29)  is  more  general  than  Kent's  given  by 

Equation  (22)  where  Kent's  A^  =  is  explicity  equal  to  £  tk 
A  . ,  that  is,  equal  to  the  average  effect  of  all  the  fragments  on  the  Jlth 

XI 

component  when  launched  from  a  particular  location.  Kent  estimated 
component  (PK/H)£  from  experiments  and  used  A  ^  to  estimate  hit 

probability.  The  association  of  A  ^  with  (P^yp)  £  to  form  A^  does 

not  eliminate  the  need  to  estimate  (Pk/h)  £  experimentally  and  certainly 
does  not  imply  that  the  conditional  kill  probability  of  the  Jlth  component 
(Pk/h^ji  =  Av)l/A  ^  can  be  reduced  merely  by  increasing  A^,  that  is, 

by  simply  placing  the  sensitive  components  in  a  bigger  box!  As  is  clear 
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from  the  discussions  above,  A  „  or  A  or  A  .  are  "sort  of  areas". 

v£  vi.i  vi 

They  are  perhaps  better  described  as. effectiveness  parameters  which  are 
defined  to  have  the  dimensions  of  an  area  by  the  role  they  play  in 
Equation  (28) .  We  can  just  as  well  speak  of  total  or  component 

as  estimated  from  experimental  firings  (or  in  some  other  way)  using 
a  threshold  or  some  more  sophisticated  lethality  criterion.  This  has 
the  advantage  of  avoiding  confusion  between  "vulnerable  areas"  and 
physical  (spatial)  areas.  A  similar  confusion  can  arise  concerning 
so-called  "lethal  areas"  as  we  shall  see. 


Shortly  after,  a  British  report  by  Kendall19  gave  a  similar  comparison 
of  time  and  proximity  fuze  effectiveness  versus  aircraft,  but  for 
rockets  instead  of  gun  projectiles.  Since  launch  stresses  on  rocket 
missiles  are  not  as  great  as  those  on  gun-launched  shells,  pre-formed 
fragments  can  be  used  instead  of  relying  on  natural  fragmentation.  This 
forces  the  designer  to  choose  one  or  more  fragment  sizes  and  shapes. 
Kendall  also  included  systematic  bias  in  the  aiming  error  and  introduced 
an  important  approximation  for  the  chance  of  at  least  one  lethal  hit. 

He  adopted  the  coordinate  systems  of  Equations  (32)  and  (33)  above 
and  assumed  that  a  proximity  fuze  would  burst  approximately  on  the  cone 
specified  by  the  co-latitude  0  =  0q  (z  =  R  cos  0) .  He  included  a 

transverse  bias  error  a  by  rotating  the  system  about  the  z  axis  until 
the  projection  of  the  burst  point  on  the  x,  y  plane  lay  on  the  x-axis 
with  mean  position  (-a,  0)  instead  of  (0,0)  as  for  no  bias.  Since 
x  =  R  sin  0q  cos  <t>  and  y  =  R  sin  0  sin  <J>,  then  Equation  (33)  becomes 


P(a,a)  = 


1  /R 


Hxp 


27TCJ  Jo  Jo  L  2o 


D  (R,0  )  sin  6  R  dR  d<(> 
o  '  o 
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2  /  -0  OT  02r. 
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D(R,0o)  IQ  (bR)  RdR 


(34) 


where  IQ  (bR)  =  J ^  e^R  cos  ^  d<|>  is  a  Bessel  function  with 

b  =  a  sin  B^/a2.  For  D  in  Equation  (33)  he  proposed  the  approximation 

D  =  1  -  .  e-r2^2R2>  =e-R2  sin  2eo''(2R2>.  (35) 


Kendall  called  the  adjustable  constant  R2  the  probability  radius  and 
pointed  out  that  for  N.  =  K  e~CT^T  ~  1 

L 


19Z?.  G.  Kendall ,  "The  Chances  of  Damage  to  Aircraft  from  A.  A.  Rockets 
Fitted  with  Time  or  Proximity  Fuzes'\  A.  C.  2435 j  28  Jul  42. 
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.632 


(36) 


e-2/^ 


so 


R-2  s  1.044  r  «  r 


(37) 


Thus  can  be  thought  of  as  defining  a  circle  in  the  x>  y,  plane 

centered  on  the  target  within  which  we  expect  at  least  one  lethal  hit. 
This  is  comparable  to  Kent's  radius  of  overhitting  in  Equation  (1) 
above,  so  that  irR^  is  a  100%  "lethal  area"  in  a.  probability  sense.  If 

we  use  Equation  (35)  in  Equation  (34)  with  no  bias  error  (a  =  0) ,  we 
p(o)  .jfjo  e  2  (7  *  3)  [l-e-  ?(4  *  £)  -o'} 


'38) 


2  2  2  2 

since  r  =  R  sin  6q  gives  r  dr  =  sin  6q  RdR .  Equation  (34)  with 

Equation  (35)  can  also  be  integrated  in  closed  form  for  a  =£  0,  provided 
we  neglect  the  integral  from  rQ  to  In  this  case  we  find 


p(a,o) 


(38a) 


o2+R22  ®  2(o2+R22) 

Kendall  did  not  choose  R-  on  the  basis  of  how  closely  Equation  (35) 

^  “cR/R  j 

could  match  the  original  function  1-exp  (-Ke  '  Jfor  given  K  and  c. 
Rather  he  chose  R»  on  the  basis  of  how  well  Equation  (38)  for  given  a 
and  r  could  be  matched  to  p(o)  computed  numerically  using  the  original 

function.  Kendall  also  considered  the  problem  of  optimizing  a  in 
Equation  (38a)  to  compensate  for  an  unavoidable  a,  as  well  as  the 
probability  of  damage  from  a  salvo  of  rockets  equipped  with  proximity 
fuzes.  His  treatment  of  the  subject  was  later  repeated  by  Carlton20 
and  his  approximation  for  D  in  Equation  (35)  above  came  to  be  called 
the  Carlton  approximation. 


For  a  time  -  fuzed  warhead  without  bias  errors  but  with  equal  range 
and  transverse  random  errors  (a  =0),  Kendall  used  symmetry  to 
integrate  over  $  and  assumed  thlt  taking  6  =  would  describe  the 

side  spray  well  enough,  so  Equation  (32)  above  becomes 

Z^A.  G,  Carlton,  "The  VT  Fuze  in  Anti-Aircraft  Gunfire "3  TM122,  Aug  453 
the  Johns  Hopkins  Applied  'Physics  Laboratory. 
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neglecting  the  integral  from  Rq  to  <*>  and  using  the  approximation 

D  ,  e-R2/(2Rj2)- 


(40) 


analogous  to  Equation  (35).  The  factor  ^  sin  0q  can  be  buried  in 

the  probability  radius  R^  which  defines  a  kind  of  "lethal  volume". 

Kendall  also  briefly  considered  the  cases  in  which  o ^  a  and  in  which 

non-zero  bias  occur.  In  a  later  report21  he  applied  these  ideas  to  the 
prediction  of  an  optimum  size  and  shape  fragment  in  a  pre- formed  fragment 
warhead  attacking  a  high  altitude  bomber.  He  concluded  that  the  frag¬ 
ment  should  be  square  on  the  two  faces  which  intially  form  part  of  the 
shell  wall  and  should  weigh  from  1/8  oz  to  3/8  oz,  assuming  no  change 
in  fragment  launch  speed  from  shells  in  service.  He  and  his  co-workers 
pointed  out  that  this  conclusion  was  quite  sensitive  to  the  threshold 
criterion  used  to  define  the  perforating  power  of  a  lethal  fragment 
and  recommended  that  much  more  experimental  information  was  needed  to 
relate  perforating  power  to  aircraft  damage  before  much  confidence 
could  be  placed  in  these  conclusions. 

A  little  after  Kendall’s  paper  appeared  in  1942  Professor  N.  F.  Mott 
published  a  report  on  this  subject22  in  which  he  discussed  methods  of 
choosing  the  optimum  size  fragment  to  be  used  in  a  controlled  fragmen¬ 
tation  warhead  attacking  aircraft.  For  example,  against  dive  bombers, 
damage  to  controls  or  one  engine  may  stop  the  attack  even  though  the  air¬ 
craft  can  return  to  base.  If  attrition  of  a  high  altitude  bomber  fleet 
is  desired,  then  damage  to  all  engines  or  hits  on  the  fuel  may  cause 
a  crash  short  of  base.  To  achieve  these  different  objectives,  different 
numbers  and  sizes  of  fragments  might  be  desired.  Since  adequate 
estimates  of  the  vulnerability  of  vital  components  to  fragment  attack 
were  required  but  not  available,  Mott  contented  himself  with  an  outline 
of  a  method  while  pointing  out  the  weakness  of  the  experimental  base  on 
which  it  rested.  The  proof  that  available  vulnerability  estimates 
based  on  controlled  firings  against  aircraft  were  inadequate  he  took 
to  be  the  fact  that  predictions  based  on  these  component  vulnerability 
estimates  did  not  agree  with  the  results  of  combat  data: 

Z1A.  G.  Walters,  D.G.  Kendall  and  L.  Rosenhead ,  "Controlled  Fragmentation 
of  the  3"  (Parallel  Sided )  Rocket  Shell >  Fitted  with  a  Time  Fuze ,  in  the 
H.  A,  A.  Roles  The  Optimum  Mass  and  Shape  of  the  Fragments  when  they 
are  all  of  Equal  Size",  A,  C.  4467,  26  Jul  43 . 

ZZN.  F.  Mott3  "Damage  to  Aircraft  by  H.  A.  A.  Fire ",  A.  C.  2480 ,  Aug  42. 
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"Estimates  of  the  vulnerability  of  two  engined  bombers  obtained 
by  adding  these  components  do  not  agree  with  results  obtained  in  combat. 
These  can  be  deduced  from  the  results  of  H,  A.  A.  shooting  in  this 
country,  where  casualties  can  be  compared  with  the  estimated  accuracy, 
and  from  a  comparison  between  the  number  of  hits  on  Bomber  Command 
aircraft  which  return  and  the  estimated  number  shot  down  by  (German) 

H.  A.  A.  In  both  cases  the  numbers  shot  down  are  at  least  three  times 
as  great  as  they  should  be  on  the  assumption  that  damage  to  controls, 
pilot  or  both  engines  is  normally  necessary  to  cause  a  crash". 

Mott  devoted  half  of  his  report  to  a  discussion  of  the  shortcomings  of 
estimated  component  vulnerable  areas  in  the. light  of  combat  data.  Of 
course  he  was  aware  of  the  limitations  of  combat  data,  but  he  also 
knew  the  limitations  of  predictive  models  based  on  inadequate  experiments. 
It  is  not  surprising  that  a  man  of  Mott’s  ability  (see  Appendix  I) 
showed  a  healthy  skepticism  toward  his  own  computational  estimates  and 
measured  the  degree  of  their  predictive  ability  by  comparisons  with 
real-world  combat  data. 

In  section  six  of  his  report,  Mott  assumed  that  a  shell  burst 
projected  a  uniform  random  distribution  of  N  identical  fragments  each 
travelling  the  same  straight  line  distance  Rg  before  it  became  ineffec¬ 
tive.  Each  aircraft  in  a  fleet  of  T  aircraft  was  looked  upon  as  a 
sphere  which  presents  a  vulnerable  area  A  and  the  poisson  probability 
-N  2 

(1-e  L)  with  Nl  =  A/R  was  used  as  the  probability  of  at  least  one 
lethal  hit.  This  estimate  for  N  is  more  like  Kent's  Equation  (3) 

Li 

above  since  it  does  not  include  the  angular  dependence  of  the  fragment 
spray  as  described  in  Equation  (28).  The  use  of  the  poisson 
approximation  he  attributed  to  Cunningham  without  citing  a  reference. 

Tlie  number  of  casualties  or  aircraft  which  receive  at  least  one  lethal 
hit  from  a  burst  somewhere  in  the  fleet  is 

c  =  Tv/0E(1'eNlj)  4 it  R2dR  =  4tt  Tv  £re3/3  -^E  ~Nn  A/r2  R2dR]  (41) 

where  T^  =  T  /  it  Rg3^  is  t^ie  number  of  targets  per  unit  volume  of 

air  space.  Equation  (41)  is  a  three  dimensional  analog  of  Equation  (1) 
above,  but  differs  from  it  in  that  the  use  of  the  poisson  expression 
for  at  least  one  lethal  hit  avoids  the  necessity  of  introducing  the 
idea  of  an  overkill  radius. 

For  vulnerable  components  Mott  used  the  symbol  A1  with 

^1-e  »  N*  =  aVr2  <  <  1 ,  so  that  the  fraction  of  target 

components  of  a  given  type  (say  an  engine)  which  receive  at  least  one 
lethal  hit  becomes 
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(42) 


Y  =  V 


(j- 


4irdR  = 


3  N„  A1/^2 


since  and  A1  were  not  taken  to  be  functions  of  R. 


In  the  case  of 


redundant  components  as  in  a  two-engine  bomber,  Mott  recommended  using 
‘  ‘  2 


M2  -  (n02- 


In  section  seven  of  his  report  Mott  discussed  the  question  of 
optimum  fragment  size.  He  adopted  Equation  (27)  above  and  solved  it 
for  R,  obtaining, 


R  =  (CD  pa  k)_1  m1/3  in  (vQ/v)  .  (43) 

Next  he  adopted  a  threshold  lethality  criterion  based  on  the  momentum 
needed  to  perforate  a  given  thickness,  T,  of  wood 

mv  =  GTA  =  GTk  m2/3  (44) 

where  G  is  a  constant  and  the  average  presented  area  was  adopted  from 

Pearson  as  in  Equation  (26)  above.  Mott  used  Equation  (44)  to  define 

a  minimum  mass,  m 

m’ 

m1/3v  =  m  1/3  v  =  GTk  (45) 

mo  v 

which  depends  on  the  launch  speed,  v  ,  as  well  as  the  target  thickness, 

T,  and  the  shape  factor,  k.  He  then  used  Equation  (45)  to  eliminate 
vq/v  in  Equation  (43)  to  obtain  an  expression  for  the  effective  range 

Rg  =  (CD  pa  k)'1  m1/3  in  (m/raj  1/3  (46) 

a  form  he  attributed  to  Welch  without  reference.  Finally  he  noted 
that  if  M  is  the  mass  of  the  shell,  then  the  number  of  fragments  will 
be  N  =  M/m,  while  the  chance  of  a  lethal  hit  on  a  target  component  will 
be  proportional  to  the  product  NR^  as  in  Equation  (21)  above, 

NR£  =  j  M  (CD  pa  k)'1  m~ 2/3  In  (m/mj  (47) 

from  Equation  (46) .  This  function  of  m  has  a  maximum  for 
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(48) 


m/m  =  e3/2  =  4.48 
m 

so  that  all  fragments  should  have  this  mass.  If  there  are  two  vulnerable 
components,  one  of  area  A^  vulnerable  to  fragments  of  minimum  mass  m^ 

and  a  second  of  area  A  ,  vulnerable  to  a  minimum  mass  m  >  m. , 

&  f  2  1 

then  the  chance  of  a  lethal  hit  will  be  proportional  to 


m-2/l  j"  ^  ln  +  ln  (in/ni2)J  (49) 

by  analogy  with  Equation  (47)  above.  This  expression  has  a  maximum 
with  respect  to  m  when 


m 


■[ 


m. 


m„ 


3/2  (A 


1  + 


1/A:  +  A2) 


(50) 


Finally,  Mott  returned  to  the  evaluation  of  Equation  (41)  in  the 
2 

case  where  A/R  is  not  negligible  compared  to  unity  and  plotted 

C/CTyA)  versus  (m/mm)  since  the  upper  limit  of  the  integral  in 

Equation  (41)  depends  on  this  ratio  by  Equation  (46) .  He  evaluated 
the  integrals  numerically  both  for  a  singly  vulnerable  component  using 


ing  (l-e  Ni) 


and  for  a  multiply  vulnerable  component,  using  yi-e  L/  "  as 
explained  above.  As  expected,  these  measures  of  effectiveness  exhibited 
maxima  as  a  function  of  (m/m  )  with  a  broader  maximum  near  10  to  15  for 

singly  vulnerable  component  and  with  a  narrower  lower  maximum  near  4  to 
6  for  two  redundant  components. 


In  an  appendix  Mott  briefly  considered  the  case  of  natural  shell 
fragmentation.  Adopting  a  suggestion  of  Welch  which  he  described  in 
more  detail  in  a  later  report2*,  Mott  wrote  for  the  number  of  fragments 
with  masses  between  m  and  m  +  dm 


dN  =  N(m)  dm  =  Nq  e“(m/ra0)1/3  d  J\ra/mo)1/3]  (51) 

where  Nq  and  mQ  are  chosen  to  fit  experimental  mass  distributions. 
Instead  of  using  the  value  in  Equation  (41) ,  he  used  an  average 
value,  R^,  computed  by  weighting 


23 n,  F.  Mott  and  E.  E.  Linfoot ,  " A  Theory  of  Fragmentation ">  A.  C.  3348 > 
11  Jan  43. 


Rg  according  to  the  distribution 


51 


given  in  Equation  (51) ,  namely. 


^  =  /  Re  d(N/No)  =  ±  (cD  pa  k/mm1/3)  7;  yd"  /)  e^dy  (52) 

1/3  1/3 

where  y  =  (m/m  )  =  (m/m  )  ,  if  we  design  our  naturally  fragmenting 

o  m 

shell  so  that  mQ  =  mm.  The  infinite  upper  limit  of  Equation  (52)  is 

an  adequate  approximation  when  almost  all  of  the  fragments  are  very 
small  compared  to  the  original  shell  mass.  Since  is  also  a  function 
of  (ra/m^) ,  Mott  could  have  plotted  C/ (T VA)  versus  (m/mm)  for  naturally 

fragmenting  shell  examples  also,  but  he  did  not. 

In  a  later  report21*  Mott  compared  the  effectiveness  of  natural  and 
controlled  fragmentation  and  concluded  that  "the  optimum  average  size 
with  natural  fragmentation  is  so  small  that  if  all  fragments  had  this 
size  they  would  scarcely  be  large  enough  to  damage  the  target".  In 
addition  he  remarked,  "the  chance  of  damage  is  very  insensitive  to 
mean  fragment  size,  and  cannot  be  improved  by  more  than  20  per  cent  by 
finer  uncontrolled  fragmentation".  In  his  calculations  he  assumed  an 
aircraft  to  have  a  vulnerable  area  of  ten  square  feet  with  resistance 
to  damage  equivalent  to  the  resistance  to  perforation  offered  by  two 
inches  of  wood.  As  before  he  used  an  average  over  all  aspects  of  the 
target.  As  we  can  see  from  Equation  (41)  above,  the  vulnerable  area 
was  taken  to  be  a  fixed  number  for  a  given  target  and  shell  which 
entered  into  the  estimation  of  the  number  of  lethal  hits  N  =  N  A/R  . 

L 

Of  course  N.  =  (R)_while  fragment  size  and  speed  entered  through  R^ 

in  Equation1, (46)  Lor  R  in  Equation  (52),  Again  his  treatment  of 

b 

vulnerable  area  is  more  like  Kent's  than  that  given  in  OSRD  738. 

Meanwhile,  at  the  Ordnance  Office  in  Washington, D.C .  similar  studies 
were  being  carried  out  by  Prof.  Marston  Morse  and  his  associates. 

Prof.  Morse  had  already  collaborated  in  OSRD  738  and  continued  his 
interest  as  revealed  by  a  series  of  reports  from  1943  to  1945.  For 
damage  to  materiel  targets  Morse  preferred  an  energy  threshold  criterion 

2  2 

mv  =  K  T  A  =*  constant  (53) 

instead  of  Equation  (44)  above.  Here  again  T  is  the  target  thickness 
(wood  or  steel  plate  equivalence) .  However,  Morse  did  not  simply  use 
the  average  area  presented  by  a  fragment  in  flight  as  the  striking 


21*iV.  F.  Mott ,  " Optimum  Fragment  Size  for  A.  A.  Shell ",  A.  C.  3366} 
15  Jan  43. 
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area.  Instead  he  preferred  the  maximum  area  or  some  quantity  related 
to  it,  at  least  when  deliberately  trying  to  overestimate  the  protection 
needed  by  an  aircraft  from  its  own  bomb-fragments  in  very  low-level 
runs25.  If  we  solve  Equation  (53)  for  v  (given  m)  then  we  have  an 
expression  for  the  minimum  striking  speed  needed  to  perforate  the 
target  thickness,  T,  namely,  the  limit  speed 

vL  =  K  (TA)1/2  m"1/2  (54) 

which,  as  Morse  points  out,  is  a  particular  case  of  the  De  Marre 
formula 


1  -1/2 

vL  =  K  (TA)  1  m  ' 


This  later  served  as  the  basis  for  the  formula 


C  a  B  y 

vL  =  10  (TA)  m  (sec  0) 


(55) 


(56) 


used  in  Project  Thor26  where  C^,  c^,  3^  and  Y  are  empirical  constants. 

Morse  proposed  checking  the  British  hypothesis  of  random  rotation  in 
flight  leading  to  A  in  Equation  (26)  above  by  using  photography.  He 
also  discussed  a  method  of  determining  the  mass  distribution  of  a 
naturally  fragmenting  bomb  or  shell  from  observations  of  the  frequency 
distribution  of  fragment  hole  areas  in  witness  panels. 


Since  proximity  fuzes  had  become  available  by  1943  a  number  of 
reports  appeared  which  discussed  the  optimum  height  of  burst  for  a 
shell  or  bomb  used  against  personnel  on  the  ground.  Birkhoff  and 
Lewy27  compared  the  results  of  three  such  reports  with  experimental 
data  and  found  satisfactory  agreement  on  the  preferred  height.  The 
first  report  they  cited  was  by  Mott26  for  the  British  3.7  inch  shell, 
but  it  is  now  unavailable.  The  second  was  by  Chandrasekhar29  and 
treated  a  105mm  shell.  Chandrasekhar  used  a  coordinate  system  with 


M .  Morse,  R.  Baldwin  and  E.  Kolchin,  " Report  on  the  Uniform  Orienta¬ 
tion  and  Related  Hypotheses  for  Bomb  Fragments,  with  Applications  to 

Retardation  and  Penetration  Problems ",  T.  D.  B.  S.  2,  20  Jan  42. 

26 

Project  Thor ,  TR47 ,  "The  Resistance  of  Various  Metallic  Materials  to 
Perforation  by  Steel  Fragments;  Empirical  Relationships  for  Fragment 
Residual  Velocity  and  Residual  Weight ",  Ballistic  Analysis  Laboratory , 
The  Johns  Hopkins  University,  Apr  61.  (AD# 322781) 

SHght  of  Burst  for  los,m  shM"‘ 

2QN.  F.  Mott ,  Army  Operations  Research  Group  Memo  No.  18,  18  Feb  42.  n 
23S.  Chandrasekhar ,  "Optimum  Height  for  the  Bursting  of  a  105mm  Shell, 
BRL  MR  129 ,  15  Apr  42. (AD#492801) 
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the  x,  y  plane  defining  the  ground  plane  and  z  =  h  representing  the 
height  of  burst.  For  a  horizontal  shell  with  its  axis  parallel  to 
the  x  axis  he  took  the  equatorial  side  spray  of  fragments  to  be 
symmetric  about  the  y,  z  plane  with  an  angular  width  of  2A  (A  =“  7.5°). 
For  prone  (."P°ker  chip")  men  of  area  A  the  area  presented  normal  to 
the  slant  range  R  is  Ah/R  so  the  expected  number  of  hits  is 

=  Nq  Ah/R  .  Chandrasekhar  used  Equation  (51),  which  he  also 

attributed  to  Welch,  for  the  mass  distribution  of  fragments,  but 
he  implicitly  assumed  a  uniform  spatial  distribution  within  the 
narrow  side  spray  so  that  is  an  angular  step  function  which 

depends  on  average  effective  range  as  determined  by  the  mass  distri¬ 
bution,  average  shape,  air  drag  and  lethality  threshold  parameters 
expressed  in  Equation  (52)  above.  His  expression  for  the  number  of 
casualties  was  similar  to  Kent's  Equation  (1) 


C  = 


*  \  Jfdxdy  \tl[’  n 

r  v  *-  min 

1  I  II 


1 


N(m)dm J  Ah/R 


(57) 


where  region  I  was  a  rectangular  area  of  overhitting  instead  of  Kent's 
circular  area  and  T  is  the  number  of  targets  per  unit  area  of  ground. 

n 


One  side  of  the  target  area  was  defined  by  the  extent  of  the  side  spray 
while  the  other  was  defined  by  the  average  effective  range.  The  mass 
integral  in  Equation  (57)  was  divided  by  the  angular  width  2A  553  2  sin  A 
as  well  as  by  a  factor  of  two  because  only  half  the  fragments  were 
projected  at  the  ground.  Chandrasekhar  smoothed  and  extrapolated 
experimental  fragment  mass  data  in  order  to  tabulate/”  N(m)dm  versus  m. 

J  m  _ 

He  then  used  Equation  (45)  as  his  lethality  criterion  and  tabulated  R„ 

b 


in  Equation  (52)  versus  m.  Next  he  tabulated  the  integrand  in  region 
II  of  Equation  (57)  for  various  values  of  h  (50  ft  <  h  <  120  ft)  and 
/~2  2~ 

+y  ,  the  radial  distance  on  the  ground  measured  from  directly  below 
the  burst  point.  Since  A  was  taken  to  be  a  constant  it  did  not  enter 
into  the  tabulation  of  the  integral,  so  that  C/(AT^)  could  be  found  as  a 
function  of  h  by  summing  over  values  tabulated  for  various  distances 

/x2+y2  since  R  =  v4i2+  }t2+y2  .  He  concluded  that  the  optimum  value  for 

2  2 

h  was  about  75  feet  with  little  sensitivity  to  A  =  2.5  ft  or  5  ft  or 
A  =  7.5°  or  15°.  At  this  height  he  estimated  that  each  shell  burst 
would  cover  an  area  about  20  ft  by  150  ft,  indicating  that  such  bursts 
should  be  placed  in  a  grid  pattern  with  these  dimensions  for  optimum 
effect.  Obviously  Equation  (57)  predicts  zero  casualties  for  h  =  0 
since  prone  "poker  chip"  men  are  safe  from  a  burst  at  ground  level  in 
this  model,  so  this  formula  does  not  reduce  to  Kent's  Equation  (1)  which 
envisions  standing  targets.  It  does  however  include  a  mass  distribution 
of  fragments  as  well  as  air  drag  effects  which  vary  with  mass  and  an 
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average  shape,  although,  like  Kent's  formula,  it  is  limited  to  a  partic¬ 
ular  shell  orientation. 

Soon  after,  Morse  and  co-workers30  discussed  the  same  problem  but 
included  foxhole  protection.  In  this  model  a  man  is  considered  barely 
safe  if  the  slant  height  R  =  h/sin  £  where  £  is  measured  in  the  R,  h 
plane  between  R  and  the  ground  plane.  If  £  =  30°  we  have  a  so-called 
”30  degree"  foxhole  with  the  maximum  effective  distance  on  the  ground 
measured  from  a  point  directly  below  the  burst  r  =  /3  h.  Like  Mott, 
Morse  used  the  poisson  expression  for  the  probability  of  at  least  one 
lethal  hit.  For  the  500  pound  general  purpose  bomb  he  constructed 
detailed  contour  plots  of  this  probability  on  the  ground  plane  for 
two  angles  of  fall  (45°  and  60°),  four  burst  heights  (0,  30,  50  and 
100  feet)  and  two  types  of  target.  For  personnel  targets  he  used  the 
energy  threshold  criterion  of  58  foot  pounds,  while  for  truck  targets 
his  threshold  criterion  for  lethal  fragments  was  their  ability  to 
perforate  1/8  inch  of  mild  steel  plate.  For  personnel  he  introduced  a 
terrain  cover  function  (l-r/r^)  multiplied  by  the  presented  area.  He 

computed  the  expected  number  of  casualties  by  placing  men  or  trucks  at 
grid  points  in  a  square  target  area,  for  example,  400  men  ten  feet 
apart,  evaluating  the  probability  of  at  least  one  lethal  hit,  1-  exp 

[-Nj,(x,y)]  ,  at  each  grid  point  and  summing  over  all  grid  points. 

This  eliminated  the  need  to  evaluate  integrals.  He  concluded 

that  a  burst  height  of  about  50  feet  would  be  optimum  against  shielded 

personnel  or  trucks  (shielding  each  other) . 

Morse  soon  followed  this  report  with  another  on  the  effectiveness 
of  the  105mm  shell31.  He  used  the  same  methodology  and  concluded  that 
the  optimum  burst  height  varied  from  30  feet  to  65  feet  as  shielding 
increased  from  natural  terrain  only  to  terrain  plus  45°  foxholes. 

Morse  compared  his  results  with  those  of  Chandrasekhar  mentioned  above 
and  enumerated  the  following  differences  between  the  two  models: 

(1)  Chandrasekhar  assigned  a  vulnerable  area  to  a  prone  man  (in  the 
ground  plane)  while  Morse  represented  a  man  by  a  sphere  which  always 
presents  the  same  area  to  the  radius  vector  R,  Morse  also  added  a 
representation  of  terrain  shielding;  (2)  Chandrasekhar  smoothed  and 
extended  experimental  pit  results  for  the  mass  distribution,  while 
Morse  used  the  results  without  modification;  (3)  Morse  used  the 
experimental  angular  distribution  of  fragments32  while  Chandrasekhar  used 
a  uniform  distribution  over  a  narrow  equatorial  spray;  (4)  Morse  used 


30AL  Morse ,  W.  Transue  and  R.  Kuebler,  " The  Optimum  Point  of  Burst  of  a 
500  lb.  G.  P.  Bomb  Equipped  with  a  Proximity  Fuze",  T.  D.  B.  S,  7, 

22  Apr  43. 

31M.  Morse,  W.  Transue  and  R.  Kuebler,  "On  the  Optimum  Height  Of  Burst 
of  a  105mm  H.  E.  Shell,  Ml,  Used  against  Personnel ",  T .  D.  B.  S. 

11,  1  Jun  43. 

32tf.  A.  Tolch,  "Number,  Angular  Distribution,  and  Velocity  of  Fragments 
from  105mm  H.  E.  Shell  Ml",  BRL  R388,  28  Jul  43.  (AD#491869) 
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an  energy  threshold  of  58  foot  pounds  for  his  lethality  criterion  while 

Chandrasekhar  used  v  =  210  (pound)  ft/sec;  and  (5)  Morse  used 
a  summation  technique  over  a  grid  while  Chandrasekhar  integrated  with 
an  assumed  constant  target  density,  T  In  addition  it  might  be  men¬ 
tioned  that  Morse  included  non-zero  angles  of  fall. 

As  mentioned  above,  Birkhoff  and  Lewy  pointed  out  that  the  con¬ 
clusions  of  Morse  and  Chandrasekhar  concerning  optimum  burst  height 
were  in  reasonable  agreement  with  each  other  and  with  the  results  of 
experiments  performed  for  the  Field  Artillery  Board  at  Fort  Bragg,  all 
employing  the  105mm  shell.  In  addition,  Mott's  conclusions  concerning 
the  British  3.7  inch  shell  were  also  about  the  same. 

"All  studies  agree  that  a  variation  of  ±  20%  from  optimum  in  the 
height  of  burst  will  decrease  the  effectiveness  by  at  most  10%.  Within 
these  limits,  the  results  of  the  four  studies  are  in  agreement  with 
each  other,  in  spite  of  numerous  differences  in  the  mathematical 
formulations  used"2? 

The  underlining  of  the  words  quoted  above  by  Birkhoff  and  Lewy  seems 
to  indicate  their  pleasant  surprise  that  a  variety  of  simple  models 
could  lead  to  the  prediction  of  the  same  useful  result,  confirmed  by 
experiment.  The  experimental  result  was  about  45  feet  for  optimum 
burst  height  (Appendix  to  Chandrasekhar's  report)  while,  model  estimates 
ranged  above  and  below  this  by  as  much  as  78%.  However,  the  particular 
point  of  agreement  noted  above  was  of  a  relative  rather  than  of  an 
absolute  nature,  namely,  that  effectiveness  for  h  values  on  either  side 
of  the  optimum  height  of  burst  did  not  change  very  much  from  the  maxi¬ 
mum  value  of  effectiveness.  Encouraged  by  this  result,  Lewy  later 
published  a  report  on  the  optimum  burst  height  for  demolition  bombs 
which  employ  a  blast  rather  than  fragment  lethal  mechanism33.  Lewy  also 
made  suggestions  on  how  to  handle  terrain  shielding34.  In  addition,  he 
published  a  report  on  the  effectiveness  of  a  300  pound  bomb  against 
personnel,  following  the  methods  of  Kent  and  Chandrasekhar35.  Soon 
after,  Huntoon36  remarked  that  a  dispersion  in  height  of  burst  which 
is  inevitable  for  a  proximity  fuze,  could  actually  be  an  advantage  in 
view  of  the  variety  of  target  conditions  actually  encountered. 

3'3ff.  Lewy }  "On  the  Optimum  Height  of  Explosion  of  Demolition  Bombs' '> 

BRL  MR185j  29  Jun  43.  (AD#493725) 

34tf.  Lewy,  "On  the  Shielding  Effect  of  Small  Elevations  of  Terrain 

from  Explosions  of  Small  but  Variable  Height ",  BRL  MR195 ,  14  Jul  43. (AD#493415) 

35tf.  Lewy ,  "Determination  of  the  Number  of  Casualties  Caused  by  the 
Explosion  of  a  300  lb.  Bomb ",  ,BRL  R403 ,  11  Sep  43. (AD#491874) 

36if.  D.  Huntoony  "Effects  on  Bomb  Damage  of  Dispersion  in  Height  of 
Burst" t  O.  S.  R.  D.  1867 ,  27  Sep  43. 
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C.  Further  Development. 

During  the  next  year  or  so  a  number  of  reports  appeared  which 
provided  additional  background  for  the  further  development  of  our 
subject.  Here  we  will  mention  some  examples.  Thomas^7  pointed  out 

1/2 

that  Equation  (51)  above,  or  its  analog  involving  (m)  ,  can  be 

obtained  by  simple  general  arguments  and  need  not  be  based  on  fracture 
models  involving  cubes,  planes  or  lines.  This  implies  that  the  ability 
to  fit  data  with  a  formula  of  this  type  is  not  an  argument  in  favor  of 
a  particular  breakup  mechanism.  For  example,  there  is  no  need  to 
assume  simultaneous  fracturing  in  any  regular  manner,  since  any  kind 
of  random  sequential  breakup  leading  to  a  decrease  in  average  size 
and  an  increase  in  total  fragment  number  will  do  quite  as  well  as 
any  other . 

Gurney30  used  an  energy  balance  to  derive  expressions  for  the  launch 
speed  of  fragments  from  cylindrical  or  spherical  shells  filled  with 
high  explosive.  When  the  metal  casing  breaks  it  has  expanded  about 
50%  or  so  to  a  radius  r  =  a  and  allowed  the  detonation  products  to 
drop  in  density  to  a  value  p.  If  we  assume  that  the  radial  speed  of 
expansion  of  these  products  is  a  linear  function  of  r,  namely, 
v  =  (r/a)vQ,  with  the  products  in  contact  with  the  case  having  the 

launch  speed  of  the  case  fragments,  vq,  then,  if  we  neglect  all  other 

energy  sinks  such  as  metal  or  gas  heating,  we  can  partition  the  detona¬ 
tion  energy  into  two  parts,  shared  by  the  case  and  the  detonation 
product  gases  as  follows 

EcC  ■  £  Mvo2  *  ff  fl  °c  (rVo/a)2  ]  (4*  r2  dr)  (58) 

where  E  is  the  energy  per  unit  mass  released  by  a  spherical  explosive 
C  4  3 

charge  of  mass  C  =  p  (—  Tra  ) ,  detonated  at  its  center.  Here  M  is  the 

C  j 

mass  of  the  spherical  shell  encasing  the  charge.  If  we  carry  out  the 
integration  in  Equation  (58)  and  solve  for  vq  we  obtain 


VQ  =  1  (2E)  /  ((M/C)  +  .6]  (59) 

If  we  carry  out  the  same  sort  of  analysis  for  an  infinitely  long 

2 

cylindrical  charge  of  mass  per  unit  length  C  =  p  (rca  )  detonated  along 
its  axis  with  M  equal  to  the  mass  per  unit  lengtft  of  the  cylindrical 

? 

metal  shell,  using  volume  element  (2iTr  dr)  instead  of  (4tt  r  dr)  in 
Equation  (58),  we  obtain  Equation  (59)  with  .5  instead  of  .6  in  the 

37 l.  H.  Thomas,  "Comments  on  Mott's  Theory  of  the  Fragmentation  of 
Shells  and  Bombs ",  BRL  R398,  4  Sep  43.  (AD# 36152) 

30fl.  W.  Gurney ,  "The  Initial  Velocities  of  Fragments  from  Bombs , 

Shells ,  Grenades ",  BRL  R405,  14  Sep  43. (AD# 36218) 
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brackets.  In  this  case  vq  approaches  /4Ec  for  a  bare  cylindrical  charge 

(M/C  -*■  o)  or  slightly  less  than  this  for  a  bare  spherical  charge.  If 
this  approximation  were  correct  in  this  limit,  we  would  expect  to  be 
able  to  use  measured  values  for  the  heat  of  detonation  for  Ec  (cal/gm) . 

Gurney  noted  that  for  TNT  filled  shells  with  .18  <  M/C  <  16.67  an 
appropriate  value  for  Ec  was  715  cal/gm,  leading  to  /2E =  8,000  ft/sec, 

which  is  considerably  less  than  the  measured  heat  of  detonation. 

Similarly  he  noted  that  we  cannot  expect  such  an  approximate  formula  to 
agree  with  observation  in  the  limit  as  M/C  becomes  arbitrarily  large. 
Experimentally  v  ->-0  well  before  M/C  00  since  a  heavy  walled  vessel 
will  completely  Sontain  the  detonation  of  a  small  amount  of  explosive 
inside  it.  Still,  over  the  range  of  military  interest  Equation  (59) 
and  its  analog  for  a  cylinder  can  serve  as  a  useful  approximation. 

Gurney  also  noted  G.  I.  Taylor's  discussion  of  the  expansion  velocity 
of  a  thin  walled  cylindrical  bomb  detonated  at  one  end  rather  than 
axially,  but  did  not  compare  the  two  models39. 

Thomas**0  applied  the  cylindrical  analog  of  Equation  (59)  to  a  90mm 
shell  divided  for  computational  purposes  into  ten  cross  sectional  parts 
perpendicular  to  its  axis.  In  this  way  the  actual  variations  of  the 
M/C  ratio  along  the  length  of  the  shell  could  be  approximated  by  using 
Gurney's  formula  with  a  different  M/C  for  each  section.  Thomas  also 
attempted  to  account  for  the  lack  of  parallelism  between  the  axis  and 
individual  sections  of  the  wall  by  simply  doubling  the  actual  angles, 
since  he  assumed  that  the  shell  would  expand  to  about  twice  its  diameter 
before  it  fragmented.  However,  in  order  to  obtain  agreement  with 
panel  data,  he  was  forced  to  assume  a  distribution  about  each  angle 
(normal  distribution  with  standard  deviation  ±  7  degrees) .  In  addition, 
he  fitted  Equation  (51)  to  panel  and  pit  data,  using  the  half  power 
form.  He  was  obviously  not  satisfied  with  this  angular  analysis, 
since  about  six  months  later  he  published  a  report  showing  the 
relation  between  the  Taylor  and  Gurney  models  of  a  cylindrical 
H,  E.  filled  shell1*1.  By  using  an  asymptotic  solution  method  to  take 
account  of  the  transverse  velocity  of  the  detonation  products  in  Taylor's 
theory  he  was  able  to  derive  Gurney's  formula  in  the  final  expanded 
state.  Taylor's  prediction  of  the  angle  of  departure  of  the  fragments 
in  this  limit  remained  the  same,  so  that  this  angle  could  be  added  to 

S^G.  I.  Taylor,  R.  C.  193.  Also  see  " Analysis  of  the  Explosion  of  a 
Long  Cylindrical  Bomb  Detonated  at  One  End”,  Paper  No.  30  (1941)  in 
The  Scientific  Papers  of  Sir  Geoffrey  Ingram  Taylor,  Vol.  Ill,  Aero¬ 
dynamics  and  the  Mechanics  of  Projectiles  and  Explosions,  Ed.  by 
G.  K.  Batchelor,  Cambridge  University  Press ,  1963. 

**UL.  H.  Thomas,  "Analysis  of  the  Distribution  in  Mass ,  in  Speed  and  in 
Direction  of  Motion  of  the  Fragments  of  the  M?1  (90mm)  A.  A.  Shell, 
when  Filled  with  TNT,  and  when  Filled  with  Ednatol ”,  BRL  R434,  3  Dec  43. 

1+1L.  H.  Thomas,  " Theory  of  the  Explosion  of  Cased  Charges  of  Simple 
Shape”,  BRL  R47S,  18  Jul  44.  (AD#491945) 
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the  original  angle  of  the  shell's  double  ogive  shape  to  obtain  a 
better  estimate  of  the  net  angle  of  departure  for  each  section.  This 
type  of  approximation  is  still  in  use  today. 


A  more  exact  theory  which  retains  the  radial  component  of  casing 
velocity  and  leads  to  results  similar  to  Taylor's  has  been  published  by 
Allison  and  Schriempf42. 

Late  in  1943,  Kent  and  co-workers43  used  the  methods  so  far  developed 
to  compare  the  standard  90mm  A.  A.  shell  effectiveness  with  that  of 
several  proposed  modifications  involving  time  plus  contact  fuzing  as 
well  as  filling  with  shapnel  balls  or  eight  ounce  cubes.  The  advantage 
of  using  balls  or  cubes  dispersed  by  spin  and  air  blast  rather  than  by 
a  large  explosive  charge  would  be  their  filling  of  the  region  forward 
of  the  shell  where  the  target  was  likely  to  be  in  contrast  to  the 
expanding  annular  pattern  characteristic  of  standard  shell  design. 

The  disadvantage  was  the  smaller  number  and  slower  speed  of  the  frag¬ 
ments.  They  concluded  that  none  of  the  proposed  designs  would  rival 
the  potential  of  a  proximity  fuzed  shell.  This  report  is  also  note¬ 
worthy  for  its  nine  appendices  in  which  the  model  assumptions  then 
current  were  reviewed.  In  their  Appendix  VII  they  summarized  the 
available  vulnerable  area  estimates  for  aircraft  and  estimated  that 
their  ratios  for  various  fragments  and  aircraft  were  correct  to  about 
±  50%,  while  they  might  be  absolutely  correct  to  within  a  factor  of 
two.  In  their  final  appendix  they  repeated  Mott's  derivation  of  an 
optimum  size  fragment  mass  given  in  Equations  (47)  and  (48)  above, 
apparently  without  being  aware  of  Mott's  work  over  a  year  before  their 
own.  Shortly  after,  Thomas  and  Birkhoff44  employed  similar  methods  to 
conclude  that  a  10%  to  30%  advantage  could  be  gained  by  filling  the  90mm 
A.  A.  shell  with  ednatol  explosive  instead  of  TNT. 

The  personnel  area  target  problem  continued  to  receive  attention  in 
the  years  from  1943  to  1945.  Morse  and  co-workers45  compared  the 
effectiveness  of  the  500  pound  bomb  loaded  with  explosive  in  ten 
different  ways  and  subjected  to  pit,  panel  and  velocity  tests  at 
Aberdeen  Proving  Ground.  In  order  to  do  this  they  employed  two  indices, 
a  space  index,  I,  describing  effectiveness  against  a  collection  of 
targets  surrounding  a  burst  in  mid  air  (A.  A.  Shell  versus  bomber 


42F.  E.  Allison  and  J.  T.  Sehr>iempf3  ", Explosively  Loaded  Metallic 
Cylinders ,  IT",  J.  Appl.  Phys.  31_  (5)  3  846  (I960). 

43<3.  Burkhoffj  R.  R.  Kent  and  D.  R.  Inglis3  " Comparative  Effectiveness 
of  Different  90mm  Shell  Against  Aircraft ",  BRL  R3103  27  Oct  43. 

44L.  H.  Thomas  and  G.  Burkhoff ,  " Use  of  Ednatol  or  TNT  in  A.  A.  Shell ", 
BRL  R4413  23  Dec  43. 


4%.  Morse3  W.  Transue  and  R.  Kuebler3  " Fragment  Damage  Indices.  A 
Comparison  of  Fragment  Effects  of  Different  Loadings  of  the  500  lb. 
G.  P.  Bomb3  AN-M643  and  the  20  lb.  Fragmentation  Bomb3  AM-M41"3 
T.  D.  B.  S.  163  27  Sep  43.  , . 
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squadron),  and  a  plane  index. J,  describing  effectiveness  against  targets 
distributed  in  a  plane  (bomb  or  artillery  shell  versus  troops  on  the 
ground).  Both  indices  were  computed  by  summing  the  trajectory  lengths 
over  which  each  fragment  is  considered  to  be  effective  according  to 
various  threshold  criteria.  For  the  index  J  each  element  in  the  sum 

2 

was  weighted  as  r  rather  than  r  to  account  for  the  slower  increase 
in  number  of  targets  in  a  plane  than  in  a  sphere.  Both  indices  were 
evaluated  for  effective  trajectory  lengths  greater  than  some  radius 
of  overhitting  and  a  more  complete  index,  K,  was  obtained  by  adding 
the  number  of  targets  inside  the  circle  or  sphere  of  overhitting46* 11 7 . 
Morse's  ideas  of  using  fragmentation  data,  postulating  threshold 
criteria  for  effectiveness,  accounting  for  air  drag  losses,  and  so 
forth  were  little  different  from  those  of  his  contemporaries.  He 
differed  from  them  mostly  in  the  quantities  he  chose  to  sum  together. 
Later  in  1945,  Morse's  method  of  making  his  conclusions  about  the 
superiority  of  one  or  another  type  of  shell  was  criticized  by  Scheffe46. 
Scheff4  admitted  that  he  was  in  the  uncomfortable  position  of  raising 
difficulties  without  having  anything  better  to  offer.  Still  his 
criticisms  are  worth  considering  since  some  of  them  apply  to  all 
writers  on  this  subject,  including  those  working  in  the  field  today. 

He  pointed  out  that  assuming  that  the  possession  of  a  certain  energy 
or  power  to  perforate  a  given  thickness  of  steel  or  wood  implies  the 
ability  to  inflict  a  certain  type  of  damage  on  a  man  or  materiel 
target  requires  a  leap  in  our  chain  of  reasoning.  How  are  we  to 
justify,  for  example,  our  assumption  that  ability  to  perforate  1/8  inch 
of  mild  steel  correlates  with  the  complex  phenomenon  of  damaging  a 
plane  or  truck  in  a  certain  way?  To  attempt  to  do  this  by  calculation 
seems  not  only  very  complicated  but  also  highly  uncertain.  To  accom¬ 
plish  it  experimentally  would  seem  to  involve  prohibitive  expenditures 
of  time  and  materiel.  Even  apart  from  these  uncertainties  Scheffe 
pointed  out  that  uncertainties  in  the  experimental  fragmentation  data 
upon  which  all  such  calculations  are  based  carry  over  into  any 
conclusions  which  are  drawn.  He  advanced  some  general  considerations 
which  indicated  that  a  difference  in  the  Morse  indices  should  not  be 
regarded  as  statistically  significant  unless  it  exceeded  30  per  cent. 
Since  the  differences  calculated  by  Morse  were  almost  always  smaller 

46Af.  Morse,  W.  Transue  and  R.  Kuebler ,  "The  Choice  of  Bombs  According 
to  the  Type  of  Fragment  Damage  Required.  New  Methods  for  Comparison 
of  Different  Missiles  as  Well  as  the  Same  Missile  with  Different 
Loadings' T.  D.  B.  S.  18 ,  23  Dec  43. 

47M.  Morse ,  W.  Transue  and  M.  Reins ,  "The  Expected  Number  of  Effective 
Fragment  Hits  from  Ground  Bursts  of  Bombs  Taking  Account  of  the  Angle 
of  Fall  and  Remaining  Velocity  of  the  Bomb,  and  the  Average  Height  of 
Targets",  T,  D.  B.  S.  49,  5  Jan  45. 

48fl.  Scheff&i;  "Some  Comments  on  the  Morse  Theory  of  Fragment  Damage", 
OSRD  6426,  3  Dec  45. 
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than  this,  the  implication  was  that  the  conclusions  drawn  concerning 
the  superiority  of  one  shell  type  over  another  were  not  meaningful. 

For  example,  the  most  widely  fluctuating  measurement  which  entered 
into  the  calculations  was  the  number  of  fragments  observed  in  the  pit 
experiments  which  led  to  a  standard  deviation  of  13%.  If  at  least  two 
standard  deviations  were  required  to  distinguish  the  effectiveness 
of  one  explosive  fill  compared  to  another,  we  see  that  the  calculated 
differences  were  not  usually  statistically  significant. 

I 

Sheffe,  pointed  out  how  experimental  uncertainties, namely,  a  lack 
of  reproducibility  in  a  certain  measurement  technique,  lead  to  uncer¬ 
tainties  in  calculated  effectiveness  parameters.  He  also  noted  that 
the  sand  pit  technique  itself  might  contain  a  built  in  bias  because 
of  secondary  break-up  and  referenced  a  report  by  Gurney  and  Lewy1*9 . 
These  authors  remarked  that  the  fragmentation  of  bombs  had  been  re¬ 
investigated  using  a  set  up  in  which  the  sand-pit  was  covered  with 
four  inches  of  celotex  or  with  fourteen  inches  of  sawdust  since  it  was 
feared  that  in  earlier  experiments  many  of  the  high  velocity  fragments 
had  been  further  broken  up  by  impact  on  sand,.  The  results  indicated 
that  for  the  500  pound  bomb  the  fragments  were  more  than  twice  as 
heavy  and  less  than  half  as  numerous  as  in  the  earlier  sand  pit 
experiments,  while  for  the  smaller  100  pound  bomb  the  fragments  were 
about  twice  as  heavy  and  half  as  numerous  as  before.  Again  the 
question  of  secondary  breakup  in  celotex  or  sawdust  could  be  raised, 
but  these  authors  did  not  do  so.  Instead  they  commented  that  the 
observed  changes  in  fragmentation  did  not  seem  to  affect  calculations 
of  effectiveness  very  much  since  the  reduction  in  the  number  of  frag¬ 
ments  is  compensated  by  an  increase  in  their  range.  For  example, 
for  the  100  lb  bomb,  separate  calculations  showed  that  doubling  the 
average  mass  and  halving  the  number  of  fragments  produced  a  change 
of  fifteen  per  cent  in  the  effectiveness  index.  The  question  of 
breakup  in  softer  recovery  media  was  studied  later  by  Bentz50  who 
concluded  that  pit  data  gives  a  number  of  fragments  at  least  several 
times  the  true  value. 

The  British  also  continued  an  active  interest  in  estimating  the 
the  effectiveness  of  anti-personnel  weapons.  In  1944,  Pearson  and 
co-workers51  introduced  the  name  ’’lethal  area"  and  pioneered  the  use 

W.  Gurney  and  H.  Lewy3  "On  the  Comparison  of  the  Efficiency  of 
the  20  lb.  Fragmentation  Bomb3  the  23  lb.  Parachute  Borib3  the  100  lb. 
G.  P.  Bombj  and  the  500  lb.  G.  P.  Bomb  Singly  and  in  Clusters  for 
Attack  by  Grvund  Bursts-on  Aircraft  Dispersed  on  Open  Ground "3  BEL 
MR289 3  24  Apr  44.  (AD#492802) 

50 W.  H.  Bentz,  "Secondary  Fragment  Break-up  in  Cane  Fiber 3  Cardboard 3 
and  Sandy"  BEL  E7043  3  Aug  49.  (AD#492317) 

51£.  S.  Pearson3  N.  L.  Johnson  and  D.  F.  Mills 3  "The  Calculations  of 
Lethal  Areas  for  Anti-personnel  Weapons3  with  Special  Eeference  to 
the  25-pr.  Shell"3  A.  C.  58853  9  Mar  44. 
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of  a  lethality  criterion  which  was  not  a  simple  step  function  like 
possession  of  58  ft-lbs  of  energy  or  ability  to  perforate  a  one-inch 
thick  board.  The  lethal  area,  A  ,  they  interpreted  as  follows.  If 

Li 

T  targets  are  distributed  on  average  uniformly  over  an  area  S»A  , 

such  that  =  T/S  is  the  target  density,  then  on  average  it  is  to  be 

expected  that  AT  T.  =  A.  T/S  targets  will  be  incapacitated.  The 
LA  L 

lethal  area  was  computed  by  integrating  the  damage  function  D(x,  y) 
over  all  possible  ground  positions.  Thus 

Al  -JJ ' D(x,  y)  dx  dy  =  j  dx  dy  (60) 

using  the  poisson  distribution  as  explained  above.  We  see  that  this 
notion  is  similar  to  Kent's  Equation  (1)  above  if  we  rewrite 
Equation  (1)  in  the  form  C  =  T.  A. .  The  major  conceptual  difference 

-N  A  L 

is  the  use  of  (1-e  )  for  the  damage  function  instead  of  adding  an 

area  integral  to  a  circle  of  100%  effectiveness  as  in  Equation  (1) . 

We  have  seen  that  Kent  himself  eventually  adopted  the  poisson  form 

(see  Equation  (13)  above) .  From  Equation  (60)  and  the  assigned 

interpretation  C  =  A.  T. ,  it  is  clear  that  a  "lethal  area"  may  have 

L  A 

the  dimensions  of  an  area,  but  it  is  not  a  spatial  region  on  the 
ground  which  can  be  paced  off  or  staked  out.  Rather  it  is  an 
effectiveness  parameter  or  casualty  index  which  when  multiplied  by 
the  number  of  targets  per  unit  area  gives  the  number  of  casualties  to 
be  expected  on  average. 

Pearson  also  suggested  using  the  "Zuckerman  man"  instead  of  the 
"1. 0-inch  wood  man",  by  which  he  meant  a  vulnerable  area  or  kill 
probability  which  varied  continuously  from  a  minimum  to  a  maximum 
value.  However,  he  proposed  no  functional  form  for  this  variation, 
remarking  that  it  should  be  determined  experimentally.  In  practice  he 
used  the  "1.0- inch  wood  man". 

Carlton20,  in  the  report  already  mentioned  above,  introduced  a 

different  meaning  for  "lethal  area",  considering  instead  an  attack  on 

a  single  (aircraft)  target  by  a  proximity  fuzed  shell.  As  mentioned 

before,  the  burst  occurs  on  the  cone  0=0  at  a  distance  z  =  R  cos  6 

o  o 

in  front  of  the  aircraft  dictated  by  the  sensitivity  of  the  fuze. 

On  the  assumption  of  symmetry  about  the  z  axis,  Carlton  discussed  the 
conditional  damage  probability  D(r)  for  fuzes  with  radial  miss  r 
found  by  averaging  D(r,  z)  over  the  values  of  z  which  occur  in  the 
burst  pattern.  He  then  defined  the  "lethal  area"  or. 'Effective  size" 
of  the  target  as  the  integral  of  D(r)  in  the  x,  y  plane,  namely, 
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(61) 


ALC  =  2ir/o  D(r)  T  dr  ' 

In  addition,  he  discussed  the  probability  of  damaging  the  target, 
given  functioning  fuzes,  as 

-  r 

P(a)  -  K  /  D(r)  e  2ff  r  dr  (62) 

a  Jo 

which  becomes  Equation  (38)  above  if  we  use  Equation  (35)  for  D(r) . 

Here  Equation  (61)  is  Equation  (62)  when  all  miss  distances  are 

equally  probable  instead  of  following  a  Gaussian  distribution  so  that 

in  Equation  (61)  we  are  assuming  a  uniform  distribution  of  burst  point 

projections  on  the  x,  y  plane  which  contains  a  single  target  at  the 

origin,  whereas  in  Equation  (60)  above  we  are  assuming  a  uniform 

distribution  of  targets  in  the  x,  y  plane  with  a  single  burst  point 

projected  onto  the  origin.  Carlton's  type  of  "lethal  area"  has  been 

labeled  A  „  to  distinguish  it  from  A  in  Equation  (60).  The  difference 
LL  Li 

is  conceptual  rather  then  computational.  In  either  case  we  are  dealing 
with  a  "casualty  index  or  "effectiveness  parameter"  which  has  a 
meaning  only  in  a  probabilistic  sense.  Because  of  this,  the  results 
of  a  single,  real  burst  may  differ  widely  from  the  estimate  computed, 
even  if  the  encounter  is  modeled  in  great  detail.  Such  estimates  as 
we  are  dealing  with  have  significance  only  for  the  strategist  who 
considers  repeated  encounters  of  the  same  type.  This  was  pointed 
out  by  Cunningham  in  1940  in  the  quotation  we  cited  above1  .  Obviously 
we  should  not  use  Equation  (61)  when  all  miss  distances  are  not  equally 
probable  as  usually  occurs  when  a  single  target  is  being  aimed  at. 

Only  if  a  is  much  larger  than  some  finite  effectiveness  limit  will 
Equation  (61)  be  a  valid  approximation,  provided  the  upper  limit  is  ^ 
instead  of  infinity.  The  Carlton  type  of  lethal  area  calculation  for  a 
single  target  is  still  widely  used,  sometimes  without  realizing  that  it 
assumes  either  a  single  burst  near  an  infinite  plane  of  uniformly  dis¬ 
tributed  targets  or  a  single  target  near  an  infinite  plane  of  uniformly 
distributed  bursts.  At  least  it  gives  a  lower  bound  for  an  estimate  of 
relative  effectiveness. 

In  1945  Lewy  and  Gurney52  generalized  Kent's  formula,  Equation  (1), 
retaining  the  idea  of  a  radius  of  overhitting,  but  considering  a  burst 
height,  hi,  above  the  ground.  In  addition,  they  approximated  the  angular 
distribution  of  fragments  from  a  shell  by  an  equatorial  side  spray  of 
angular  width  2A  superimposed  on  a  spherical  spray,  each  spray  with  its 
own  constant  number  of  fragments  per  steradian,  N^.  The  effective¬ 
ness  of  the  spherical  spray  is,  of  course,  independent  of  the 

52fl.  Lewy  and  R.  W.  Gurney ,  "Rapid  Methods  for  the  Prediction  of 
Damage  by  Fragments  in  its  Dependence  on  the  Variables  of  the  Tactical 
Situation ",  BRL  R5Z03  14  Feb  45. (AD#PB165304) 
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attitude  of  the  projectile  at  burst.  However,  the  effectiveness  of  the 
side  spray  depends  not  only  on  its  width  and  the  height  of  burst  but 
also  on  the  angle  of  fall  to  as  defined  in  Equation  (8)  above.  A 
symmetric  side  spray  is  bounded  by  the  two  nappes  of  a  right  circular 
cone,  each  of  half  angle  (tt/2-A),  which  is  the  co-latitude  measured 
from  the  axis  of  the  shell  (or  tangent  to  its  trajectory  at  burst). 

The  solid  angle  occupied  by  the  side  spray  is  4ir  minus  twice  the  area 
of  spherical  cap  intercepted  by  one  nappe  of  the  cone.  From  Equation 
(2)  this  is  ft  =  4*  -  4tt  [1-cos  (tt/2-A)]  =  4ir  sin  A.  If  we  use  a  result 
obtained  by  Menaechmus  and  Appolonius  in  the  fourth  and  third  centuries 
B,  C.,  we  note  that  the  intersection  of  the  cone  with  the  ground  plane 
forms  a  conic  section  which  is  an  ellipse  of  u  >  ir/2  -  A,  a  parabola 
if  u)  *  tt/2  -  A  and  a  hyperbola  if  u  <  ir/2  -  A.  Since  the  hyperbola  is 
a  common  case  it  is  shown  in  Figure  2.  Also  shown  in  Figure  2  is  a 
circle  of  radius  r2  which  is  the  maximum  radius  of  lethality  on  the 

ground.  Lewy  and  Gurney  considered  "SO0  foxhole"  shielding  so  that 
r2  =  /3h.  In  Figure  2  this  circle  is  shown  intersecting  both  branches 
of  the  hyperbola  and  limiting  the  potentially  lethal  area  of  the 
side  spray  to  the  region  without  grid  lines.  Of  course  this  circle 
might  intersect  only  one  branch  of  the  hyperbola  or  not  intersect  it 
at  all. 

Let  us  consider  a  right-handed  coordinate  system  in  Figure  2  with 
origin  at  the  burst  point,  z^-axis  pointing  along  the  projectile  ve¬ 
locity  vector  at  the  moment  of  burst  and  y^-axis  pointing  up  from  the 
Figure.  Then  the  equation  of  the  cone  is 


(x1)2  +  (y1)2  =  (z1)2  tan2  (01)  =  (z1)2  /  tan2 A  (63) 

1  12 

since  the  co-latitude  0  =  it/2- A .  If  we  add  (z  )  to  both  sides  of 

Equation  (63)  the  cone  equation  becomes 


f  I'i^  -  1.2  ,  1.2  f  1.2  .  .  2  . 

(x  )  +  (y  )  +  (Z  )  =  (z  )  7  sin  A  . 


(64) 


We  may  express  this  equation  in  anotherjCoordinate  system  obtained  by 
a  counterclockwise  rotation  about  the  y  -axis  through  an  angle  u+ir/2, 
namely 


x^  =  -  x  sin  to  +  z  cos  to 


y1  =  y 

z*  =  -  x  cos  to  -  z  sin  to 


(65) 


Now  the  equation  of  the  cone  is 


44 


In  this  coordinate  system  the  equation  of  the  ground  plane  is  z  =  -h 
and  the  equation  of  the  conic  section  mentioned  above  is  obtained  by 
letting  z  =  -h  in  Equation  (66).  If  we  do  this,  divide  by 
2  2 

(cos  to  -  sin  A)  and  add 


1 

. 


sin  to  cos  to 

2  2 
(cos  to- sin  A) 


] 


(67) 


to  both  sides,  the  equation  of  the  conic  section  is 


u-v' 


2 

sin  A 


2  2 
(cos  to-sin  A) 


sin  A  cos  A 

2  2 
(cos  to-sin  A) 


=  a 


(68) 


Equation  (66)  is  the  same  as  Equation  (9)  above  since  0  =  tt/2  -  A  so 
cos  6  =  sin  A.  The  conic  section  of  Equation  (68)  is  centered  at  the  pint 
(xQ,0,-h).  If  we  divide  Equation  (68)  by  a2,  we  put  it  into  standard 
form,  namely 


with 


1 _ 

2 

a 


C69) 


2  2  2 
cos  A  /  (cos  to  -  sin  A) 


(70) 


Here  a,  defined  in  Equation  (68) ,  is  the  serai-transverse  axis  and  b 
is  the  semi- conjugate  axis.  The  eccentricity  is 


e=  ^  a^+b^  /  a  =  cos  u)  /  sin  A  (71) 

and  its2relation  to  the  discriminant  6  of  Equation  (10)  above  is 
6  =  4(e  -1).  As  stated  before,  for  oj  <  (tt/2-A)  we  have  a  hyperbola 
(e>l),  for  ui  =(/t/2-A)  we  have  a  parabola  (e=l),  and  for  u)>^tt/2-A)  we 
have  an  ellipse  (e<l)  which  becomes  a  circle  if  to  =  tt/2  (e=o) ,  the 
case  of  a  vertical  fall. 


For  a  vertical  fall  (ui=tt/2),  xq  =  0  in  Equation  (65),  and 
Equation  (68)  is 


x 2  +  y2  -  (h/tan  A) 2 


(72) 


so  the  bounding  conic  section  is  a  circle  of  radius  h/tan  A  centered 
directly  below  the  burst  point.  No  side  spray  fragments  fall  inside 
this  circle  and  no  casualties  are  possible  due  to  this  spray  unless 
h/tan  A  <  r  .  If  h  -*■  o  this  circle  of  exclusion  shrinks  to  a  point 
and  we  return  to  the  case  considered  by  Kent.  The  circle  of  exclusion 
also  shrinks  to  a  point  as  A  ir/2  since  this  represents  a  spherical 
fragment  distribution.  In  this  case  it  is  meaningless  to  distinguish 
between  an  equatorial  side  spray  and  a  uniform  spherical  background. 

For  a  horizontal  fall  (an  unlikely  event) ,  w  =  o  and  the 
eccentricity  e  =  1/sin  A  >  1  for  A  <  7t/2,  so  we  have  a  hyperbola 
with  xQ  =  0  from  Equation  (b7) .  The  side  spray  is  concentrated 

directly  below  the  burst  point.  We  can  re-write  Equation  (68)  as 


/ 1  \./  h2 ) 

i  r  i  i  h2 

(73) 

(e2-1)  (e2-1  ) 

1  2  2  2 
[tan  A (esc  A- 1)  e  -1 

when  o)  =  o.  Now  if  h  -*■  0,  the  boundary  degenerates  into  two  straight 
lines  given  by 


y  = 


(74) 


If  A  -*■  tt/2  (uniform  spherical  spray)  for  w  =  o,  the  branches  of  the 
hyperbola  move  far  apart  since  a  «  in  Equation  (68) ,  so  the  circle 
of  lethality  of  radius  r?  fixed  by  shielding  becomes  the  only  limiting 
factor. 


Lewy  and  Gurney  introduced  the  approximation 


Nfl(R)  =  exp  (-b  -  cR)  (75) 

for  the  number  of  effective  fragments  per  steradian  where  the  constants 
a  and  b  were  evaluated  by  fits  to  experimental  data  over  the  range 
h  <  R  <  2  h  as  discussed  by  Thomas53.  For  prone  men  they  used  the 

2 

"poker  chip"  approximation  (Ap/R  )  (h/R)  for  the  solid  angle  subtended 

*  2 

by  each  target  For  unshielded  standing  men  they  used  (Ap/R  )(r/R)„ 

In  the  case  of  a  ground  burst  r  =  v4P-h2  -►  R  and  Kent’s  integral  in 
Equation  (1)  above  was  found  to  involve  a  tabulated  function,  if  we 
use  Equation  (75)  for  N^,  namely,  with  t  =  cr. 


53L.  H.  Thomas j  " Computing  the  Effect  of  Distance  on  Damage  by 
Fragments' '3  BRL  R468,  18  May  44.  ( AD#492590) 
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2tt  Ap  e 


A,  <•- 


/t)  dt  =  2?r  Ap  e  bEx  (cr^ 


(76J 


where  Ej  is  an  exponential  integral.  Similarly  for  air  bursts  against 
unshielded  prone  men,  we  have 


2tt  h 


Ap 


c  e 


’A"  ( 

/ cR  .  t 
J  min 


(e'Vt2)  dt 


[ 


■] 


2i  h  AD  e  °/R  .  E_  (cR  .  )  (77) 

P  min  2  v  min' 


using  t  =  cR.  Again,  E  is  an  exponential  integral.  For  shielded 
prone  men  Equation  (77j  has  a  finite  upper  limit  and  numerical  integra¬ 
tion  must  be  used.  Clearly  Equation  (77)  or  the  more  general  case 
including  shielding  vanishes  as  h  ■+•  0  implying  no  casualties  if  all 
targets  are  prone  or  in  foxholes  when  a  ground  burst  occurs.  For  given 
id  and  A,  the  number  of  casualties  increases  as  h  increases,  passes 
through  a  maximum  and  decreases  as  h  continues  to  increase  since 
more  fragments  become  ineffective  due  to  air  drag  and  more  fragments 
fall  outside  the  circle  of  radius  P.  =  /3h  determined  by  shielding. 

Lewy  and  Gurney  illustrated  their  modified  version  of  Equation  (1) 
for  a  500  lb  bomb  with  remaining  velocity  of  1000  ft/sec  bursting  with 
angle  of  fall  id  =  55°  over  men  in  30°  foxholes  and  found  an  optimum 
burst  height  near  h  =  60  feet.  They  plotted  the  number  of  casualties 
divided  by  the  target  density,  namely,  C/T^,  versus  h,  a  quantity 

which  they  did  not  name,  but  which  we  have  seen  had  already  been 
dubbed  a  "lethal  area"  by  Pearson  (see  Equation  (60)  above) .  It  is 
a  casualty  index  which  applies  to  a  single  burst  over  a  plane  contain¬ 
ing  a  uniform  distribution  of  targets. 

We  will  close  this  section  of  our  review  with  a  brief  account  of 
a  report  by  Whitcomb54  who  modified  Kent's  original  formulation  expressed 
in  Equation  (4)  above  in  order  to  broaden  the  discussion  of  optimiza¬ 
tion.  Since  Equation  (3)  set  equal  to  unity  defines  the  "radius  of 
overhitting1,'  R  ,  we  have 


R^  =  n  ^  -  N  a  -  M/.m.  a 

1  "ft  AP  ft  “p  ft  AF 


(78) 


where  the  total  number  of  fragments  N  is  equal  to  the  mass  of  the 
shell  case,  M,  divided  by  the  mass  of  a  fragment,  m,  when  all  fragments 
are  alike.  From  Equation  (43)  we  have 


R  =  a  m 


1/3 


m 


(79) 


0HS.  E.  Whitcomb 3  " Fragmentation  Efficiency  of  Bombs ",  BRL  R53Z, 
5  Mar  45. (AD#492432) 
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“  (m/nnf) 


1/3 

with  the  minimum  mass,  m'  ,  able  to  perforate  and 

m 


where  y  = 

damage  a  target' is  given  by  a  threshold  criterion  like  Equation  (56) 


-  I 


.-3/4  =  10ci  Aal,  so  the  target 


with  d,  =  3/4,  3^  =  -  — ,  =  0  and  y 

steel  equivalentthickness  is 

T  1/3  4/3  .  1/3  4/3 

T  =  y  m  v,  =  y  m  v  (80) 

L  m  m  k  ' 

defining  and  the  corresponding  maximum  striking  speed,  v^.  The 

parameter  x  is  defined  from  Gurney’s  cylinder  formula,  the  analog  of 
Equation  (59)  above,  namely,  with  .5  instead  of  .6. 


_/\>V .  V  . 

~\vm/  ~  2E  C+2M 


(81) 


so 


kJLm 

i+x 


~~  =  where  W  is  the  warhead  mass  the  sum  of  the 


explosive  charge  mass,  C,  and  the  case  mass,  M.  The  number  of  frag¬ 
ments  can  now  be  written 


N  =  «  = 

m 


M  /  W  \  jn  =  _W 

m  1 M+C  In?  irt 

\  /  m  m 


W 


(82) 


y5(l+x)/(l~)C)  V 


w 

so  —  =  M  w/m.  If  we  divide  the  square  of  Equation  (79)  by  (78) 

and  use  Equation  (82)  for  we  find 


wz  (83) 
P 


where  z 


=  [y  In  (xly 


1/2  .  3/4\  I2 


y  '  ^  and  P  =  W  Ap  /  [fimm5/3  a2]  . 


Whitcomb 


used  as  his  measure  of  efficiency  Equation  (4)  above  divided  by  Ap 
(W /Jn  )  which  gives  an  efficiency 

2 

c/  [TA  AP  wv]=  ywiT)  [‘  * ln  ( 57 )Z] 

ir(H/ra)  Ap/fl  p  1  „  p  -j 

=  Ap  M  w/m  •  L1  +  ln  M  =  *  L1  + 


(84) 
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In  this  treatment  he  is  considering  ground  bursts  only  so  r  =  R  (since 
h  =  o)  as  well  as  vertical  angle  of  fall  only.  However,  by  varying 
the  parameters  w,  z  and  P  in  Equation  (84)  he  tried  to  discuss  the 
effect  of  size  and  design  changes.  Since  P  is  proportional  to  weight 
he  could  use  it  to  compare  the  efficiency  of  two  bombs  of  the  same 
design  but  different  size.  However,  both  w  and  z  depend  on  both  the 
C/M  ratio  (through  x)  and  the  fragment  size  (through  y) ,  so  his 
formulation  does  not  lend  itself  to  a  clear  discussion  of  the  effect 
of  changing  these  two  design  parameters  even  for  vertical  fall  and 
ground  burst.  By  plotting  a  number  of  curves  with  only  one  parameter 
varied  at  a  time,  he  did  however  conclude  that  for  an  optimum  bomb 
weight  relative  to  a  given  type  of  target,  there  should  exist  both 
an  optimum  C/M  ratio  and  an  optimum  fragment  size. 

III.  EFFECTIVENESS  STUDIES  AFTER  WORLD  WAR  TWO 

A.  The  Refinement  of  Analytical  Methods. 

Soon  after  World  War  Two  a  large  number  of  aircraft  became 
available  for  experimentation  and  hundreds  were  expended  for  this 
purpose.  Assessments  were  made  of  various  levels  of  damage  done 
during  controlled  test  firings  and  these  judgments  were  expressed  as 
component  vulnerable  areas55760.  Much  of  this  work  was  directed  toward 
recommending  families  of  antiaircraft  and  artillery  weapons61-62, 

K.  Weiss  and  A.  Stein ,  " Airplane  Vulnerability  and  Overall 
Armament  Effectiveness , "  BRL  MR462 ,  May  47. (AD#9093) 

5SJ.  M .  Garmousakisj  "The  Effect  of  Blast  on  Aircraft ,  " 

57H.  X.  Weiss  and  A.  Stein ,  "Vulnerability  of  Aircraft  to  75mm  Air- 
burst  Shelly"  BRL  MR  482t  Sep  48.  (AD#41431) 

584,  Stein  and  H.  Kostiak ,  "Damage  by  Controlled  Fragments  to  Aircraft 
and  Aircraft  ComponentSy"  BRL  MR4873  Feb  49.  (AD#53127) 

59X.  K.  Weiss ,  J.  Christian  and  L.  M.  Peters ,  ''Vulnerability  of  Aircraft 
to  105mm  and  75mm  HE  Shelly"  BRL  R687y  Mar  49.  (AD#54768) 

604.  Stein  and  H.  Kostiaky  " Methods  for  Obtaining  the  Terminal  Ballistic 
Vulnerability  of  Aircraft  to  Impacting  Projectiles  with  Application  to 
Generic  Jet  Fighter3  Generic  Jet  Bombery  F-47  Piston  Fighter  and  B-50 
Piston  Bombery"  BRL  R768y  Jun  51.  (AD#130730) 

6V.  H.  Frazer y  H.  P.  Hitchcocky  R.  H.  Kenty  F.  G.  King3  J.  R.  Laney 
T.  E.  Stem  and  H.  K.  Weiss ,  "A  Study  of  a  Family  of  Antiaircraft 
Weapons BRL  TNU93  Sep  49.  (AD# 378732) 

62f.  E.  Grubbs y  R,  H.  Kenty  J ,  R,  Lane  and  H.  K.  Weiss ,  "A  Family  of 
Field  Artillery BRL  R771y  Jul  5%.  (AD#377190) 
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The  vulnerability  of  fuel  and  ammunition  was  given  particular  atten¬ 
tion6  3  »61*.  Since  similar  work  was  being  done  elsewhere,  a  conference 
on  aircraft  vulnerability  was  held  in  1949. 65 

In  1946,  Fano66  had  concluded  that  reliable  effectiveness  estimates 
could  not  be  made  "because  of  insufficient  knowledge  of  the  vulnerability 
characteristics  of  targets".  By  1949  methodology  studies  began  to 
appear  once  more.  Weiss67  discussed  the  Lewy  and  Gurney  approximation 
in  Equation  (73)  above.  He  used  the  1/2  power  form  of  Equation  (51)  to 
obtain  for  the  number  of  fragments  with  mass  greater  than  m 


N 


o-C'V 1/2  d  r(m/n,0)1/2' 


N  e-(m/,,o)1/2 


He  also  assumed  a  threshold  lethality  criterion 


(85) 


m  v*  =  m  v  *  =  L.  (86) 

mo  1  v  ' 

where  y  and  L.  are  constants  and  m  is  the  minimum  mass.  If  *  3  1  we 
l  m 

have  a  momentum  threshold  like  Mott's  in  Equation  (44)  above,  or  if 
Y:  “  2  we  have  Morse's  energy  threshold  as  in  Equation  (53).  In 
addition  he  made  the  usual  assumption  as  in  Equation  (27)  above,  that 


v 


(87) 


If  we  let  m  =  m  in  Equation  (85) ,  we  have  N  =  N 
m  mo 

dividing  Equation  (85)  by  we  obtain  the  ratio 


(rti  /m  ) 
e  v  m  o 


1/2 


and 


63i4.  Stein  and  M.  G.  Torsch ,  "Effectiveness  of  Incendiary  Ammunition 
Against  Aircraft  Fuel  Tanks,"  BEL  ME484,  Oct  48,  (AD#42383) 

64A:.  S.  Jones ,  " A  Comparative  Study  of  the  Vulnerability  of  Warheads 
Loaded  with  TNT  and  with  Composition  B,"  BEL  ME486,  Jan  49. (AD#802306) 

6 5 "Report  on  First  Working  Conference  on  Aircraft  Vulnerability ," 

BEL  ME488,  Mar  49. (AD# 377 164) 

66V.  Fano j  "Discussion  of  the  Optimum  Characteristics  of  Weapons  for 
Most  Efficient  Fragmentation, "  BEL  E594,  Jan  46.  (AD#492600) 

67tf.  K.  Weiss,  "Justification  of  an  Exponential  Fall-off  Law  for 
Number  of  Effective  Fragments , "  BEL  R697,  Feb  49.  (AD#801777) 


N/N  =  e 
m 


/m_\1/2  /.N1''2 

v-)  Mi 


~/m  \1/2r/--  \1/2 
elnr)  1st) 

\  o /  \  m/ 


-1 


(88) 


If  we  put  Equation  (87)  into  Equation  (86) ,  take  the  square  root  then 

1/3 

the  logarithm,  then  multiply  byf — J  ,  we  find 
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where  the  second  equality  is  approximately  true  for  1  < 
within  10%.  From  Equation  (89)  we  find  m— 
which  we  may  substitute  in  Equation  (88)  to  obtain 


i  15  to 
-1/3 


m 


(90) 


where  C  =  .45  ay  m  m0  If  we  divide  Equation  (90)  by  to 

obtain  =  N/Q,  we  see  that  Lewy  and  Gurney's  approximation  in 

Equation  (73)  above,  based  on  particular  experimental  data,  should 
hold  approximately  for  any  shell  to  which  Equations  (85),  (86)  and 
(87)  can  be  applied,  at  least  in  the  range  of  (m/m  )  specified.  Thus 

m  b 

Weiss  was  able  to  generalize  Equation  (73)  ,  relating  eu  and  c  to  the 
parameters  of  the  problem.  The  same  could  be  done  using  the  1/3  power 
law  of  Equation  (51) or  any  other  power.  It  is  interesting  to  note  the 
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expected  number  of  lethal  hits  on  vulnerable  area  Av(0)  according 
to  Equation  (90)  or  (73),  namely 


1  *  nh  V«vr2 


N  A  (0)  e"cR 
m  v*-  _ 

2 

ft  R^ 


^  e_CR/R2 


(91) 


which  has  the  same  form  as  Equation  (31)  above  if  we  neglect  the 
angular  dependence  of  K(i|>)  and  take  to  be  constant.  As  we  recall, 

the  authors  of  OSRD  738  arrived  at  this  form  by  using  Equation  (30)  to 
approximate  the  behavior  of  Av  in  Figure  1  while  Weiss,  Lewy  and 
Gurney  obtained  it  as  an  approximation  to  N  or  N  .  In  either  case  the 

-CR 

dependence  on  R  is  clear.  The  factor  e  involves  loss  of  effectiveness 

because  of  air  drag  (ot=  C  p  k)  as  well  as  a  damage  threshold  criterion 

(™  )  and  a  fragment  mass  distribution  (N  and  m  )  .  The  factor  — ■„  enters 
m'  v  o  o'  r2 

because  of  the  assumption  that  the  fragments  emanate  from  a  point  in 
space,  the  burst  point.  The  solid  angle  ft  may  be  4ir  for  full  spherical 
symmetry,  or  2tt  for  a  hemisphere  as  in  Equation  (31)  or  some  lesser 
amount  for  a  more  restricted  fragment  spray. 

Later  in  the  same  year  Weiss  published  another  report60  detailing 
modifications  to  the  methods  already  developed  in  World  War  Two.  For 
a  uniform  side  spray  of  half  angle  A  thrown  forward  to  co-latitude 
6jj  by  the  remaining  velocity  of  the  shell  Weiss  used  ft  =  4ir  sin  A 
sin  0  .  For  very  high  altitude  bursts  against  aircraft,  Weiss  pointed 

out  that  C  0  in  Equation  (91)  so  that  Equation  (32)  becomes 
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B8//.  K.  Weiss,  "A  Method  for  Computing  the  Probability  of  Killing  a 
Multiply  Vulnerable  Aircraft  Target  with  'N'  Rounds  of  Fragmenting 
Shell,"  BRL  MR495,  12  Sep  49.  (AD#24574) 
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for  fixed  fuze  angle  0  ,  neglecting  the  integral  from  Rq  to  00 . 


The 


final  result  in  Equation  (92]  can  be  found  in  integral  tables69.  Here 

2 


=  K2  R2  and  b2  =  K^. 


Air  drag  was  included  by  using  a  series  expansion  of  e 
term  by  term. 


and  integrating 


Weiss  also  considered  various  other  questions.  For  instance,  how 
are  we  to  combine  component  kill  probabilities  in  the  case  of  duplicated 
or  multiply  vulnerable  components?  A  difficulty  arises  for  example 
when  we  consider  two  pilots  sitting  close  together  so  that  a  round 
bursting  close  to  one  is  also  near  the  othet.  Because  of  this  lack  of 
position  independence,  the  probability  of  killing  pilots  computed  by 
averaging  over  all  burst  positions  will  not  be  the  same  as  that  obtained 
by  computing  the  probability  of  killing  each  one  independently  and  then 
combining  the  probabilities.  In  an  appendix  he  concluded  that  the 
approximation  of  treating  duplicated  components  as  if  they  were  independent 
is  good  enough  for  practical  purposes.  Weiss  also  considered  the  com¬ 
bined  mechanical  time  and  point  detonating  fuze  for  which  the  probability 
of  a  direct  hit  should  be  evaluated  and  then  combined  with  the  probability 
of  fragment  kill.  Finally,  he  discussed  the  correlation  between  succes¬ 
sive  burst  positions  when  systematic  errors  exceed  random  errors, 
especially  in  the  case  of  a  maneuvering  target. 

In  1951,  Taylor  and  Kravitz70  modified  the  work  of  Lewy  and  Gurney 
by  changing  their  integration  variables  to  two  angles  instead  of  an 
angle  and  a  distance.  In  addition,  they  approximated  the  damage 
probability  function  as 


instead  of  Equation  (35)  above.  Here  a  and  3  are  purely  empirical 
constants  obtained  by  fitting  over  an  appropriate  range  for  R.  This 
enabled  them  to  integrate  over  one  angle  in  terms  of  Bessel  functions, 
but  numerical  methods  were  required  to  integrate  over  the  second  angle. 


6yJ.  M.  Ryshik  and  I.  S.  Gradstein ,  Tables  of  Series,  Products  and 
Integrals,  Moscow  63. 

70 W.  C.  Taylor  and  S.  Kravitz ,  UA  Method  for  Computing  an  Estimate  of  the 
Antipersonnel  Effect  of  an  Airburst  Shell,  n  BRL  MR5S5,  13  Aug  51 . (AD#802074) 
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About  the  same  time,  Trauring71  in  an  appendix  elaborated  upon  Lewy  and 
Gurney's  numerical  method  of  evaluating  their  distance  integral. 

The  question  of  shaping  the  nose  of  a  fragmenting  bomb  or  shell 
received  some  attention  also  at  this  time.  Kent72  reasoned  that  if 
personnel  targets  could  be  considered  to  be  uniformly  distributed  on 
the  ground,  then  it  should  be  advantageous  to  shape  the  warhead  so  that 
a  uniform  distribution  of  fragments  on  the  ground  would  result.  Kent 
considered  the  case  of  a  bomb  falling  vertically  with  negligible  remaining 
velocity  and  bursting  at  height  h.  He  assumed  that  the  fragments  are 
projected  normal  to  the  metal  case  and  wrote  2tt  a  x.ds  for  the  number  of 
fragments  in  a  ring  of  width  ds  on  the  bomb  surface  located  a  distance 
*  from  the  axis,  where  a  is  the  constant  number  of  fragments  per  unit 
surface  area.  If  0  is  the  co- latitude  measured  from  the  bomb  axis  to 
the  radius  from  the  bomb  center  to  a  point  on  the  ground,  then  r  =  h 
tan  0  is  the  distance  on  the  ground  to  this  point  from  directly  below 
the  burst  point.  The  fragments  from  the  ring  of  width  ds2on  the  bomb 
fall  on  a  ground  ring  of  area  2irrdr  =  2tt  (h  tan  0)  (h  sec  0  d  0) ,  so 
the  number  of  fragments  per  unit  area  of  ground  is  the  ratio 

_ 2i  a  x  ds _ 

2  2 
2ir  h  tan  0  sec  0  d  0 


where  we  wish  Z  to  be  constant.  Since  ds  = 

2 

tan  0'  =  dy/dx  with  d(tan  9)  =  sec  0  d  0,  Equation  (87)  becomes 


x  ix .  (h2 (£)  d(£)  C9S) 

The  first  integral  is 

(£)-  4  (96) 

where  and  9  for  x  =  0.  The  second  integral  gives 

a  curve  through  the  origin 


=  Z  (94) 


r  2 

/l+ (dy/dx)  dx  and 


7  Trauring ,  "The  Effectiveness  of  Various  Weapons  Used  in  Air  Attack 
on  Ground  Troops , "  BRL7S4 3  17  May  SI. (AD#377188) 

72R.  H.  Kent3  "^ie  Shape  of  a  Fragmentation  Bomb  to  Produce  Uniform  Frag¬ 
ment  Densities  on  the  Grounds"  BRL  R7623  IS  Jun  51.  (AD# 377098) 
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(97) 


2  \  V2 
!-) 


a  solution  which  Kent  attributed  to  a  Professor  Barnett.  Since  x  and  y 
increase  without  limit  in  Equation  (97),  there  obviously  must  be  a  change 
in  shape,  perhaps  to  a  cylinder,  at  some  finite  x  value.  In  other  words. 
Equation  (97)  describes  only  a  nose  shape.  It  is  somewhat  of  a  curiosity 
since  it  does  not  consider  such  problems  as  non-vertical  angles  of  fall, 
the  presence  of  nose  fuzes  or  the  question  of  how  to  initiate  the 
explosive  charge  in  order  to  produce  the  assumed  normal  launch  direction 
for  the  fragments.  Since  aerodynamic  considerations  impose  some 
requirements  on  the  design  of  a  bomb  or  shell  nose  shape,  one  might 
also  inquire  about  the  flight  characteristics  provided  by  Equation  (97) 
as  well  as  the  effect  of  remaining  velocity.  Still,  the  whole  idea  of 
shaping  warheads  for  optimum  effect  embodied  in  Equation  (97)  stimulated 
a  number  of  people  to  pursue  the  topic.  Using  the  same  assumptions, 

Weiss 7 3  considered  more  general  methods  of  determining  optimum  angular 
distributions.  His  conclusion  was  that  the  number  of  fragments  per 
unit  solid  angle  should  increase  with  angle  measured  from  the  nose, 
reach  a  maximum  and  then  vanish  as  the  equator  is  approached.  Of 
course,  a  bomb  or  shell  with  a  main  equatorial  side  spray  thrown  forward 
by  a  large  remaining  velocity  will  have  this  type  of  dynamic  fragment 
distribution.  Other  workers  published  reports  considering  the  same 
problem7*1*  75 . 

In  1952  Weiss76  published  a  fairly  comprehensive  discussion  of  the 
methodology  in  use  at  the  time.  He  included  an  explicit  discussion  of 
the  addition  of  the  static  fragment  launch  velocity  v  ,  and  the  remain¬ 
ing  velocity  of  the  shell,  v  , to  obtain  the  resultant  launch  velocity  v, 
with  components 


v  sin  9  =  v  sin  8 
o  o 

v  cos  0  =  v  cos  0  +  v  (98) 

o  o  r  v 

7  iH.  K.  Weiss ,  "Optimum  Angular  Fragment  Distributions  for  Air-Ground 
Warheads ,"  BRL  R829y  Sep  52.  (AD#377182) 

74 "A  Theoretical  Basis  for  the  Optimization  of  Munition  Shape  and  Munition 
Distribution  Raving  a  Generalized  poisson  Formulation,"  Aug  62,  Honeywell 
Ordnance  Division,  Hopkins ,  Minnesota. 

75L.  F.  Nichols  and  W.  R.  Benson,  "Optimum  Fragment  Distributions  Using 
the  Calculus  of  Variations, "  Picatinny  Arsenal  TM121S,  Jun  63. 

7GH.  K.  Weiss,  "Methods  for  Computing  the  Effectiveness  of  Fragmentation 
Weapons  Against  Targets  on  the  Ground,"  BRL  R800,  Jan  52.  (AD# 397181) 
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from  which  we  can  obtain  the  dynamic  launch  speed 
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v  +  2  v  v 
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and  launch  angle 
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arc  tan 


cos 


0  + 
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(100) 


where  0  and  0  are  measured  from  the  shell  axis  which  is  taken  to  be 
parallel  to  v  . 

Since  effectiveness  depends  on  both  fragment  number  and  speed, 

Weiss  reasoned  that  for  a  given  total  weight,  W  =  M  +  C,  there  ought  to 
be  an  optimum  value  of  C/M  for  a  given  target  damage  criterion  since 
larger  C  favors  larger  v,  while  larger  M  favors  larger  N.  With 
Equation  (86)  as  damage  criterion  determining  m^,  the  number  greater 

than  this  minimum  mass,  Nm>  can  be  written  from  Equation  (85).  If  we 

divide  this  number  by  the  case  mass,  we  obtain  the  number  of  potentially 
lethal  fragments  per  unit  case  mass,  namely. 


since  the  average  mass  can  be  calculated  as 


(101) 


d  (%)  =n^-Xdx  =  v2^ 


(102) 


Equation  (101)  is  a  function  of  mQ  which  has  a  maximum  for 


ra  =  (2  m  )  =  m  /2  (103) 

o  m 

so  choosing  a  design  which  results  in  natural  fragmentation  characterized 

by  ra  =  ra  /4  in  Equation  (85)  will  maximize  the  number  of  lethal  frag- 
o  m 

ments  obtainable  per  unit  case  mass.  Choosing  an  average  mass  equal 
to  half  the  minimum  mass  as  in  Equation  (103)  for  natural  fragmentation 
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is  in  sharp  contrast  to  Mott's  choice  for  the  mass  of  a  controlled  or 
preformed  fragment  equal  to  4.48  as  expressed  in  Equation  (48) 

above.  If  we  rely  on  natural  fragmentation,  the  best  we  can  do  is 
make  the  average  mass  ineffectively  small  according  to  Equation  (103) . 
If  we  put  Equation  (103)  in  Equation  (101)  we  find 


-  v  Y 

(N  /M)w  =  (2/m  )e  =  .27/m  =  .27  -f- 

m  Max  Jtr  '  m  Lj 


(104) 


In  the  case  of  zero  air  drag  we  should  choose  a  controlled  fragment 

size  of  in  ,  so  that  (N  /N)  =  1/m  ,  illustrating  that  natural  fragmen- 
m  ram 

tation  can  only  be  27%  as  efficient  by  Equation  (104).  Of  course, 

if  we  took  air  drag  into  account  as  Mott  did,  we  would  choose  a  frag¬ 
ment  size  larger  still. 


For  zero  air  drag,  Weiss  went  on  to  assume  that  fragment  size 
optimization  had  been  carried  out  by  adopting  Equation  (103),  and 
discussed  the  maximum  number  of  potentially  lethal  fragments  as  a 
function  of  the  ratio  M/W.  In  other  words,  he  first  determined  an 
optimum  fragment  size  and  then  an  optimum  launch  speed,  v  . 


He  realized >that  the  problem  should  strictly  speaking  not  be 
treated  in  a  stepwise  fashion  since  all  parameters  should  be  optimized 
simultaneously.  However,  he  hoped  to  gain  some  insight  by  this 
simplified  procedure.  From  Gurney's  formula  for  a  cylinder  we  have 

^vq2/4E^  =  1/  ^+2  M/C^  and 


M  M  M/C 

W  “  C+M  1+M/C 


l-(vo2/4E) 
1+ (vq2/ 4E) 


(105) 


as  in  Equation  (81)  above. 


If  we  multiply  Equations  (104)  and  (105) 


CVH)Max  W»)  -  (.27  voVl^1  -(v02/4e)]/[i*(v02/4e)]  .  (106) 

For  the  dynamic  problem  we  should  use  Equation  (99)  above  for  the 
launch  speed  instead  of  the  static  launch  speed  v  .  If  we  neglect  v^ 

and  maximize  Equation  (106)  with  respect  to  vq  we  find  this  occurs  for 
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i  +  (2/y)2  -  (2/y) 


(107) 


2/4E^  =  V 


r 

which  i-s  0.414  for  y..=  .:2.. giving  ^  =  1.414  in  Equation  (105).  Similarly, 
y  =  3/2  gives  ^  =  1.000  and  y  =  1  gives  ^  =  .618.  Strictly  speaking, 
this  does  not  give  us  (N  /W)  because  of  the  stepwise  procedure. 

m  M3.X 

Corrections  for  air  drag  would  tend  to  raise  the  optimum  value  of 
average  mass  and  lower  the  optimum  value  of  launch  velocity. 


Weiss  also  discussed  the  desirability  of  using  a  function  other 
than  the  simple  step  function  implied  by  Equation  (86)  to  express  the 
conditional  kill  probability  of  a  target  which  has  been  hit,  Pk 

Equation  (86)  implies  that  for  L  =  m  v* 

L  >  =  1 •  1^  would  seem  more  reasonable  for  P. 

from  zero  to  unity  over  a  range  of  values  for  L. 
linear  relation 


<  Ll»  PK/H 


K/H“ 

=  0  and  for 

to  increase 


K/H 

Weiss  considered  a 
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(108) 


for  <  L  <  L2>  with  P^ 

If  we  use  N  =  N  e-  ^%/mo^ 
m  o 

obtain 


=  0  for  L  <_  L. 

1/2  1 

and  mi  = 
m 


and  =  1  for  L  >_  L^. 

_  y 

Vq  in  Equation  (90) , 


we 


N  =  Nq  exp  |  -  mo‘1/2  [lj1/2  vQ'y/2  +  ^ .45  aV  1/6vq‘Y/6)r]  j  (109) 

for  the  number  of  potentially  effective  fragments  as  a  function  of  R. 
Weiss  used  L  instead  of  Lj  in  Equation  (109)  and  argued  qualitatively 

that  the  appearance  of  L  to  fractional  powers  in  Equation  (109), 
especially  the  1/6  power  in  the  coefficient  of  R,  would  not  change 
the  R  dependence  of  N  or  N.  very  much  even  if  a  relation  like  Equation 
(108)  were  used.  Consequently,  he  convinced  himself  that  it  would  not 
be  worthwhile  to  consider  this  additional  complication  in  spite  of  the 
fact  that  it  seems  reasonable  to  do  so.  However,  he  gave  no  quantitative 
arguments  to  justify  this  position. 

It  seems  that  Weiss  found  himself  in  this  difficulty  because  of  his 
method  of  arriving  at  Equation  (90)  and  so  Equation  (91)  .  If  he  had 
adopted  the  viewpoint  of  0SRD  738  as  expressed  in  Equation  (28)  above 
in  which  decrease  in  effectiveness  due  to  air  drag  is  incorporated 
into  A  rather  than  N,  then  he  would  not  have  found  himself  in  this 
position.  We  can  rewrite  Equation  (28)  as 
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(110) 


NAy(0) 

2irR2 


K  W) 


PK/H^ 


where  the  target  kill  probability,  given  a  hit,  is  PV/L1(R)  =  A  (R)  / 

Av(0)  with  Ay(0)  <_  Ap.  Now  Equation  (108)  or  some  other  functional 

form  dependent  on  m  v*(R)  can  be  used  in  Equation  (110)  to  estimate 

-N 

Nl  and  the  damage  function  D  =  1  -  e  L  .  We  shall  see  an  example 
of  another  functional  form  for  in  Equation  (113)  below. 

In  a  later  part  of  his  report,  Weiss  suggested  another  approxima¬ 
tion  for  the  damage  function,  namely, 


-Nt 


D  =  1  -  e 


N. 


-CR  „ 
K  e  /R" 


(111) 


and  D  =  1  for  R  <  R^.  Here  is  Kent's  radius  of  overhitting  in  a 

three-dimensional  context  and  with  air  drag  and  target  vulnerability 
contained  in  the  argument  of  the  exponential  in  the  coefficient  G. 

This  approximation  not  only  unifies  the  idea  of  a  radius  of  overhitting 
with  the  poisson  expression  for  the  probability  of  at  least  one  lethal 
hit,  but  also  simplifies  many  integrals  involving  the  product  of  D 
with  other  factors  like  aiming  error  functions  or  ground  target  cover 
functions.  Weiss  evaluated  a  number  of  such  integrals  for  special  cases 
in  closed  form,  but  required  numerical  integration  for  general  cases. 

He  used  this  approximation  to  estimate  ground  lethal  areas  in  later 
reports77*"  .  Weiss  advocated  the  use  of  integrable  formulations  whenever 
possible  since  these  more  easily  contribute  to  our  general  understanding 
of  a  problem.  In  particular,  he  pointed  out  that  aim  error  could  be 
neglected  when  calculating  the  effectiveness  of  an  area  weapon  since 
it  was  smaller  than  the  extent  of  the  target.  However,  he  was  careful 
to  include  it  when  its  standard  deviation,  a,  was  comparable  to  or  larger 
than  either  the  extent  of  the  target  or  the  extent  of  the  weapon's 
lethal  mechanism.  His  approach  was  in  sharp  contrast  to  some  of  the 
purely  empirical  formulations  in  use  at  the  time79.  The  Weiss 

r,H.  K.  Weiss,  "Description  of  a  Lethal  Area  Computation  Problem ,  " 

BRL  MR72S ,  Sep  S3. (AD#21133) 

7BH.  K.  Weiss j  " Methods  for  Computing  the  Effectiveness  of  Area  Weapons , " 
BRL  R879 ,  Sep  S3.  (AD#24478) 

79L.  N.  Enequist,  "Rapid  Estimation  of  Lethal  Areas  for  Mortar  and 
Similar  Thin  Walled  H.  E.  Projectiles,",  BRL  MR635 ,  Dec  52.  (AD#3766) 
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approximation  in  Equations  (110)  and  (111)  can  be  compared  with  the 
Kendall-Carlton  approximation  in  Equation  (35)  and  the  Taylor-Kravitz 
approximation  in  Equation  (93) . 

For  many  years  various  authors  had  mentioned  the  desirability  of 
improving  the  threshold  lethality  criterion  for  disabling  personnel, 

1  2 

usually  taken  to  be  —  m  v  =58  ft-pounds.  In  1956  Allen  and 

Sperrazza80  reviewed  previous  criteria  and  proposed  a  new  criterion  based 
on  extensive  firing  tests.  They  adopted  a  purely  empirical  formula 


and  found  that  =  3/2  seemed  to  fit  the  data  best.  The  other  three 
positive  parameters,  a,  b  and  n,  assumed  various  values  (with  n  not  far 
from  0.5)  when  this  formula  was  used  to  describe  probability  of 
incapacitation  in  various  military  stress  situations.  Incapacitation 
within  five  minutes  for  a  defender  was  assigned  a  different  probability 
than  for  an  attacker,  given  a  random  hit.  If  this  formula  is  used  instead 
of  Equation  (108)  or  the  simpler  threshold  criterion,  then  it  is  clear 
that  simply  integrable  expressions  for  lethal  area  become  much  more 
difficult.  Since  Equation  (113)  is  a  purely  empirical  formula,  it  may 
be  desirable  to  find  another  formula  to  fit  the  data  which  allows 
integrations  to  be  carried  out  in  closed  form. 

In  the  next  decade  lethal  area  methodology  underwent  various  refine¬ 
ments,  but  no  basic  changes.  In  1963  Myers8*  codified  and  updated  the 
basic  methods  discussed  by  Weiss  to  include  more  complicated  cover 
functions82  and  the  effects  of  blast.  He  also  used  Equation  (113)  for 
his  casualty  criterion,  treating  each  mass  group  separately  and  summing 
over  all  groups  to  obtain  the  net  for  the  distribution  centered 

at  a  particular  polar  angle  0  (or  latitude  if>) .  In  addition,  he  treated 
fragment  angular  dependence  by  using  average  values  of  K(^)  at  five 
degree  intervals  (instead  of  2.5°  intervals  as  in  OSRD738) .  The 
effects  of  air  drag  were  estimated  through  v(R)  in  Equation  (113)  for 
which  he  used  the  usual  formula.  Because  of  these  complications  he 
always  had  to  use  numerical  integration  schemes.  This  1963  report  of 
Myers  has  become  the  standard  for  the  industry03  and  is  still  used 
Allen  and  J.  Sperrazza t  "Nett  Casualty  Criteria  for  Wounding  by 
Fragments , "  BRL  R996,  Oot  56.  (AD#137681) 

81tf.  A.  Myers ,  " Lethal  Area  Description 3"  BRL  TN1510 ,  Jul  63.  (AD#612041) 

82fl.  W.  Harris  and  K.  A.  Myers ,  " Cover  Functions  for  Prone  and  Standing 
Men  Targets  on  Various  Types  of  Terrain tn  BRL  MR1203 3  Mar  59.  (AD#309048) 

8 ^ Reference  2  in  Joint  Munitions  Effectiveness  Manual y  "Lethal  Areas  of 
Selected  U,  S.  Army,  U.  S.  Navy,  and  U.  S.  Marine  Corps  Surface  to 
Surface  Weapons  Against  Personnel  and  Materiel  Targets , "  Jul  73. 
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today  with  only  minor  modifications.  It  has  the  virtue  of  including 
in  some  way  or  other  all  of  the  factors  which  are  felt  to  be  important. 
However,  like  most  complicated  codifications  it  has  the  disadvantage 
of  obscuring  the  experimental  and  theoretical  approximations  which 
are  buried  in  the  code  and  lull  the  unwary  user  into  believing  that 
the  answers  are  very  accurate  or  that  minor  differences  in  the  answers 
are  significant.  Since  complexity  impedes  understanding,  it  is  not 
helpful  as  a  design  tool  either,  and  in  fact  is  generally  used  merely 
to  evaluate  designs  already  arrived  at  by  other  means.  Myers  himself 
seemed  to  realize  the  approximate  nature  of  his  description  and  discuss¬ 
ed  in  particular  the  use  of  the  poisson  approximation  for  the  damage 
function  instead  of  the  Binomial  which  some  people  felt  to  be  more 
correct.  As  he  pointed  out,  "the  Binomial  itself  must  be  considered 
only  as  an  approximation  since  fragments  projected  from  the  same  shell 
cannot  really  be  considered  as  'independent  events'.  He  added  that: 

"To  the  author's  knowledge,  no  extensive  comparisons  of  lethal 
areas  computed  using  the  Binomial  and  poisson  forms  for  the  probability 
of  incapacitation  have  been  made.  However,  in  some  isolated  cases 
where  such  comparisons  have  been  made,  the  difference  did  not  appear 
to  be  significant,  particularly  in  light  of  other  assumptions  inherent 
in  the  lethal  areas  concept".8^ 

B.  The  Adoption  of  Monte  Carlo  Methods  for  Air  Defense  Problems. 

In  order  to  treat  the  complexities  of  the  air  defense  problem,  a 
different  (Monte  Carlo)  technique  has  been  adopted  instead  of  the 
analytical  techniques  which  were  used  until  the  early  1950's.  In 
1950,  F.  G.  King84  described  his  method.  Although  he  mentions  a 
previous  partial  formulation  of  the  problem  by  Cunningham  in  England, 
his  report  is  clearly  the  ancestor  of  the  end-game  analyses  which  are 
still  in  use  today.  Although  today's  analyses  are  very  complicated 
and  detailed.  King's  original  reason  for  using  this  method  was 
simplicity: 

"The  lotto  method  proposed  in  the  present  report  seems  very  simple 
when  compared  with  the  formulation  of  the  problem  in  mathematical  sym¬ 
bolism.  The  basic  vulnerability  data  is  known  for  the  most  part  only 
to  order  of  magnitude  and  at  best  only  to  ten  or  twenty  per  cent,  so 
that  a  simple  method  is  in  order.  If  the  basic  data  were  known  within 
one  per  cent,  the  lotto  method  would  still  be  justified  as  a  short  cut 
to  an  approximate  answer  but  would  be  inefficient  for  getting  the 
answer  to  one  per  cent  accuracy." 

King's  lotto  or  monte  carlo  method  is  a  mathematical  experiment  in 
which  the  events  of  an  anti  aircraft  attack  are  acted  out.  First,  for 
a  given  engagement  range  a  burst  is  positioned  with  respect  to  the 
target  by  drawing  a  card  at  random  from  a  box  which  represents  the 
distribution  of  aiming  errors  at  the  range  in  question.  King  used 

G.  King,  " Lotto  Method  of  Computing  Kill  Probability  of  Large 
Warheads , «  BEL  MRSSO ,  Dec  50.  (AD#802148) 
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one  thousand  cards  as  a  discrete  representation  of  a  trivariate  normal 
distribution  about  the  aim  point.  Second,  the  distance  and  angle 
from  the  burst  point  to  each  vulnerable  component  was  measured,  using 
a  model  airplane  on  a  target  stand  and  a  model  burst  on  a  warhead 
stand  to  a  scale  of  1:72.  A  component  was  taken  to  be  shielded  if  at 
least  its  center  of  mass  was  shielded.  Third,  using  approximate  treat¬ 
ments  of  missile  and  target  velocity,  air  drag  at  a  given  altitude, 
warhead  fragmentation  characteristics  and  component  vulnerability  data, 
the  kill  probability  of  each  vital  component  was  estimated  as  a  function 
of  miss  distance  (and  angle  if  angular  data  was  available).  Unfortun¬ 
ately  King  gives  no  detailed  description  of  how  this  crucial  step  was 
actually  carried  out,  but  it  is  safe  to  assume  that  it  involved  the 
techniques  we  have  discussed  previously.  Fourth,  the  kill  probability 
of  each  component  computed  in  step  three  for  a  particular  miss  distance 
was  turned  into  a  zero  or  a  one  by  reading  down  a  list  of  random  numbers 
between  zero  and  one.  If  the  next  random  number  read  was  less  than 
the  computed  kill  probability  the  component  was  considered  killed, 
otherwise  not.  Sometimes  this  step  was  repeated  several  times  and  the 
results  averaged  before  going  on.  Fifth,  for  each  burst  it  was  decided 
whether  or  not  the  aircraft  was  killed,  subject  to  the  requirement 
that  a  certain  number  of  engines  or  pilots  must  be  killed  or  that  fuel 
be  hit  or  that  the  burst  occur  within  the  blast  kill  envelope.  Sixth, 
steps  one  through  five  were  repeated  several  hundred  times.  Seventh, 
the  single  shot  kill  probability  at  the  given  engagement  range  was 
found  by  averaging  the  single  burst  kill  probabilities  obtained  in 
each  step  six.  A  two-shot  kill  probability  could  be  determined  by 
taking  an  average  over  a  number  of  pairs  of  bursts,  considering  cumula¬ 
tive  damage,  and  so  on.  Eighth,  since  engagements  between  an  aircraft 
and  a  rapid  fire  gun  or  battery  of  guns  consist  of  single  or  multiple 
shots  at  successively  different  ranges,  a  new  box  of  cards  was  introduced 
in  step  one  and  steps  one  through  seven  were  repeated  for  a  number  of 
representative  ranges.  This  eighth  step  can  also  be  repeated  an 
appropriate  number  of-times  to  improve  the  accuracy  of  the  prediction. 
However,  there  is  no  point  to  repeating  it  too  many  times  since  overall 
accuracy  is  already  limited  by  the  accuracy  of  the  experimental  frag¬ 
mentation  and  vulnerability  data.  An  application  of  part  of  this 
method  to  a  particular  system  was  later  made  by  Sacks  and  King85. 

In  1953  Juncosa  and  Young86  computerized  King's  method  after  replacing 
parts  of  it  with  mathematical  idealizations.  In  particular,  they  adopted 
a  suggestion  which  they  attributed  to  J.  Von  Neumann  that  the  airplane 
and  its  blast  kill  envelope  be  replaced  by  ellipsoids  resembling  the 

at> S.  Sacks  and  F.  G.  King,  "Single-shot  'Probability  of  an  A-Kill  on  the 
B-29  Type  Bomber  by  Nike  I,"  BRL  MR  55Z,  Jul  51.  (AD#377282) 

86M.  L.  Juncosa  and  D.  M.  Young ,  "A  Mathematical  Formulation  for  Ordvac 
Computation  of  the  Probability  of  Kill  of  an  Airplane  by  a  Missile 
BRL  R867j  May  52.  (AD#17267) 
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the  fuselage,  wings,  empennage  and  engines.  In  addition,  they  replaced 
all  vulnerable  components  by  points  to  simplify  shielding  by  the  air¬ 
craft  parts  represented  by  ellipsoids.  Third,  they  restricted  their  con¬ 
siderations  to  a  particular  angle  of  attack,  and  finally,  they  evaluated 
three  dimensional  integrals  over  the  product  of  aim  and  damage  functions 
by  a  monte  carlo  method  or  by  numerical  quadratures  instead  of  averaging 
over  the  results  of  several  hundred  randomly  chosen  burst  points. 

In  1960,  Stiegler87  modified  the  method  of  Juncosa  and  Young.  He 
allowed  random  (monte  carlo)  selection  of  straight  line  approaches  by 
the  missile  toward  the  target  for  elevation  angles  other  than  zero  degrees 
and  azimuths  other  than  45  degrees.  He  also  computed  a  missile  homing 
point  which  is  the  center  of  the  target  "glitter  points"  (radar  .  reflect¬ 
ing  points)  rather  than  the  center  of  the  aircraft.  He  also  computed  the 
intersection  of  "glitter  points"  with  the  fuze  cone  and  took  direct  hits 
into  account.  Stiegler' s  method  was  later  modified  by  Monks88-90  and 
various  versions  of  this  technique  are  still  in  use  today. 

As  can  be  seen,  the  development  of  the  monte  carlo  method  since 
1950  has  followed  the  same  trend  as  the  development  of  lethal  area 
methodology  and  similar  remarks  can  be  made.  More  and  more  detail  has 
been  included  so  that  no  feature  be  left  out  of  the  model.  Unfortunately, 
the  resulting  complicated  codes  and  the  precision  to  which  the  final 
answer  is  expressed  tend  to  make  the  unwary  user  place  more  credence  in 
the  answers  than  is  justified.  By  and  large  the  methodology  cannot  be 
used  as  a  design  tool,  and  is  only  used  to  evaluate  weapon  designs 
arrived  at  in  other  ways. 

Although  monte  carlo  methods  have  dominated  the  evaluation  of  air 
defense  effectiveness,  there  have  been  some  reports  on  the  development 
of  analytical  methods  in  this  area.  For  example,  in  1966, 

Banash91  proposed  using  the  Kendall -Carlton  approximation  in  Equation  (40) 

above  with  R^2  =  1.25  K  so 

a/i4.  D.  Stiegler,  "A  Mathematical  Formulation  for  ORDVAC  Computation 
of  the  Single-Shot  Kill  Probabilities  of  a  General  Missile  Versus  a 
General  Aircraft /  BRL  MR1306 ,  Nov  60.  (AD#249957) 

88U.  G.  Monks,  "A  Mathematical  Simulation  of  the  SA-2  Missile  for  Evalua¬ 
tion  of  its  Lethality /  BRL-J4R1690 ,  Aug  65 .  (AD#368987L) 

89U.  G.  Monks ,  w Aberdeen  Blast-Fragmentation ,  Fixed  Angle  Fuze  Missile 
End  Game  Simulation  (AMEGS), "  AMSAA  TR66 ,  Apr  73. 

9uU.  G.  Monks,  "A  Focused  Blast-Fragmentation  Warhead  Evaluation  Model," 
AMSAA  TR70 ,  Feb  73. 

9 C.  Banash,  " An  Analytic  Procedure  for  Evaluating  Blast-Fragmenta¬ 
tion  Warheads  in  Air  Defense,"  BRL  R1317 ,  Apr  66.  (AD#485264) 
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(114) 


k/r2  Rl_  r 

D  =  1  -  e  VK  *  e  „  2  =  e  2.5K 

2R3 


where  air  drag  is  being  neglected.  Equation  (114)  is  a  fair  approxi¬ 
mation  for  miss  distances  in  the  range  0  <  K  <  and  eventually 
improves  again  as  R  ■>  00  and  D  ->  0 .  He  compared  kill  probabilities  cal¬ 
culated  by  his  analytical  method  with  those  obtained  by  the  monte  carlo 
method  in  four  examples  involving  three  hypothetical  warheads  versus 
missile  and  aircraft  targets.  In  most  cases  the  agreement  was  within  6%. 

C.  Current  Vulnerability  Methodology.  Warhead  Optimization. 

Whatever  the  method  used  to  estimate  effectiveness,  there  is  a 
continuing  need  for  vulnerable  area  information  or  kill  probability, 
given  a  hit,  since  targets  and  weapons  are  constantly  changing.  For 
aircraft  the  current  state  of  the  art  has  been  summarized  in  a  recent 
report92.  In  Section  II  above  we  quoted  professor  Mott22  who  remarked 
on  the  wide  discrepancy  between  combat  data  and  the  effectiveness  esti¬ 
mates  of  his  day.  He  would  be  encouraged  to  hear  that  matters  have 
improved  somewhat.  Combat  damage  data  has  been  used  to  develop  vul¬ 
nerable  areas  for  at  least  three  helicopters. 


"In  all  three  cases,  the  average  vulnerable  areas  obtained  from 
the  combat  damage  data  and  generated  in  vulnerability  studies  are 
consistent."93 


While  this  agreement  is  not  for  overall  effectiveness,  still  it  is 
encouraging  in  that  combat  data  can  be  used  to  provide  guidance  where- 
ever  possible.  As  these  authors  remark,  "Data  on  aircraft  damage 
which  occurs  in  combat  are  very  valuable  since  they  originate  in  the 
'real  world'  conditions  of  military  conflict....  thus  the  body  of  air¬ 
craft  combat  damage  data  provided  by  the  collections  from  Southeast 
Asia  serves  as  the  main  source  of  information  on  aircraft  reactions  to 
the  effects  of  component  damage  while  airborne.”  They  also  remark  on 
use  of  such  data  for  vulnerability  reduction.  Still  in  many  cases  not 
only  is  combat  damage  data  not  available,  but  not  even  controlled 
experimental  information  on  specific  threat  performance  or  target 
vulnerability  is  available. 


y2#.  W.  Mowrer3  R.  E.  Walther ,  R.  D.  Mayerhofer  and  R.  N.  Schumacher 3 
"Aircraft  Vulnerability  Assessment  Methodology.  Volume  I  -  Generalj" 

BRL  R1796j  Jul  75.  See  also  Volumes  II  to  XII  for  particular 

subsystems.  Volume  I  also  appeared  as  JTCG/AS-76-V-0043  Jul  77. (AD#B005878L) 

9 ^Reference  92  3  p  72. 
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’’The  availability  and  validity  of  these  data  greatly  influence  the 
reliability  of  the  results  of  any  study  based  on  these  data.  It  is 
generally  difficult  to  quantify  this  influence  in  terms  of  confidence 
limits,  but  some  general  comments  are  germane.  A  completely  sufficient 
body  of  experimental  data  is  seldom,  if  ever  available  for  use  in  a 
study....  In  fact,  in  most  studies  these  considerations  require  that 
the  analyst  make  a  number  of  assumptions  to  fill  the  information  gaps. 
Thus,  certain  errors  may  be  introduced  into  the  study  by  the  necessity 
of  making  these  assumptions.  Furthermore,  all  of  the  analysis  proce¬ 
dures  require  certain  simplifying  assumptions  and  approximations,  which 
introduce  additional  errors...." 


"The  sensitivity  of  the  complex  models  used  in  vulnerability  studies 
(particularly  in  computer  programs)  to  the  assumptions,  approximations 
and  various  data  inputs  is  neither  apparent  nor  readily  determined."^ 


Still,  some  few  sensitivity  studies  have  been  carried  out  in  related 
areas  and  according  to  these  authors  the  final  answers  are  highly 
sensitive  to  all  of  the  parameters  examined,  namely,  fragment  residual 
weight,  speed  and  shape,  shotline  location  and  grid  size,  and  component 
detail  and  values.  Other  factors  obviously  introducing  error 

whose  magnitude  has  not  yet  been  examined  were  also  discussed.  Some 
of  these  are  neglect  of  all  but  the  largest  residual  fragment  pieces 
as  well  as  target  pieces  in  subsequent  perforations,  neglect  of  ricochets, 
the  need  to  use  single  plate  results  repetitively  to  describe 
multiplate  perforations,  the  need  to  estimate  equivalent  thicknesses 
of  non-plate  components  and  the  necessity  of  considering  only  a  small 
number  of  target  aspects. 


This  report  does  not  give  details  on  how  component  PK/H  values 

are  estimated.  Often  such  information  is  kept  classified.  However, 
some  unclassified  reports  discuss  the  matter  in  a  general  way.  For 
example,  Johnson95  points  out  that  for  a  given  fragment  mass  striking 
a  particular  component  at  a  certain  angle  the  kill  probability  is 
represented  by  a  multi-step  function  of  versus  striking  speed,  v. 

There  is  a  threshold  speed  below  which  PK/H  =  0  and  the  various  steps 
reflect  the  cumulative  damage  of  sub  components.  However,  PK/H  does 

not  approach  unity  no  matter  how  high  the  speed  in  Figure  3.  This 
seems  to  reflect  the  definition  used,  namely,  PK^  =  \/&p>  the 

vulnerable  portion  of  the  area  divided  by  the  measured  presented  area. 


^ Reference  923  pp  80-81. 

95V.  P.  Johnson 3  "An  Application  of  the  Weibull-Gnedenko  Distribution 
Function  for  Generalizing  Conditional  Kill  Probabilities  of  Single 
Fragment  Impacts  on  Target  Components , "  BRL  R17893  May  75.  (AD#B004591L) 
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STRIKING  VELOCITY  x  10‘2(m/s) 

Figure  3.  Typical  Kill  Probability  versus  Striking  Speed  Curve  for  a  Giyen 
Fragment  Mas s  and  Kill  Criterion  (from  BRL  R1789") 


GIVEN:  FRAGMENT  MAS5  AND  KILL  CRITERION 


with  <  Ap  because  parts  of  the  component  are  non-critical  dead  space. 
This  literal  interpretation  of  Ay  as  an  area  whose  maximum  value  is  the 
sum  of  a  number  of  critical  areas  less  than  Ap  can  lead  to  unnecessary 


semantic  difficulties.  Since  is  not  allowed  to  approach  unity  in 

this  interpretation,  it  can  be  objected  that  it  should  not  be  called 
a  probability.  This  is  a  valid  objection  since  one  requirement  for 
a  probability  distribution  is  that  the  sum  of  the  probablities  of  all 
its  subsets  be  equal  to  unity.  The  objection  can  be  avoided  by  calling 
p  a  damage  estimator  or  something  similar  rather  than  a  probability. 
A^simple  way  of  avoiding  the  problem  is  to  deal  only  with  Ay  as  in 

the  first  form  of  Equation  (28)  above.  We  can  of  course  multiply 


A  by  any  form  of  unity  we  desire  without  changing  the  result.  We 
cXuld  multiply  by  Ap/Ap  with  (Av)Max  <  Ap,  but  the  ratio  Ay/Ap  could 

not  be  called  a  probability.  Instead  we  could  meet  the  objection  by 
multiplying  by  some  A/A  such  that  A  =  (Av^Max  t0  obta^n  a  Probability 

Ay/A.  In  Equation  (28)  A  was  chosen  to  be  Ay(0)  so  that  Ay(R)/Ay(0) 

in  Equation  (29)  would  be  unity  for  R  =  0.  As  we  recall,  the  assump¬ 
tion  was  that  only  fragments  capable  of  damaging  the  target  at  short 
range  (R  =  0)  would  be  considered.  Since  R  =  0  corresponds  to  v  =  vq  by 

Equation  (27),  we  can  write  for  the  target  conditional  kill  probability 


PK/H  =  VR)/Av(0)  -  (11S) 

with  subscripts  added  if  components  are  being  considered.  Another 
alternative  would  be  to  let  A  be  equal  to  the  sum  of  the  presented 
areas  of  the  vulnerable  subcomponents  (in  the  case  of  a  component)  or 
the  sum  of  the  presented  areas  of  all  such  subcomponents  (for  the 
whole  target).  In  this  way  parts  of  the  target  which  are  not  even 
potentially  vulnerable  to  the  weapon  in  question  are  excluded  from  con¬ 
sideration.  Thus  the  difficulty  introduced  by  using  an  Ap  which 

includes  both  potentially  vulnerable  and  invulnerable  parts  can  be 
eliminated. 

Johnson  suggests  using  a  function  of  the  form 


PK/H  1 


exp  [-a  (mv  -  b)n] 


(.117) 


although  he  multiplies  it  by  a  factor  (Pj^/^)Max  <  1  since  he  accepts 
Figure  3  as  a  probability.  Here  we  have  let  (pKj/^Max  ”  1  so  Pk/h  oan  l 
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be  a  probability.  Equation  (117)  as  it  stands  is  the  well-known 
Weibull- Gnedenko  distribution  and  would  be  the  same  as  Equation  (113) 
above  if  the  variable  were  mv^  instead  of  mv.  There  are  analytical 
advantages  to  using  a  form  like  Equations  (117)  or  (113)  instead  of 
a  multitude  of  graphs  like  Figure  3.  In  addition,  it  can  help  to 
unify  our  description  of  both  personnel  and  materiel  vulnerability. 

Another  functional  form  has  been  proposed  by  Capillo96,  namely. 


1  +  BpCV  ^  ” 1  (.118) 

where  A,  B,  p  and  c  are  constants.  In  this  formulation  a  separate  set 
of  constants  is  required  for  each  mass  (and  angle)  so  it  is  not  as 
succinct  as  Equation  (117) .  In  addition,  it  does  not  resemble 
Equation  (113) .  Later  in  his  report  Capillo  gives  a  rather  succinct 
verbal  summary  of  the  state  of  the  art: 

"Vulnerability  assessment,  as  it  now  stands,  is  an  empirically 
based  quasi-science  drawing  knowledge  and  interpretation  from  historical 
data,  experimental  tests,  limited  analytics,  and  engineering  experience 
and  judgment." 

In  some  of  the  reports  w.e  haye  considered  s.o  far  w.e  have,  seen  that 
one  or  another  parameter  like  burst  height  or  fragment  mass  has  heen 
selected  and  varied  in  order  to  optimize  effectiveness.  This  tradition 
has  continued,  although,  when  numerical  methods  are  used,  the  range  of 
the  variable  selected  is  usually  rather  narrow  and  only  a  few  curves 
over  this  range  are  compared.  Occasionally  a  more  general  treatment 
using  analytical  methods  has  appeared.  In  1950  Dougherty97  discussed 
optimizing  the  fragment  angular  distribution,  K(ifi),  and  concluded 
that  a  cone  pattern  was  desirable  for  a  proximity  fuze.  He  also  gave 
a  derivation  of  the  damage  function  which  started  with  certain  assump¬ 
tions  and  employed  summation  techniques  instead  of  limiting  processes 
as  in  Equations  (12)  to  (18)  above.  He  let 


Nl  =  NP  (119) 

where  p  =  Av(R)/Av(0)  =  PK/H  and  N  =  NK(j|i)  Av(0)/(2frR2)  is  the 
expected  number  of  hits  so  that  is  the  expected  number  of  lethal 


3b6'.  C.  Capillo,  "Fighter  Aircraft  Vulnerable  Area  Data  Patterns: 
Modeling .  Data  Reduction ,  and  Interpolation  Schemes. "  BRL  MR2648.  Aug  76. 

(AD# BO 1334 8 L) 

37 C.  B.  Dougherty ,  "Optimization  of  Warhead  and  Fusing  Parameters , " 


Rand  Corporation3  RM349t  Mar  50. 
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hits.  He  then  assumed  that  the  chance  of  a  hit  follows  a  poisson 
distribution  with  expected  value  N,  while  the  chance  of  a  kill  obeys 
a  binomial  distribution 


G  (y;  n) 


(?)  Py  o-P)n->’ 


(120) 


which  is  the  chance  of  exactly  y  fragments  killing  the  target  out  of 
n  hits  (or  trials).  The  chance  of  at  least  one  fragment  killing  the 
target  is  then 


n  n 

£  G(y;n)  =  £  G(y;n}  -  G(p;n)  =  1-  (l-p)n  (121) 

y=i  .  y=o 

n 

since  G(y;n)  =  1.  If  the  hit  probability  is  assumed  to  follow  a 

y^o 

poisson  distribution,  then  the  joint  probability  of  at  least  one 
fragment  hitting  the, target  and  at  least  one  of  these  killing  it  is 


(122) 


which  is  the  usual  damage  function,  but  arrived  at  in  a  way  which 
brings  out  better  its  nature  as  a  joint  probability. 

It  may  be  objected  that  the  hit  probability  also  would  be 
described  better  by  a  binomial  distribution  especially  if  the  number 
of  fragments,  N,  is  not  very  large.  In  this  case  we  can  modify 
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Dougherty's  derivation  by  letting  the  probability  that  one  fragment 
hit  the  target  be 


*h=K0I0  Av(0)  /  (2ttR^ 

such  that  N  =  Np^  as  before.  Now  the  joint  probability  is 

d = (i  ■  p»)  N'n  [i= (i-p)n] 

*  (1-Ph)N  -  £  (n)  PH"  (1-Ph)N'"  (1-p)n 


(.123) 


=  1- 


Hit  era '-■] 


=  1- 


(i-e»y 


i  + 


d-V 


(124) 


*  l'  (1-PPh)N  -  1  -  (1-PK/H  " 

N 

where  we  have  used  the  identity  (1+X)N  =  £  (  jXn  in  the 

n=o  \n  / 

last  line,  then  rearranged  the  factor  in  brackets  amd  cancelled  (1-P,,)1 

5  H 

An  approximation  to  Equation  (124)  is 


next  to 


N  _  x_e  N  In  (1-PPH) 

*  1  -  exp  [-NPHP]  =  1  -  exp  [-NP]  =  l-e"NL  (125) 

which  is  the  usual  damage  function.  The  second  form  in  Equation  (125) 
is  merely  an  identity  while  the  third  form  is  an  approximation  for 
P^P  <<  1 .  As  we  see,  the  result  is  equivalent  to  using  the  poisson 
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approximation  to  the  binomial  for  the  hit  probability  in  Equation  (122). 
(Recall  Equation  (15)  above  also) .  Now  however,  we  see  that 


P 


H 


P 


P  P 

H  K/H  3 


KOJO  Av(o) 
2nR2 


A  (R) 
vv  J 

A  (o) 


«  1 


(126) 


is  close  to  unity)  or  by  having 
Thus  even  burst  close  to  a  tough  target  can  be  reasonably 


can  be  satisfied  by  having  PH  small  (even  if  P^H  is  close  to  unity) 
or  by  having  small  (even  if 

both  small, 

represented  by  the  approximation  in  Equation  (125) .  A  burst  close  to 
a  weaker  target  might  require  Equation  (124).  We  must  remember,  of 
course,  that  any  of  these  representations  are  model  idealizations  and 
presume  event  independence.  They  also  presume  a  model  geometry  which 
treats  the  burst  as  a  point  as  in  Equation (2)  above.  Thus  even  if 


2  2 

A  (R)  3  Ap  =  4tt  sin  (0/2)  R  ,  we  see  that  its  largest  value  is 

v  *  2 

A  (R)  =  2ttR  ,  since  0.,  =  tt/2  .  That  is,  the  burst  is  exterior  to 

v  Max 

the  target.  More  often  AV(R)  <  Ap. 


In  1964  Sewell98  presented  an  interesting  discussion  of  optimizing 
the  charge-to-metal  ratio,  x  =  C/M,  subject  to  weight  or  volume  con¬ 
straints.  He  did  not  attempt  to  maximize  kill  probabilities,  lethal 
areas  or  other  measures  of  effectiveness.  Instead,  he  looked  only 
at  the  warhead,  not  at  the  target  or  engagement  conditions  and 
maximized  the  launch  function,  E  =  M  voY,  for  y  =  1,  3/2  and  2.  In 

other  words  he  postulated  that  a  first  cut  in  warhead  design  might  be 
made  by  considering  the  maximum  total  momentum,  energy  (or  something 
in  between*)  which  could  be  imparted  to  the  warhead  case.  Since  M  = 
M(x) ,  then 


dE  _  3E  dM  9E_  =  0 

dx  3M  dx  +  0v  dx 


(127) 


o 

will  give  us  a  condition  for  a  maximum  of  E  with  respect  to  x  =  C/M. 
For  example,  for  a  fixed  weight  W  =  C  +  M  =  M(l+x),  so  M  =  W/(l+x)  and 


57  =  -  W  /  (1+x) 2  (128) 

We  may  rewrite  Gurney's  cylinder  formula,  the  analog  of  Equation  (59) 
above,  as  vq  =  (2E)  '  (1+ .  5x )  "  V2  xl/2  ^  s0 


98i?.  G.  S.  Sewell  j  "Fixed-Weight  and  Fixed-Volume  Constraints  on 
Optimum  Charge -to -Metal  Ratios  in  Warhead  Vesign} "  ROTS  TP  3430 , 
Mar  64. 
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.5  (2E) 


(129) 


dv 

o 

dx 


i/2  x“ 1/2  (1+-5X)-3/2 


and  for  E  =  M  vY  we  have  =  v  Y,  while  =  y  Mv  Y_1.  If  we 

o  3.M  o  3  v  o 

put  these  relations  in  Equation  (127)  we  find 


x  =  C/M  =V(y/2)2+1  +  (y/2  -  1)  (130) 

as  the  optimum  value.  Thus,  if  y  =  2,  x  =1.414  and  if  y  =  3/2,  x  = 

1.000  while  if  y  =  1,  x  =  ,618  exactly  as  in  Equations  (105)  to  (107) 
above.  It  is  not  hard  to  show  that  putting  Equation  (107)  into 
Equation  (105)  and  solving  for  C/M  yields  Equation  (130),  so  the 
equivalence  holds  for  every  Y.  This  is  not  surprising  since  optimiza¬ 
tion  of  E  with  respect  to  M/W  (or  C/M)  is  theoretically  independent 
of  fragment  size  in  this  view. 

Sewell  repeated  the  same  procedure  for  a. sphere  of  fixed  weight, 
finding  somewhat  smaller  values  of  C/M  to  be  optimum.  In  addition,  he 
discussed  both  geometries  subject  to  a  fixed  volume  constraint,  namely, 
constant  V  =  C/cg  +  M/pm,  for  which  the  optimum  C/M  values  will  depend 
on  the  ratio  of  explosive  and  metal  densities,  PE/Pm*  In  1976 
Zulkoski"  modified  Sewell's  treatment  by  postulating  various  factors 
to  represent  the  experimental  fact  that  cylinders  of  finite  length  do 
not  obey  the  Gurney  formula  except  at  one  cross  section,  say  at  the 
middle.  Fragments  launched  from  sections  near  the  ends  will  have 
different  (lower)  speeds  because  of  differences  in  explosive  gas 
product  confinement.  His  approach  otherwise  was  that  of  Sewell. 

Perhaps  the  most  comprehensive  attempt  to  optimize  warhead  charac¬ 
teristics,  at  least  in  the  case  of  area  targets  on  the  ground,  is  a 
1967  report  by  Scherich  and  Kitchen100.  The  parameters  they  considered 
were  warhead  height  of  burst,  angle  of  fall  and  remaining  velocity  as 
well  as  upper  and  lower  limits  of  the  side  spray  and  up  to  ten  fragment 
classes.  They  used  a  numerical  integration  scheme  to  evaluate  the 
lethal  area  integral  and  the  well-known  method  of  steepest  ascent  to 
seek  out  at  least  a  local  maximum.  To  avoid  missing  larger  values  of 
the  lethal  area  they  also  made  provision  for  searches  over  wider  ranges 
of  the  parameters.  However,  because  of  the  complexity  of  the  method 

^y!F.  Z ulkowskiy  "Development  of  Optimum  Theoretical  Warhead  Design 
Criteria , "  NWC  TP5892 ,  Aug  76. 

100E.  L.  Scherich  and  T.  Kitchen ,  "Weapon  Optimization  Techniques , " 
Martin  Mariettay  AFATL-TR-67-l28y  Oct  67. 
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used,  it  appears  to  be  difficult  to  gain  insight  into  the  optimization 
process. 


IV.  CONCLUSION 

In  this  report  we  have  chosen  to  view  our  subject  from  an  histori¬ 
cal  perspective  because  of  a  conviction  that  much  can  be  learned  in 
this  way.  Hopefully,  the  repetition  or  continuation  of  mistakes  and 
dead  ends  can  be  avoided,  while  hints  or  suggestions  for  new  ;  directions 
can  be  gathered  for: the  future.  In  addition,  the  meaning  of  what  we 
are  now  doing  might  be  clarified  by  seeing  how  it  arose  from  the  uni¬ 
fication  of  a  number  of  different  viewpoints.  Unfortunately,  we  have 
only  seen  the  early  British  and  the  American  Army  (BRL)  views  of  the 
subject.  More  could  be  learned  from  a  broader  perspective,  but  this 
report  is  already  long  enough.  A  certain  amount  of  new  material  has 
also  been  woven  into  the  text  to  assist  in  clarifying  our  ideas  and  to 
lay  the  groundwork  for  future  advances.  As  was  mentioned  in  the 
introduction,  a  later  report  in  this  series  will  consider  methods 
of  simplifying  descriptions  of  effectiveness  so  that  analytical  tools 
for  optimizing  warhead  designs  will  be  more  accessible.  Still  another 
report  will  be  devoted  to  the  origins  and  possible  improvement  of  the 
methodology  which  has  evolved  to  describe  the  effectiveness  of  anti¬ 
armor  weapons.  Since  direct  hits  by  single  lethal  mechanisms  (massive 
penetrators)  are  generally  required,  the  methodology  should  be  simpler 
in  some  ways.  However,  since  direct  hits  are  required,  we  must  always 
include  aiming  errors  in  our  overall  effectiveness  analysis.  Lethal 
areas  will  not  do.'  In  any  case,  there  are  enough  similarities  to 
include  such  a  discussion  in  the  present  report  series,  but  enough 
differences  to  warrant  a  separate  report. 
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APPENDIX  II.  FRAGMENT  PRESENTATION  AREA  AND  SHAPE  FACTOR 


In  their  paper16  mentioned  above,  Pearson  and  Bishop  discussed 
the  results  of  gun  firings  of  various  sizes  and  shapes  of  simulated 
fragments  (square  plates  ranging  from  .02  oz  to  .50  oz) .  They  observed 
that  the  loss  in  speed  with  distance  was  approximately  a  linear 
function  of  speed,  namely. 


dv/dR 


-  $  v 


(A-l) 


so  that 


dv/dt 


-  3  v 


(A-2) 


where  3  =  -  Cn  p  A/m,  if  Equation  (25)  above  holds.  However,  the 

^  a  -1/3 

experimental  data  indicated  that  6  was  proportional  to  m  instead 

of  m  1  as  might  b£  expected  if  A  were  independent  of  mass,  m.  For 
this  to  be  true,  A  would  have  to  be  proportional  to  m 2/3  as  in  Equation 
(26)  above.  This  is  clearly  true  in  the  case  of  a  sphere  for  which 

m  =  p  (^r)  d3where  p  is  its  density  and  d  is  its  diameter.  The  presented 

area  of  a  sphere  is,  of  course,  A  =  ^d2  =  ~  (— )  2/^3  m2/^3  which  is 
2/3  4  4  V  3 

proportional  to  m  .  For  a  cube  in  face-on  flight,  m  =  pd  and  A  = 

,2  =  (— )2//3  m2^3,  while  for  edge-on  flight,  A  =  (— )2/^3  m2^3.  In  1940 

dp  p 

Taylor101  remarked  that  rectangular  blocks  of  dimensions  a,  b  and  c, 

when  fired  from  guns,  would  probably  rotate  so  that  all  aspects  of 

presented  area  relative  to  the  flight  path  would  occur  with  about 

equal  frequency.  If  such  a  rectangular  parallelepiped  tumbled  randomly 

in  flight,  the  average  presentation  area  would  be 


A  = 


(ab  +  be  +  ca) 


(ab  +  be  +  ca)  2/3 
2p2/3  (abc )2/3  ” 


k  m 


2/3 


(A-3) 


when  put  in  the  form  of  Equation  (26) .  Pearson  and  Bishop  then  post¬ 
ulated  that  arbitrary  shapes  such  as  occur  for  fragments  of  bombs  or 

—  2/3 

shells  would  also  rotate  randomly  with  A  =  k  m  .  By  approximating 
real  fragments  from  a  3.7  inch  shell  with  rectangular  parallelepipeds, 
they  estimated  an  average  value  for  the  shape  factor,  k,  in  the  case 


1016.  J.  Taylor ,  A.  C.  79  (1940). 
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of  steel  fragments.  This  basic  approach  was  repeated  in  a  number  of 
later  British  reports  which  employed  a  greater  variety  of  fragment 
sizes,  shapes  and  initial  speeds. 

In  this  country,  Schwarzschild  and  Sachs102  used  high  speed  ' 
cinematography  of  bomb  detonations  to  determine  the  launch  speeds  and 
retardation  of  fragments  in  air.  However,  fragment  orientations  in 
flight  were  not  determined  because  of  the  great  distance  between  the 
camera  and  the  burst  point.  They  mentioned  that  spark  photography 
experiments  on  real  fragments  launched  from  a  gun  by  Charters  seemed  to 
confirm  the  random  rotation  hypothesis.  However,  a  report  published 
soon  after  by  Charters  and  co-workers10 3  was  more  modest  in  its 
conclusions.  They  presented  their  data  in  terms  of  drag  coefficients 
based  on  maximum  and  minimum  area  as  well  as  average  area  calculated 
by  approximating  fragment  shapes  with  ellipsoids,  as  suggested  by 
Schwarzschild  and  Sachs.  Although  rotation  was  observed,  its  exact 
nature  (random  orientation  or  not)  was  difficult  to  observe. 
Schwarzschild  and  Sachs  also  proposed  using  the  1.8  power  of  the  speed 
instead  of  the  2.0  power  on  the  right  side  of  Equation  (A-2)  above  for 
subsonic  speeds,  in  which  case 

Vo1/5  -  BR/5^  5  (A-4) 

The  question  of  whether  or  not  fragments  from  real  warheads 
actually  follow  the  random  rotation  hypothesis  at  least  over  distances 
of  tens  or  hundreds  of  meters  is  difficult  to  answer  satisfactorily 
in  an  experiment.  If  one  is  dealing  with  a  large  number  of  fragments 
of  various  sizes  and  shapes,  we  might  also  ask  whether  or  not  it  is 
worth  knowing  a  characteristic  shape  factor  for  each  fragment,  since 
we  are  interested  in  the  average  effects  of  the  whole  group  of  fragments. 
For  this  purpose  an  average  shape  factor  characteristic  of  the  group 
should  serve  just  as  well.  This  has  the  further  advantage  of  providing 
an  easy  way  to  estimate  the  striking  area  of  a  fragment  as  it  hits  the 
target,  a  factor  which  is  needed  if  we  are  to  use  penetration  equations 
based  on  relations  like  Equation  (56)  above.  The  random  orientation 
hypothesis  also  seems  to  have  appealed  to  a  number  of  people  because 
of  its  simplicity,  apart  from  the  question  of  describing  a  large 
collection  of  objects  statistically.  At  the  suggestion  of.  ^erse  and 
co-workers104  various  investigators  began  to  measure  average 

1112 M.  Schwarzschild  and  R.  G.  Sachs,  "Properties  of  Bomb  Fragments ,  " 

BRL  R347,  7  Apr  43. 

103V.  F.  Braun,  L.  H.  Thomas  and  A.  C.  Charters,  "Retardation  of 
Fragments,"  BRL  R425,  15  Nov  43.  (AD#PB22113) 

104M.  Morse,  W.  Transue  and  M.  Heins,  "The  Theory  of  the  Presentation 
Areas  of  Fragments  and  the  Icosahedron  Area  Gage,"  T.D.B.S.  44, 

20  Sep  44. 
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presentation  areas,  accepting  random  rotation.  For  example,  Simpson 
and  Bushkovitch* adopted  Morse's  suggestion  and  built  an  apparatus 
to  measure  twenty  uniformly  distributed  orientations  by  projecting 
fragment  shadows  in  directions  corresponding  to  the  twenty  faces  of 
an  icosahedron.  While  this  was  done  only  for  a  limited  number  of 
faces  and  relatively  small  number  of  selected  fragments,  it  was  a 
systematic  way  to  determine  experimentally  an  average  shape  factor, 
k.  These  authors  mentioned  parallel  efforts  in  England.  This  technique 
was  refined  in  succeeding  years10a”l°9  and  the  proportionality  between 
—  2/3 

A  and  m  has  been  well  established  experimentally  for  fragments 
recovered  from  a  variety  of  weapons  with  steel  cases.  Presumably  for 
other  metals  we  need  only  employ  a  ratio  of  densities  to  the  2/3  power 
as  indicated  by  Equation  (A-3)  to  calculate  a  form  factor  for  drag  or 
penetration  calculations,  provided  we  have  similar  geometrical  shapes. 

It  is  interesting  to  note  that  Equation  (A-3)  states  that  the 
average  presentation  area  of  a  rectangular  block  is  equal  to  one-fourth 
of  its  surface  area,  namely,  one-fourth  of  2(ab  +  be  +  ca) .  This  is 
also  clearly  true  for  a  sphere  and  other  simple  shapes.  The  result 
can  be  generalized  to  any  irregularly  shaped  convex  body  by  a  simple 
argument.  Choose  cartesian  and  spherical  coordinate  systems  with 
a  common  origin  such  that  the  radius  vector,  R,  is  parallel  to  the 
normal  vector,  cHCs, which  represents  an  element  of  surface  area.  The 
magnitude  of  the  projection  of  3/Cs  on  the  z  axis  is  dAs  cos  6  which 
is  also  the  size  of  the  shadow  this  element  of  area  makes  on  the  x-y 
plane.  If  we  vary  the  angle  <j> ,  holding  9  constant,  the  size  of  this 
shadow  stays  the  same.  If  we  vary  9,  holding  4>  constant,  the  shadow 
size  varies  from  dAs  to  zero  as  6  varies  from  0  to  n/2.  For  0  >  tt/2, 
the  element  of  area  is  in  the  half-space  below  the  x-y  plane,  is 
obseured  by  the  rest  of  the  body,  and  its  shadow  size  is  taken  to  be 
zero.  Consequently,  the  average  shadow  size  or  presentation  area, 
giving  equal  weight  to  all  aspects  is 


105M.  H.  Simpson  and  A.  V.  Bushkovitch ,  Fragment  Contour  Projector 
and  the  Presentation  Areas  of  Bomb  and  Shell  Fragments BEL  E501, 

8  Nov  44.  (AD# 715954) 

106/1.  V.  Bushkovitch ,  "Presentation  Areas  of  Shell  Fragments ," 

BEL  E536,  Apr  45. 

107 J.  E.  Shai) j  "A  Measurement  of  the  Drag  Coefficient  of  High  Velocity 
Fragments ,  "  BEL  E744,  Oct  50.  (AD# 80 1550) 

10aV.  E.  Porter,  J.  L.  Machamer  and  W.  0.  Ewing,  "Electro-Optic 
Icosahedron  Gage, "  BEL  E877 ,  Sep  53.  (AD#21197) 

109D.  J.  Dunn,  Jr.,  and  W.  E.  Porter,  "Air  Drag  Measurements  of 
Fragments,"  BEL  ME915,  Aug  55.  (AD#77240) 
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dAs  cos  0  sin  0d0d<}> 
sin  6  d  6  d  I 
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cos  0  sin  0d0 


4tt 


dA 

=  4~"(A-5) 


since  the  integral  over  0  from  ir/2  to  tt  is  zero  because  dA  is  obscured 

s 

by  the  rest  of  the  body  as  already  mentioned.  This  result  may  be 
repeated  for  each  element  of  surface  area,  although  generally  speaking, 
a  translation  of  the  origin  will  be  required  to  make  the  radius  vector 
parallel  to  the  element  of  surface  area  in  question.  Since  the  total 
surface  area  is  simply  the  sum  of  all  its  elements,  we  have  the  desired 
result. 


A  =  1/4  Ag  (A-6) . 

This  result  is  sometimes  attributed  to  Cauchy  although  no  reference 
has  been  found.  We  can  use  this  result  to  derive  Equation  (A-3)  by 
noting  that 

A  =  KV2/3  =  K(-)2/3  (A-7) 

s  P 

where  K  is  a  dimensionless  geometric  ratio  and  V  is  the  volume  of  a 
body  of  mass,  m,  and  density,  p.  If  we  put  Equation  (A-7)  into 
Equation  (A-6),  we  obtain  Equation  A-3)  in  the  form 

r  ,  K  ^  2/3  ,  2/3 

A  =  ( — jTTJ  m  =  k  m  .  (A-8) 

4p 

Here  the  shape  factor,  k,  depends  only  on  two  quantities.  One  is  the 
K-number  defined  by  Equation  (A-7)  as  the  ratio  of  suitable  powers  of 
two  extensive,  geometric  properties  of  the  body,  its  surface  area  and 
its  volume.  The  other  is  the  negative  two-thirds  power  of  an  inten¬ 
sive,  physical  property  of  the  body,  its  density.  In  particular, 
the  shape  factor,  k,  is  independent  of  mass.  However,  this  does  not 
eliminate  the  possibility  that  certain  k  values  may  occur  more 
frequently  for  certain  mass  groups  because  of  metal  processing  or 
some  other  factor. 

While  the  meaning  of  K  in  Equation  (A-7)  is  intuitively  evident 
and  its  numerical  value  is  easily  found  for  simple  geometrical  shapes, 
it  is  worth  noting  that  a  K-number  can  be  defined  for  any  convex 
body  no  matter  how  complicated  its  shape.  It  may  be  very  difficult 
to  calculate  this  number  for  an  oddly-shaped  body,  but  we  only  need 
to  define  it  in  order  to  derive  Equation  (A-8) .  The  use  of  an 
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icosahedron  gage  or  some  other  technique  will  enable  us  to  estimate  its 
value  to  the  accuracy  we  desire.  These  assertions  can  be  clarified 
by  the  following  whittling  experiment.  If  we  start  with  a  rectangular 
parallelepiped  with  edges  A,  B  and  C  we  may  whittle  out  any  convex 
shape  we  want  by  a  series  of  slices  each  of  which  produces  a  small 
pyramid  which  is  discarded,  leaving  behind  new  vertices  where  three 
edges  meet  in  addition  to  the  old  ones.  A  curved  surface  may  be 
produced  as  accurately  as  we  please  by  making  a  large  enough  number 
of  small  enough  cuts  (an  infinite  number  in  the  geometrical’ limit) . 
Suppose,  for  example,  we  start  by  cutting  off  a  corner  with  slant 
heights  a,  b  and  c  corresponding  to  the  edges  A,  B  and  C.  If  we  take 
edge  a  as  the  pyramid  height  and  the  triangle  with  b  and  c  as  sides 
for  the  base,  we  find  the  volume  JL  ab)  c  =  i.  abc.  The  surface 

area  of  the  pyramid  is  the  sum  of  the  triangular  areas  of  its  four  faces, 

namely,  i  ab,  *  be,  L  ca  and  /s (s-a) (s-b)s-c)  where  s  =  JL  (a+b+c) . 
222.  2 


This  sum  divided  by  (*  abc) 

6 


2/3 


is  the  K-number  of  the  pyramid.  The 

surface  area  of  the  remaining  mass  is  the  original  surface  area  minus 
that  part  of  the  pyramid  surface  area  which  was  exposed  before  the 
cut  plus  the  newly  exposed  surface  area.  The  volume  of  the  remaining 
mass  is  simply  the  original  volume  minus  the  pyramid  volume.  Thus, 
the  K-number  of  the  remaining  mass  is 


2 [AB+BC+CA]  -  y(ab+bc+ca)  +  /s(s-a) (s-b) (s-c) 


K  = 


,2/3 


[ABC  -  £  abc] 2/3 


(A-9) 


If  we  slice  off  a  new  pyramid  which  contains  one  of  the  newly  produced 
vertices,  the  calculation  will  be  complicated  by  the  fact  that, 
generally  speaking,  no  right  angles  will  be  present.  Still,  the 
calculation  of  a  K-number  for  the  remaining  mass  is  analytically 
possible.  The  same  is  true  if  we  had  started  with  a  polyhedron  other 
than  a  rectangular  parallelepiped.  As  the  number  of  cuts  increases 
the  analytical  formulas  for  the  K-number  of  the  remaining  mass  become 
increasingly  cumbersome  and  only  become  simple  again  as  some  simple 
geometric  shape  is  reached.  We  will  not  pursue  this  matter  here  since 
our  only  purpose  is  to  illustrate  the  fact  that  a  K-number  can  be 
defined  for  arbitrary  convex  shapes  like  those  which  most  often  occur 
in  natural  fragmentation. 

Finally,  it  is  interesting  to  note  the  dependence  of  v  on  the 
fragment  mass,  m,  in  Equation  (27)  above,  namely,  in 


v  =  v 

o 


e 


-aRm 


-1/3 


(A-10) 


81 


where  a  =  C  pa.k  varies  from  near  zero  for  very  high  altitudes  .  ,, 

(pa  58  0)  to  about  ,6(.001  gm/cm3)  (.38  cm2/gm2/3)  »  .2  x  10~3  gm  5 / 
cm  =  .02  gm  l'3/meter  near  sea  level.  Here  we  have  used  k  =  .38  cm2/ 
gm2/3  for  a  cubical  steel  fragment,  using  a=b=c  in  Equation  (A-3) . 

If  the  c  dimension  were  0.1  or  10  times  a=b,  then  k  would  only  be  about 
twice  as  large.  For  natural  fragmentation  c  is  usually  less  than  five 
times  a=b,  giving  a  k  variation  of  only  about  25%  from  the  cubical 
value  for  a  mass  variation  of  two  orders  of  magnitude.  Because  of 
this  it  is  common  to  choose  a  single  representative  value  of  k  to 
describe  a  large  collection  of  fragments  of  different  sizes  and  shapes. 
Let  us  do  this  so  that  a  is  the  same  for  all  fragments  in  Equation 
(A-10),  neglecting  the  even  smaller  variation  of  a  with  altitude  for 
a  given  burst.  In  addition,  let  us  assume  that  v0  is  the  same  for  all 
fragments,  an  approximation  which  is  nearly  true  if  the  charge  to 
mass  ratio  is  nearly  constant  over  the  length  of  a  shell  whose  lethal 
fragments  are  mostly  found  in  a  zone  near  the  equatorial  plane.  Then 
the  ratio  of  the  speeds  of  two  masses,  mj  and  m2,  at  distance  R  is 


R1  =  exp  [-  aR  (mj-1^3  -  nu'*^3)].  (A-ll) 

If  m2  =  1  gm  then  R^  can  be  found  as  a  function  of  m-^  and  R.  Table  I 
gives  representative  values  of  R1  for  m2  =  1  gm  and  a  =  .02  gm1/3  meter 
versus  m^  (gm)  and  R  (meters) .  If  we  compare  various  masses  with 
m2  =  0.1  or  10.0  gm  we  obtain  somewhat  different  values  of  R1.  Still 
Table  I  gives  us  a  feeling  for  the  variations  to  be  expected  over  the 
variable  ranges  of  interest  and  for  the  size  of  the  corrections  we 
might  expect  if  we  were  to  employ  different  values  of  k  or  a  for 
various  mass  groups.  Generally  speaking,  such  corrections  are  not 
justified  unless  the  other  errors  in  our  model  are  much  less  than  they 
usually  are. 

TABLE  I.  Speed  ratio  of  various  masses,  m  (gm) ,  at  various  distances 
from  the  burst  point,  R  (meters),  compared  to  a  1  gm  mass,  using 
a  =  .02  gm  1/'3/meter. 
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make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name : 


Telephone  Number: 
Organization  Address: 
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